Store a member variable for density in each planet, and if the density is above a certian range, apply a first order correction (gravitational Lorentz transformation) to create some precession. Like the “anomalous perihelion precession of mercury”. If anybody has questions about the texhnical nature of what I’m talking about, I’d be super happy to talk about it