The differential-algebraic equation of the dynamic equations of constrained multibody system is a system of index three.It is a very difficult problem for solving. At present,there are no satisfactory numerical methods yet,we construct the symplectic composition Runger-Kutta methods for the differential-algebraic equations,and use mainly it for simulation in the canonical form of the equations.The results arebetter than Baumgarte's methods.