k:=0;
j:=0;
input f$,h,e,n,l;
dim x(n+1),y(n+1),s(l+1),t(l+1);
input x(1),y(1);
for i:=1 to n;
x:=x(i);
y:=y(i);
t(i):=y+h*val(f$);
a:=2;
for j:=1 to l;
a1:=a;
x:=x(i);
y:=y(i);
for k:=1 to a;
f:=val(f$);
x:=x+h/a1;
y:=y+h/a1*f;
next k;
s(1):=y;
f1:=1;
for k:=1 to j;
f1:=f1*2;
s(k+1):=s(k)+(s(k)-t(k))/(f1-1);
if abs(s(k+1)-s(k))
j1:=j+1;
for k:=1 to j1;
t(k):=s(k);
next k;
a:=a+2;
next j;
goto l2;
label l1;
x(i+1):=x(i)+h;
y(i+1):=s(k+1);
next i
label l2;
input h1;
for i:=1 to n step h1/h;
print x(i),y(i);
next i;
end;