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 yn=yn+h*f(tn,yn)+h^2/2*df(tn,yn) end ye=1/(1+100*exp(-T)) println(h," ",abs(yn-ye)) end