The implicit symplectic algorithm of Hamilton multibody system with topological tree configuration is studied. The canonical equations of multibody system and Jacobian matrix of RHS of the equations are obtained in the form of matrix.The paper presents algorithm of Jacobian matrix to raise computational efficiency. The implicit symplectic Runge Kutta algorithm is used in solving the canonical equations of multibody system. Numerical results show that the algorithm has higher computational efficiency and can keep computation stable for long time simulation.