数値計算プログラム

input f$;
input h,n,m;
input x1,y1;

m1:=n/m;
x2:=x1;
y2:=y1;

for i:=1 to m1;
for j:=1 to m;
x:=x2;
y:=y2;
k1:=h*val(f$);
x:=x2+0.5*h;
y:=y2+0.5*k1;
k2:=h*val(f$);
y:=y2+0.5*k2;
k3:=h*val(f$);
x:=x2+h;
y:=y2+k3;
k4:=h*val(f$);
f:=(k1+2*(k2+k3)+k4)/6;
x2:=x2+h;
y2:=y2+f;
next j;
next i;
end