Title: Formulation and numerical solution of the equations of constrained mechanical motion