Function vsop87::jupiter

pub fn jupiter(jde: f64) -> VSOP87Elements

Calculates VSOP87 solution for Jupiter.

This function calculates the VSOP87 solution (heliocentric ecliptic orbital elements for the equinox J2000.0) for the planet Jupiter. The parameter needed is the Julian Day (JD) for the given date. It returns the VSOP87Elements of the VSOP87 solution.

Example

Given a date in JD, we can get the orbit of the planet Jupiter. In this case, we calculate the orbit of Jupiter in December 30th, 1799.

let vsop87_elts = vsop87::jupiter(2378495.0);

assert!(vsop87_elts.a > 5.2027276672 && vsop87_elts.a < 5.2027276674);
assert!(vsop87_elts.l > 1.4820596291 && vsop87_elts.l < 1.4820596293);
assert!(vsop87_elts.k > 0.0464780412 && vsop87_elts.k < 0.0464780414);
assert!(vsop87_elts.h > 0.0116460263 && vsop87_elts.h < 0.0116460265);
assert!(vsop87_elts.q > -0.0019921307 && vsop87_elts.q < -0.0019921305);
assert!(vsop87_elts.p > 0.01123447 && vsop87_elts.p < 0.01123523);

It can then be converted into keplerian elements:

use vsop87::{KeplerianElements, VSOP87Elements};

let k_elements: KeplerianElements = vsop87_elts.into();
let convert_back = VSOP87Elements::from(k_elements);