python求解微分方程组

Python可以使用SciPy库的odeint函数来求解微分方程组。

首先，需要导入相关的库：

from scipy.integrate import odeint
import numpy as np

然后，定义微分方程组的函数，例如：

def myODE(y, t):
    dydt = np.zeros((3,))
    dydt[0] = -0.1*y[0] + 0.2*y[1]*y[2]
    dydt[1] = 0.1*y[0] - 0.2*y[1]*y[2] - 0.3*y[1]
    dydt[2] = 0.3*y[1] - 0.2*y[1]*y[2]
    return dydt

其中，y是微分方程组的解向量，t是时间。

接下来，定义初始条件和时间步长：

y0 = [1.0, 2.0, 3.0]
t = np.linspace(0, 10, 101)

最后，调用odeint函数求解微分方程组：

sol = odeint(myODE, y0, t)

sol即为微分方程组的数值解。可以使用matplotlib库将结果可视化：

import matplotlib.pyplot as plt

plt.plot(t, sol[:, 0], label='y1')
plt.plot(t, sol[:, 1], label='y2')
plt.plot(t, sol[:, 2], label='y3')
plt.legend()
plt.show()

这样就可以得到微分方程组的数值解，并可视化结果。