The error accumulation of dynamics can be eliminated by using canonical equation and symplectic integral method so that the computational accuracy can be ensured effectively.