Algebraic dynamical algorithm is used to solve numerically the artificial earth satellite motion equation for the first time, and the result is compared with Runge-Kutta algorithm and symplectic algorithm under the forth order. The result shows that algebraic dynamical algorithm has higher accuracy than the others in the long-term calculation of satellite orbit. We also discussed the quadrupole and octupole deformation of the Earth's impact on the satellite orbit.