数値計算プログラム

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

i1:=0;
i:=0;
w:=0
a:=0;
b:=0;

label li;
input w,n;
input x1,x2,x3;

if n/2=int(n/2) then a:=n-0.5:b:=n+0.5; else a:=n+0.5:b:=n-0.5;
for x:=x1 to x2 step x3;
f:=0;
for i=-a to a step 2;
i1:=i*p;
f:=fnf(x+i1/w)/i1^2+f;
next i;
for i:-b to b step 2;
i1:=i*p;
f:=f-fnf(x+i1/w)/i1^2;
next i;
next x;
input i2;

if i2=0 then l1;
end;