Function vsop87::earth_moon

pub fn earth_moon(jde: f64) -> (f64, f64, f64, f64, f64, f64)

Calculates VSOP87 solution for Earth - Moon barycenter

This function calculates the VSOP87 solution (heliocentric ecliptic orbital elements for the equinox J2000.0) for the Earth - Moon barycenter. The parameter needed is the Julian Day Efemeris (JDE) for the given date. It returns, in order, a tuple with the values a, l, k, h, q, p of the VSOP87 solution.

Examples

let (a, l, k, h, q, p) = vsop87::earth_moon(2122820.0);

assert!(a > 1.0000134925 && a < 1.0000134927);
assert!(l > 1.8519621672 && l < 1.8519621674);
assert!(k > -0.0029638176 && k < -0.0029638174);
assert!(h > 0.0168402193 && h < 0.0168402195);
assert!(q > 0.0010301900 && q < 0.0010301902);
assert!(p > -0.00005346 && p < -0.00005270);