Аннотация:The high-precise analytical theory of a nonspherical planet satellite has been developed. All perturbations proportional up to a product Jlm·Jpq·Jrs·Jkn (where Jij is an arbitrary harmonic of planet potential expansion) are calculated. The influence of the high-order perturbations on the motion of some satellite in long-term intervals has been investigated. The results of comparison of the numeric and analytical integrations of the satellite motion equations in two years' interval are as follows: for the ETALON-1 satellite the r.m.s. error is 5.6 cm; for the Phobos moon the r.m.s. error is 2.5 cm. The theory is intended to be used for processing precise laser measurements of the Earth passive geodynamic satellites in long-term intervals.