Construction and analysis of higher order variational integrators for dynamical systems with holonomic constraints
T. Wenger, S. Ober-Blöbaum, S. Leyendecker, Advances in Computational Mathematics 43(5) (2017) 1163–1195.