SciPy庫中提供了多種求解常微分方程的方法，最常用的是使用odeint函數。下面是一個示例代碼，演示如何使用odeint函數求解常微分方程：
import numpy as np
from scipy.integrate import odeint
# 定義常微分方程的函數
def model(y, t):
dydt = -0.5 * y
return dydt
# 定義初始條件
y0 = 1
# 定義時間點
t = np.linspace(0, 10, 100)
# 求解常微分方程
y = odeint(model, y0, t)
# 打印結果
print(y)
在這個示例代碼中，首先定義了一個簡單的一階常微分方程dy/dt = -0.5 * y，然后定義了初始條件y0 = 1，接著定義了時間點t，最后使用odeint函數求解常微分方程并打印結果。
如果需要求解更復雜的常微分方程，可以根據實際情況修改model函數來定義更復雜的微分方程。SciPy庫中還提供了其他求解常微分方程的函數，例如ode和solve_ivp，可以根據具體需求選擇合適的方法來求解常微分方程。
