SymPySciPy：求解具有不同变量的常微分方程组

SymPy / SciPy：求解具有不同变量的常微分方程组

如果要使用读取文件的相同脚本来解决系统（因此

systemOfEquations

可以用作全局变量），并且使用的 唯一
变量

systemOfEquations

是

x

，

y

并且可能使用，则

t

可以

dX_dt

在同一文件中定义如下：

def dX_dt(X, t):    vals = dict(x=X[0], y=X[1], t=t)    return [eq.evalf(subs=vals) for eq in systemOfEquations]

dX_dt

可用于中

odeint

。在下面的ipython会话中，我已经运行了创建

systemOfEquations

和定义脚本

dX_dt

：

In [31]: odeint(dX_dt, [1,2], np.linspace(0, 1, 5))Out[31]: array([[ 1.        ,  2.        ],       [ 1.00947534,  1.90904183],       [ 1.01905178,  1.82223595],       [ 1.02872997,  1.73939226],       [ 1.03851059,  1.66032942]]