В учебнике
Л.Д.Ландау Е.М.Лифшиц "Механика" п.14 "Движение в центральном поле" рассмотрен вывод зависимости (в полярных координатах) угла от длины радиус-вектора, с началом в центре тяжести системы и концом во втором теле:
![$\varphi=\int\frac{M/r^2dr}{\sqrt{2\cdot m[E-U(r)]-\frac{M^2}{r^2}}}$ $\varphi=\int\frac{M/r^2dr}{\sqrt{2\cdot m[E-U(r)]-\frac{M^2}{r^2}}}$](https://dxdy-04.korotkov.co.uk/f/3/9/3/393d54a6d13219cfd594d578a3072ae582.png)
(1)
где

- момент импульса,

- полная энергия системы,

- функция потенциальной энергии тела от длины радиус-вектора.
Далее показывается, что пусть даже пространство и финитно, траектория не обязана быть замкнутой, т.е. за период обращения угол

меняется и определяется, логично, как:
![$\Delta\varphi=2\int\limits_{r_m}^{R_m}\frac{M/r^2dr}{\sqrt{2\cdot m[E-U(r)]-\frac{M^2}{r^2}}}$ $\Delta\varphi=2\int\limits_{r_m}^{R_m}\frac{M/r^2dr}{\sqrt{2\cdot m[E-U(r)]-\frac{M^2}{r^2}}}$](https://dxdy-04.korotkov.co.uk/f/b/7/0/b70a0d602b603dde7a2064bfe30e5d6d82.png)
(2)
где

и

- максимальная и минимальная длина радиус-вектора соответственно.
И вот тут я не пойму. На
Википедии ничего про это не написано, а объясняется данное явления действием других тел, релятивистскими эффектами и приливными силами. В формуле же не учитывается ничего из перечисленного. Можно ли использовать формулу (2) для гравитационно связанных систем? Я считаю, что можно. Даже посчитал один частный случай: тело массы 100 кг находится на орбите с апоцентром 35000 км над поверхностью Земли и перицентром 180 км. У меня получилось, что за один оборот линия апсид сместится на 3 градуса.