where the index j runs from 1 to M. For brevity, these functions g i are grouped into an M -dimensional vector g below. The task is to solve the combined set of differential-algebraic (DAE) equations, instead of just the ordinary differential equations (ODE) of Newton's second law.