python求微分方程组的数值解

要求解微分方程组的数值解，可以使用Python中的SciPy库中的odeint函数。下面是一个求解微分方程组的例子：

假设有以下微分方程组：

y1' = y2

y2' = -y1

并且初始条件为y1(0) = 1, y2(0) = 0。

可以使用以下代码来求解：

from scipy.integrate import odeint
import numpy as np

def model(y, t):
    y1, y2 = y
    dydt = [y2, -y1]
    return dydt

# 初始条件
y0 = [1, 0]

# 创建时间点
t = np.linspace(0, 10, 101)

# 求解微分方程组
y = odeint(model, y0, t)

# 输出结果
print(y)

在上面的代码中，model函数定义了微分方程组的形式，y0定义了初始条件，t定义了时间点，odeint函数求解微分方程组，并返回数值解。最后打印结果即可。

注意，以上仅为示例，实际应用时需要根据实际情况定义相应的微分方程组和初始条件。