我正在尝试在python中绘制SIR模型。但是，我找不到下面给出的绘制模型的好方法；
dS / dt = -aSI ...易感
dI / dt = aSI-bI…已感染
dR / dt = bI ...已恢复
您应注意，第三个方程式不影响前两个方程式。结果，它可以忽略。仅使用前两个方程式，开发一个Python脚本即可绘制模型的相空间。您应该得到这样的输出。
我真讨厌：
from pylab import *
g,L = 1,1
xvalues, yvalues = meshgrid(arange(-3, 3.1, 0.1), arange(-3, 3.1, 0.1))
vx = xvalues
vy = yvalues - xvalues
streamplot(xvalues, yvalues, vx, vy)
show()
但结果是：
