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.Comment: 10 pages, 7 figures. Accepted by MNRAS, resubmitted because of a typo
in the title, added author affiliation