This IDL function solves Kepler's equation for an elliptic orbit.