Function soi_radius 

pub fn soi_radius(planet_orbit_radius: Km, mu_planet: Mu, mu_sun: Mu) -> Km

Sphere of influence radius (km) using Hill’s approximation.

r_SOI = a_planet * (m_planet / m_sun)^(2/5)

where a_planet is the semi-major axis of the planet’s orbit.