>> f=@(t,y) y^2-y^3;
>> ode23(f, [0, 2e4], 1e-4)
>> ode23s(f, [0, 2e4], 1e-4)