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 system of ordinary differential equations is integrated in time with the energy preserving average vector field (AVF) method.