Function vsop87::mars

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

Calculates VSOP87 solution for Mars

This function calculates the VSOP87 solution (heliocentric ecliptic orbital elements for the equinox J2000.0) for the planet Mars. 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::mars(2159345.0);

assert!(a > 1.5237578877 && a < 1.5237578879);
assert!(l > 4.0669853278 && l < 4.0669853280);
assert!(k > 0.0821906316 && k < 0.0821906318);
assert!(h > -0.0427917583 && h < -0.0427917581);
assert!(q > 0.0103081045 && q < 0.0103081047);
assert!(p > 0.01313608 && p < 0.01313684);