Abstract: The analytical formulas for calculating satellite orbit perturbation, which is caused by zonal spherical harmonics in the earth's gravitational field, are given in this paper. The number of the term of harmonics can be chosen as needed. The formulas are valid for orbits with any eccentricities.