y0=1/101 f(t,y)=y*(1-y) df(t,y)=(1-2*y)*f(t,y) for j=1:12 m=2^j T=10 t0=0 h=(T-t0)/m yn=y0 ynm1=1/(1+100*exp(h)) for n=1:m tn=t0+h*n yp=ynm1+2*h*f(tn,yn) ynp1=yn+h/12*(-f(tn-h,ynm1)+ 8*f(tn,yn)+5*f(tn+h,yp)) ynm1=yn yn=ynp1 end ye=1/(1+100*exp(-T)) println(h," ",abs(yn-ye)) end