KeplerOrbit

Routines for working with near-Keplerian orbits. Convert between cartesian and Keplerian orbital elements, calculate libration widths, precession rates...etc. Routines are vectorized so that you can easily work with large numbers of bodies.