% flame propagation

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

