% matlab script to illustrate stiffness % using simple flame propagation modelclose allclear all% define the right-hand side of the ODEF = inline(‘y^2 – y^3′,’t’,’y’);delta = 0.01; % not stiffticode23(F,[0 2/delta],delta);tocdisp(‘Press any key to continue.’)pause% try again, this time with stiffness …delta = 0.0001;ticode23(F,[0 2/delta],delta);tocdisp(‘Press any key to continue.’)pause% now use a stiff solver to solve the problemticode23s(F,[0 2/delta],delta);toc
，如同磁铁吸引四周的铁粉，热情也能吸引周围的人，改变周围的情况。
