This paper presents a method to construct variational integrators for timedependent lagrangian systems. The resulting algorithms are symplectic, preserve the momentum map associated with a Lie group of symmetries and also describe the energy variation.Research in part supported by grants DGICYT (Spain), PB97-1257 and PGC2000-2191-E.Peer reviewed