  In this paper, we propose a second-order energy-conserving approximation
procedure for Hamiltonian systems with holonomic constraints. The derivation of
the procedure relies on the use of the so-called line integral framework. We
provide numerical experiments to illustrate theoretical findings.
