We introduce energy-preserving integrators for nonholonomic mechanical systems. We will see that the nonholonomic dynamics is completely determined by a triple ({{\mathcal {D}}}^*, \varPi , \mathcal {H}), where {{\mathcal ...