help ode45;
pause
delta = 0.01;
F = inline('y^2 - y^3','t','y');
opts = odeset('RelTol',1.e-4);
ode45(F,[0 2/delta],delta,opts);
pause
figure
delta = 0.0001;
ode45(F,[0 2/delta],delta,opts);
opts = odeset('RelTol',1.e-6);
pause
figure
ode45(F,[0 2/delta],delta,opts);
pause
help ode23s
pause
figure
ode23s(F,[0 2/delta],delta,opts);
