Python は一次常微分方程式を解きます

1. Python を使用して一次常微分方程式を解く

1. 微分方程式を解くには、pycharm にインストールできるscipyライブラリを使用すると同時に、numpy ライブラリと matplotlib をインポートする必要があります。

2. scipy.integrate.odeint() を使用して解決します。これは通常、3 つのパラメーターを使用し、他のパラメーターをデフォルトにすることを指します。

func: 導関数 f(y,t)、つまり関数形式で表された t における y の導関数

y0: 初期条件 y0、配列形式で表現されます。

t: 関数値に対応する一連の時点を解きます。シーケンスの最初の要素は、初期条件 y0 に対応する初期時間 t0 です。時系列は単調増加または単調減少である必要があり、値の繰り返しが許可されます。

scipy.integrate.odeint(func, y0, t, args=(), Dfun=None,col_deriv=0, full_output=0, ml=None, mu=None, rtol=None, atol=None, tcrit=None, h0 =0.0, hmax=0.0, hmin=0.0, ixpr=0, mxstep=0, mxhnil=0, mxordn=12, mxords=5, printmessg=0, tfirst=False) 3.
举例

 この機能については

#导入库函数
from scipy.integrate import odeint
import numpy as np
import matplotlib.pyplot as plt 
#定义函数的导数
def dy_dt(y,t):
    return np.sin(t**2)
#给定初始值和时间范围
y0=[1]
t = np.arange(-10,10,0.01)
#使用odeint()方法：
y=odeint(dy_dt,y0,t)
#绘图
plt.plot(t, y)
plt.title("picture")
plt.show()