A two-loop sparse matrix numerical integration procedure for the solution of differential/algebraic equations: Application to multibody systems