Variational methods are a class of symplectic-momentum integrators for ODEs. Using
these schemes, it is shown that the classical Newmark algorithm is structure preserving in a
non-obvious way, thus explaining the observed numerical behavior. Modifications to variational
methods to include forcing and dissipation are also proposed, extending the advantages
of structure preserving integrators to non-conservative systems
