前回、ロジスティック方程式について、数値解法を自作の関数を作成して解いた。 akh1116.hatenablog.com 今回は、scipyライブラリのintegrateモジュールを用いてコーディングをし、解いてみる。 自作で関数を作成するより、モジュールに渡す引数に注意すれば…

今回は、Pythonによる常微分方程式を解く際によく使われるscipyライブラリのメモである。 Scipyによる常微分方程式 Pythonを用いて常微分方程式を解く際は、ライブラリを用いるのが最も手軽で、かつ高精度な結果が得られる。 常微分方程式を解くライブラリと…

常微分方程式の数値計算例題-ロジスティック方程式 ここでは、常微分方程式の数値計算の例題として、ロジスティック方程式を扱ってみる。 ロジスティック方程式 ロジスティック方程式は、生物の個体数の増加を表すモデルであり、以下のような形の方程式とし…

本稿から具体的な数値計算の手法についての話に入る。 まずは、古典力学の基本方程式であるNewtonの運動方程式を数値的に解く手法から。 古典力学と常微分方程式 古典力学は、基本的には質点の運動はNewtonの運動方程式で記述される。 Newtonの運動方程は、…