The shape and perigee of artificial satellite orbits shall continually variate due to caused by the atmospherical drag and the earth gravity.

The dynamical model included the effects of the earth non-spherical shape, atmospheric drag, solar radiation pressure as well as lunar and solar perturbations.

A new method of fast computing the spherical harmonic perturbation on satellites and the recurrence relations of the coefficients of the Earth's gravity are given.