The 5th-order approximation for the y-component of
used in this paper is
where
.
Writing x = x(0), y = y(0) for the initial conditions, we have
,
and
where
and
.
The x-component of
is then just
to 5th-order in h.