Line integral solution of Hamiltonian systems with holonomic constraints

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.