These IDL routines convert Keplerian orbit elements into Cartesian coordinates and velocities, and back