stiffex1.m

来自「matlab宝典,电子工业出版社,包含该书的源码2-4章」· M 代码 · 共 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 ];