Function vsop87::neptune

pub fn neptune(jde: f64) -> VSOP87Elements

Calculates VSOP87 solution for Neptune.

This function calculates the VSOP87 solution (heliocentric ecliptic orbital elements for the equinox J2000.0) for the planet Neptune. 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 Neptune. In this case, we calculate the orbit of Neptune in December 19th, 1499.

let vsop87_elts = vsop87::neptune(2268920.0);

assert!(vsop87_elts.a > 30.1963044187 && vsop87_elts.a < 30.1963044189);
assert!(vsop87_elts.l > 5.1088676118 && vsop87_elts.l < 5.1088676120);
assert!(vsop87_elts.k > 0.0091964091 && vsop87_elts.k < 0.0091964093);
assert!(vsop87_elts.h > 0.0031103619 && vsop87_elts.h < 0.0031103621);
assert!(vsop87_elts.q > -0.0102800265 && vsop87_elts.q < -0.0102800263);
assert!(vsop87_elts.p > 0.01148076 && vsop87_elts.p < 0.01148152);

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);