% Moler, NCM, odes.pdf, p. 19

F = @(t,y) y^2 - y^3;

%tic;
delta = 0.01;  % 0.0001
opts = odeset('RelTol',1.e-4);
ode45(F,[0 2/delta],delta,opts); grid
%[t,y] = ode45(F,[0 2/delta],delta,opts);
%[t,y] = ode23s(F,[0 2/delta],delta,opts);
%[t,y] = ode15s(F,[0 2/delta],delta,opts);
%toc 
%plot(t,y,'o-'), grid
