Y=ill(t,x)實際上是
y'=ill(t,x),
y必須是列向量。
y(1)
=
dy1/dt,
y(2)=dy2/dt
所以這條線是
dy1/dt
a*
y1*y2-b*y1,
dy2/dt=
-a*y1*y2
就是sir模型方程(y1=i,
y2=s)