HamaguRePythonで数値解法を試す先日、「パソコンで見る天体の動き」という本を買いました。名前の通り、天体の軌道をパソコンで表示するための数式やサンプルコードなどが載っています。内容はとても素晴らしいのですが、サンプルコードがBASICやFortranで書かれていたので、理解を深めるためにいくつかの解法をPythonで書き直してみました。数値解法の詳しい説明を知りたい方は、是非「パソコンで見る天体の動き」をご覧ください。以下の微分方程式をオイラー法とルンゲ・クッタ法、フェールベルグ法を使って解いてみます。関数の傾き」で増加量が求まる。