空気抵抗中の落下

速度に比例する抵抗を受けながらの落下について、Pythonを使って、Euler法、odeint、Runge-Kutta法の3種類の計算をやってみた。グラフは比較用の解析解を加えた4つをプロットしているけど、一目で分かる違いなんてない。ちなみにグラフは上が変位、下が速度である。まぁ、終端速度がありますね、というくらいの話である。

しばらくプログラミングは何もしていなかったので、タイピングがどこかぎこちなくて、打ち間違いが滅茶苦茶多かったが、ようやくマシになってきた。