Надо рассчитать движение спутника с круговой орбитой высотой 900 км. Если численно интегрировать его движение с учётом 2-й зональной гармоники геопотенциала, то модуль расстояния до спутника от центра Земли должен быть постоянным в течение оборота, поскольку вековые изменения из-за сплюснутости Земли не влияют на геометрию орбиты, а только на угловые параметры орбиты. А у меня получается изменение модуля между началом и концом оборота на +-15-20 километров. Почему?