数値計算プログラム

def fnf(x):=sin(x) */exanple/*

input a,b,n,e
s:=0
repeat
d:=b-a
x:=d/2+a
n1:=1
t:=abs(fnf(x)*d*2/3
s:=t+(abs(fnf(a))+abs(fnf(b)))*d/6
i:=2
repeat
s1:s
s:=(s-t/2)/2
n1:=n1*2)
c:=d/n1/2+a
t:=abs(fnf(c))
c1:=c
c2:=d/n1
j:=1
repeat
c1:=c1+c2
t:=t+abs(fnf(c1))
until j t:=t*d*2/(n1*3)
s:=s+t
if abs(s-b1) until i i:=i+1
label1
until ((a>b)+(n<=1))>0
end