Abstract
It is shown that linear time-dependent invariants for arbitrary multidimensional quadratic systems can be obtained from the Lagrangian and Hamiltonian formulation procedures by considering a variation of coordinates and momenta that follows the classical trajectory and defines a Noetherian symmetry transformation.