  This paper presents a method to construct variational integrators for
time-dependent lagrangian systems. The resulting algorithms are symplectic,
preserve the momentum map associated with a Lie group of symmetries and also
describe the energy variation.
