数値計算プログラム

input f$;
input h,n;

dim x(n),y(n);
input x(1),y(1);

for i=2 to n;
x:=x(i-1);
y:=y(i-1);
y(i):=y(i-1)+h*val(f$);
x(i):=x(i-1)+h;
next i;

label l1;
input x1,x2,x3;

if x1x(n) then l1;
h1:=x3/h;
h2:=(x1-x(1))/h+1:
h3:=(x2-x(1))/h+1;

for i=h2 to h3 step h1;
print x(1),y(i);
next i;

end;