We discuss systematic extensions of the standard (Störmer–Verlet) method for integrating the differential equations of Hamiltonian mechanics. Our extensions preserve the symplectic geometry exactly, as well as all Nöether conservation laws caused by joint symmetries of the kinetic and potential energies (like angular momentum in rotation invariant systems). These extensions increase the accuracy of the integrator, which for the Störmer–Verlet method is of order τ2 for a timestep of length τ, to higher-orders in τ. The schemes presented have, in contrast to most previous proposals, all intermediate timesteps real and positive. The schemes increase the relative accuracy to order τN (for N = 4, 6 and 8) for a large class of Hamiltonian systems.