R. J. Y. McLeod & J. M. Sanz-Serna, Geometrically derived difference formulae for the numerical integration of trajectory problems, IMA J. Numer. Anal. 2(1982), 357-370.