Stiff な方程式なのか、うまく解けなかったときのメモです。 連続槽型反応器の濃度変化をグラフ化したく、常微分方程式を数値的に解くことを考えました。 python scipy にあるodeint でまずは解いてみました。 pythonのodeintに渡す関数の引数の中の並びが標準で ...

scipy.integrate モジュールの中の odeint 関数や solve_ivp 関数は、ルンゲクッタ法を利用して常微分方程式を解くのに特に適しています。 solve_ivp はより新しく、柔軟なAPIを提供しており、様々な種類のルンゲクッタ法（例えばRK45, RK23など）を選択できます。

Open DOS Prompt on Windows or Terminal on Linux/Mac OSX Type "python -m pykpp model.kpp" where model.kpp is a valid KPP input file. Example below assumes Installation and that mcm_ch4.kpp has been ...

Thus far we have seen several examples defined for real values. odeint can handle complex state types, hence ODEs which are defined on complex odeint can also be used to solve ordinary differential ...