admin

3.2. 二階の常微分方程式¶　二階のODEを解くには、二元連立一階微分方程式に置き換えてやればいい。具体例でみてみよう。 (例題) $\ddot{x} + \dot{x} + x = 0$を、初期値$x(0)=0$、$\dot{x}(0)=1$で解き、$0 [...]

3. 常微分方程式(Ordinary Differential Equation, ODE)¶　常微分方程式(Ordinary Differential Equation, ODE)は物理のあらゆる場面で現れる。この微分方程式を数値計算で解いてみよう。 3.1. odeint¶　微分方程式(ODE)のアルゴリズムには、オイラー法、Runge-Kutta法などがある。一度くらいはこれらの解法の仕組みを知っておく必要があるが、それは後回しにして、SciPyに用意されているODEを解くfunctionであるodeintを使って、実際に計算してみることから始めよう。なお、ODEを解くfunctionのことをODEソルバーと言う。 [...]