% matlab script to illustrate stiffness % using simple flame propagation model close all clear all % define the right-hand side of the ODE F = inline('y^2 - y^3','t','y'); delta = 0.01; % not stiff tic ode23(F,[0 2/delta],delta); toc disp('Press any key to continue.') pause % try again, this time with stiffness ... delta = 0.0001; tic ode23(F,[0 2/delta],delta); toc disp('Press any key to continue.') pause % now use a stiff solver to solve the problem tic ode23s(F,[0 2/delta],delta); toc