Discrete Conservative Vector Fields Induced by the Trapezoidal Method

The use of symmetric schemes has revealed interesting stability properties for the long time simulation of conservative, and in particular Hamiltonian systems. By using one of the simplest symmetric formulae, namely the trapezoidal method, we show that, under certain circumstances, one can attach to the discrete system a discrete energy function and therefore obtain a discrete conservative dynamical system.