Advances in the Numerical Solution of Constrained Differential Equations