Function phase_angle 

pub fn phase_angle(planet1: Planet, planet2: Planet, jd: f64) -> Radians

Compute the angular separation between two planets at a given Julian Date.

Returns the signed angle from planet1 to planet2 in the direction of orbital motion (counter-clockwise), normalized to (-π, π].