const Vector3f GravitationalForce::CalculateForce( const Body &a, const Body &b ) const { const Vector3f vDiff = b.Position - a.Position; const float numerator = Constant * ( a.Mass * b.Mass ); const float denominator = std::pow( vDiff.Magnitude(), 2 ); return vDiff.Normalised() * ( numerator / denominator ); }