delta = 0.01;
F = inline('y^2 - y^3','t','y');
opts = odeset('RelTol',1.e-4);
subplot(2,1,1)
ode45(F,[0 2/delta],delta,opts);
subplot(2,1,2)
ode45(F,[0 2/delta],delta,opts);
axis([0 2/delta .9998 1.0002])

figure
delta = 0.0001;
subplot(2,1,1)
ode45(F,[0 2/delta],delta,opts);
subplot(2,1,2)
ode45(F,[0 2/delta],delta,opts);
axis([9800 11000 .9998 1.0002])

figure
delta = 0.0001;
subplot(2,1,1)
ode23s(F,[0 2/delta],delta,opts);
subplot(2,1,2)
ode23s(F,[0 2/delta],delta,opts);
axis([9800 11000 .9998 1.0002])