A variable-step double-integration multi-step integrator is derived using divided differences. The derivation is based upon the derivation of Shampine-Gordon, a single-integration method. Variable-step integrators are useful for propagating elliptical orbits, because larger steps can be taken near apogee. As a double-integration method, the integrator performs only one function evaluation per s...