  In this paper, we define arbitrarily high-order energy-conserving methods for
Hamiltonian systems with quadratic holonomic constraints. The derivation of the
methods is made within the so-called line integral framework. Numerical tests
to illustrate the theoretical findings are presented.
