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

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