我正在评估一组系数随时间变化的常微分方程。coefficients['gamma']* I return dSdt, dIdt, dRdt 特别是,我在一个预先计算的数组中有equations over the time grid, t.
ret = odeint(deriv,y0,t,args=(N,coefficients)) 当我运行odeint时,它也被评估超过max(t)<e
我有一个常微分方程系统,我试图在其中包含一个“误差”项,这样它就变成了一个随机常微分方程系统。 为了解决python中的ODE系统,我通常使用scipy的odeint。一个来自Scipy Cookbook的例子,涉及著名的僵尸启示录: # zombie apocalypse modelingimport matplotlib.pyplotDead Pop.; No New Bi