数値計算プログラム

def fnf(x):=sin(x);

dim m(27);
data 1,0,2,-0.577350269,0.577350269,3,-707106781;
data 0,0.707106781,4,-0.794654472,-0.187592474,0.187592474;
data 0.794654472,5,-0.832497487,-0.374541409,0;
data 0.374541409,0.832497487,6,-0.866246878,-0.422518654;
data -0.266635402,0.266635402,0.422518654,0.866246818;

for i:=1 to 27;
read m(i);
next i;

label l1;
input a,b,e,n,l;
if a>=b l1;
n1:=0;
for i:=1 to n;
n1:=n1+i;
next i;
dim n(n);
for i:=1 to n;
n(i):=m(n1+i);
next i;
a3:=(b-a)*2;
for i:=1 to l;
a5:=a;
a3:=a3*0.5;
s:=0;
for j:=1 to 2~(i-1);
a4:=a5;
a5:=a5+a3;
for k:=1 to n;
x:=(a4+a5+a3*n(k))*0.5;
s:=s+abs(fnf(x));
next k;
next j;
s:=s*a3/n;
if i:=1 then l2;
if abs(s-s1) label l2;
print i,d;
s1:=s;
next i;
i:=i-1;
labl l3;
print s,i;
end;