A new algorithm for solving differential/algebraic equations of multibody system dynamics

The second order Euler-Lagrange equations are transformed to a set of first orderdifferentiallalgebraic equations. which are then transformed to state equations by usinglocul parameterization. The corresponding discretization method is presented. and theresults can be used to implementation of various numerical integration methods. Anumerical example is presented finally