位置と速度を分離して、万有引力は個々の更新で処理しない。

万有引力は距離関係ないけどまぁ有効距離を設定したとして、一つ一つの星にかかる引力を合計して
その星の質量から加速度と速度を計算する。
で、単位時間（多分1フレーム）あたりの移動量がでるから、それを位置に加算するだけ。