pub fn calculate_position( true_anomaly: f32, heliocentric_distance: f32, argument_of_periapsis: f32, mean_anomaly: f32, ) -> (f32, f32, f32)