Формулы для численного вычисления аномального смещения перигелиев без применения тензоров.
a - большая полуось орбиты в метрах
r_0 - гравитационный радиус Солнца в метрах (2953,25).
ε - эксцентриситет орбиты (Мерк. 0,20563593, Земля 0,01671123).
Интегралы не берущиеся.
Вычисление делается в электронной таблице за 500 шагов
через дельту = 0,0125663706143592 радиана.
Интегрирование по полному обороту.
L=a \int_{0}^{1} \sqrt{1-\varepsilon^2 \cos^2 (\varphi)} ~d (\varphi)
L1=\int_{0}^{1} a * \left (1- \frac{r_0}{a}\frac{(1- \varepsilon \cos (\varphi))}{(1-\varepsilon^2)}\right)^{-3/2} }*\sqrt{1-\varepsilon^2 \cos^2 ( \varphi)} ~d ( \varphi)
Смещение перигелия за сто лет в угловых секундах:
Δφsec= = \left(\frac{L_1 - L}{L} \right)
N — число оборотов планеты за сто лет.
Проверено для Меркурия (42,9) и Земли (3,8).
Проверка была и для двойных квазаров.
Подставляя параметры звезды и орбит планет, не сложно получить и угловую скорость их вращения.
Формула получена на основе гипотезы опубликованной на странице https://beard-studio.website.yandexcloud.net/
Гипотеза по нашему времени крамольная, но мне нравится. Да и формула на её основе рабочая.