ましろのログスマートフォン用の表示で見るましろのログこの広告は、90日以上更新していないブログに表示しています。参考:Pythonで微分方程式を解きます。アルゴリズムはRungeKutta法を使います。忘れた頃に定期的に必要になってくるのでここに関数化して書いておきます。下にコードがあるのでコピペして使ってください。numpyをimportしないと使えません。ような一階微分方程式をyについて解くことができます。この微分方程式に従ってRungeKutta法でtに関して発展させて行きます。pythonは関数がオブジェクトとして引数に使えます。戻り値t発展の経過がlist型で返ってきます。ndarrayではありません。厳密解コード例RunkeKuttaを関数かしたおかげで1行でRungeKuttaができてます。これからはこれをコピペして物理の時間発展などの微分方程式を解いていこう。参考:JupyterをGistで共有する。