pub fn eccentric_anomaly(e: f64, period: f64, t_p: f64) -> f64
Calculates the eccentric anomaly in degrees