Перигелий Меркурия
Для определенности запустим по орбите вокруг центрального тела, например, протон, конструктивно оформленный в виде устойчивой структуры из двух синхронных и синфазных фотонов.
Орбитальная скорость

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

подставляем в (1) и получаем

фокальный радиус (справа)
, где
- расстояние от центра до перигелия
ускорение свободного падения на расстоянии
от тела массы 
- расстояние от фокуса с центральным телом до точки положения тела на орбите.
- время одной итерации
где
- Комптоновская длина волны протона, нейтрона, электрона
- скорость света.
- падение за одну большую итерацию в метрах
- условное падение за одну малую итерацию будет умножено на "тройку"
"Тройка" -
- зависит от расстояния между центрами фотонов, выраженное в единицах 

где
- средняя орбитальная скорость.
Сама функция плавно изменяется от
при
до
при
.


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

подставляем 

Осталось проинтегрировать
- угол смещения за половину оборота или
- угол смещения перицентра за оборот,
где 
Теперь нужно показать, как работать под арксинусом для любой 
Вещество на орбите горизонта событий черной дыры
Для проверки формулы движения перицентра и её природы запустим по орбите радиусом
тело. Пусть масса ЧД равна
массам Солнца. Орбитальная скорость, естественно, постоянна и равна
. Никаких особенностей орбита не имеет, эксцентриситет её равен нулю - это строго круговая орбита.
Теперь об ожиданиях от расчета. Поскольку при
обе итерации становятся равноправными, то за оборот они должны накопить каждая по
радиан. Но для приближения к реальности рассчитаем каждую из них раздельно.
Для большой итерации, которая работает всегда, при любой орбитальной скорости, и в случае движения перигелия Меркурия отвечает за основное орбитальное движения тела
, где
- гравитационная постоянная,
- масса черной дыры,
- Комптоновская длина волны протона,
- гравитационный, он же орбитальный радиус,
- Комптоновская частота протона,
- время полного орбитального оборота.
Для любой массы черной дыры
, что и ожидалось. Вторая часть должна быть тем, что мы называли движением перицентра 
- угол смещения перицентра за оборот,
где 
В случае ЧД формула значительно упрощается до
,
которая при подстановке значений дает 
Расчет периастра для двойного пульсара Халса-Тейлора (B1913+16)
Массы компонентов и орбитальные характеристики:
M⊙,
M⊙, эксцентриситет
,
время оборота
час =
сек, средняя скорость движения по орбите
км/с.
Из уравнения
можно выразить
- большую полуось
![Rendered by QuickLaTeX.com a=\sqrt[3]{\frac{GP^2(M_{1}+M_{2})}{4\pi ^2}}](https://photonphysics.ru/wp-content/ql-cache/quicklatex.com-854000e4af3e9b63885b7bce46400633_l3.png)
Подставляя значения в 
градусов в год.
При этом в уравнение
подставляется относительная скорость компонентов системы
км/сек.*
*Для таких малых скоростей
можно считать равной
.
Расчет периастра для двойного пульсара B2127+11C
Массы компонентов и орбитальные характеристики:
M⊙,
M⊙, эксцентриситет
,
время оборота
час
сек,
средняя скорость движения по орбите
км/с.
Пользуясь выкладками предыдущей задачи, получим
градусов в год.
Расчет периастра для двойного пульсара J1807−2500B
Массы компонентов и орбитальные характеристики:
M⊙,
M⊙, эксцентриситет
,
время оборота
суток
сек.
Пользуясь выкладками предыдущих задач, получим
угловой минуты в год.