% 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