c++ - 中心に引力のある球

Question

球の中心（惑星のように）に重力を適用する最も簡単な方法は何だろうと思います。ODEが使用するのに最適なライブラリになると思いますか？しかし、球の中心に重力を適用するためにそれを実装する方法がわかりません。

Answer 1

あなたに必要なのは ：

a = F / m = -GM r / r ^ 3

ここで、rmは質量球からの質量の半径方向ベクトルでありM、Gは重力定数です。

これを統合するために、必ずしも多くのライブラリは必要ありません。各タイムステップの加速度を計算するだけです。唯一のポイントは、各ステップで最初におおよその新しい位置を計算し、次にタイムステップの中点での位置の加速度を計算し、それを使用して決定する場合、はるかに良い結果（たとえば、閉じた軌道！）が得られることです。新しい位置と速度。