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 for n=1:m tn=t0+h*n yp=yn+h*f(tn,yn) yn=yn+h/2*(f(tn,yn)+f(tn+h,yp)) end ye=1/(1+100*exp(-T)) println(h," ",abs(yn-ye)) end