無次元化した形で設定する。すなわち として計算する(G は万有引力定数、M は質量の単位、T は時刻の単位、L は長さの単位)。
まず 2 体だけで考える。
posA -- A の位置。 posB -- B の位置。 massA -- A の質量。 massB -- B の質量。 vectAB = posB - posA -- A から B に向かうベクトルを求める。 lengthABsq = vectAB:lengthsq() -- そのベクトルの距離の二乗を求める。 normAB = vectAB:normalize() -- 単位ベクトルも求める。 forceA = normAB:mult(massA * massB/lengthABsq) -- A にかかる force を求める。 forceB = -forceA -- B にかかる force を求める。
参考:
pp.130-131, p.184