Python : 常微分方程式

odeint, BDFの利用について
Stiff な方程式なのか、うまく解けなかったときのメモです。

連続槽型反応器の濃度変化をグラフ化したく、常微分方程式を数値的に解くことを考えました。

問題反応器サイズ　1L のCSTRに、反応速度式 $${-r_A = k C_A^{1.558}}$$となる原料を 100mmol/Lの濃度で、流量 10, 3, 1.2, 0.5 L/hrの各流量変化させ