AdamsBashforth:=proc(k)
local p,i,dt,t;
  p:=interp([seq(i*dt,i=0..k)],[seq(f[i],i=0..k)],t);
  y[k+1]-y[k]=factor(int(p,t=k*dt..(k+1)*dt));     
end;