The equation of satellite orbital motion was derived in the conventional terrestrial coordinate system.

First, now in all common numerical methods for the satellite dynamics equations the most important problem that needs to be improved is stability problem.

The differential equations that describe the artificial satellite motion are very complicated, and it is impossible to get the analytic express of solution.