We propose a novel algorithmic method for constructing invariant variational schemes of systems of ordinary differential equations that are the Euler–Lagrange equations of a variational principle. The method is based on ...