数値計算プログラム

def fnf(x):=sin(x) <---exsample-->

dim m(27);
data 1,0,2,-0.577350269,0.577350269,3,-0.707106781;
data 0,0.707106781,4,-0.794654472,-0.187592474,0.187592474;
data 0.794654472,5,-0.832497487,-0.374541109,0;
data 0.374541409,0.832497487,6,-0.866246818,-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;
if a>=b then label l1;
a1:=a+b;
a2:=b-a;
i1:=1;
x:=a1*0.5;
s1:=2*abs(fnf(x));
for i:=2 to 6;
i1:=i1+i;
s:=0;
print i-1,s1;
for j:=1 to i;
j1:=i1+j;
x:=(a1+a2*m(j1))*0.5;
s:=s+abs(fnf(x));
next j;
s:=s*a2/m(i1);
if abs(s-s1) s1:=s;
next i;
i:=i+1;
label l2;
print "resalt=";s;spc$(3);"step=";i;
end