Before revision:  The resulting system of ordinary differential equations is integrated in time by the energy preserving average vector field (AVF) method.
After revision:  The resulting skew-gradient system is integrated in time by the energy preserving average vector field (AVF) method.