  We derive and present a fast and accurate solution of the initial value
problem for Keplerian motion in universal variables that does not use the
Stumpff series. We find that it performs better than methods based on the
Stumpff series.
