pub fn radius(a: f64, e: f64, period: f64, t_p: f64) -> f64
Calculates the total radius for a twobody system