数値計算プログラム

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

input a,b,n,e
s:=0
if ((a>b)+(n<=1))>0 then l1
d:=b-a
x:=d/3+a
y:=d*2/3+a
n1:=1
t:=(abs(fnf(x))+abs(fnf(b)))+d/8
for i:=2 to n
s1:=s
s:=(s-t/3)/3
n1:=n1*3
c5:=d/n1/3*2+a
c4:=d/n1/3+a
t:=abs(fnf(c4))+abs(fnf(c5))
c6:=c4
c7:=c5
c2:=d/n1
for j:=1 to n1-1
c6:=c6+c2
c7:=c7+c2
t:=t+abs(fnf(c6))+abs(fnf(c7))
next t
t:=t*d*3/(n1*8)
s:=s+t
if abs(s-s1) next i
i:=i-1
label l2
end
label l1
end