stiffex1.m

matlab宝典,电子工业出版社,包含该书的源码2-4章

function dydt = stiffex1(t,y)
dydt = [ (-0.04*y(1) + 1e4*y(2)*y(3))
         (0.04*y(1) - 1e4*y(2)*y(3) - 3e7*y(2)^2)
         3e7*y(2)^2 ];