А это сообщение специально для Petrovich_Tot , который написал,
"открывай, у меня есть модель на Математике по релятивистским эффектам в двойных системах. Сравним."
Так вот, в своей статье "Влияние скорости гравитации на изменения параметров орбит планет"
http://modsys.narod.ru/Stat/Stat_Est/Vlijanie3.html http://modsys.narod.ru/Stat/Stat_Est/Vlijanie/Vlijanie3doc.zip я смоделировал движение двойного пульсара В1913+16 с использованием различных законов тяготения и в том числе с использованием ОТО Эйнштейна. Уравнения ОТО я взял у Ландавшицей (здесь помог Дробышев на форуме SciTecLibrary), но, как показали дальнейшие исследования, эти уравнения применимы только при малых скоростях и сейчас я уже пользуюсь уравнениями Сергея Хартикова, которые он самостоятельно вывел мне на Астрофоруме. Хотя для этого пульсара пойдут и уравнения Ландавшицей, поэтому ниже я привожу таблицу 7-3 из статьи (в сокращенном виде), где даны получающиеся смещения параметров орбиты этого пульсара (долготы перигелия AlfaP и восходящего узла AlfaU, угла наклона орбиты Betta, эксцентриситета Eks, большой полуоси Rsr) при применении этих уравнений и для сравнения привожу смещения, которые получаются по другим теориям, а так же привожу и свой комментарий из статьи к этой таблице. При этом параметры системы были следующие - масса пульсара 1,42*Мсолнца, масса компаньона 1,41*Мсолнца, большая полуось 971,437 тыс. км, эксцентриситет 0,617, долгота периастра 178,9 градуса, угол наклона орбиты к расчетной плоскости 42,9 градуса.
Таблица 7-3. Смещения параметров орбиты двойного пульсара PSR 1913+16 рассчитанные при скорости гравитации равной скорости света (Vgr=Vsv) для потенциалов Гербера, Лиенара-Вихерта (Л-В) и для ОТО и при Vgr=100*Vsv для потенциалов Юдина. Все эксперименты выполнены в покоящейся системе двойного пульсара (кроме нескольких при скорости системы VY=100 км/с) и при расчете запаздывающих координат для потенциалов Лиенара-Вихерта и Юдина по методике DX2, а также с основным шагом решения P0=3,6 с, который уменьшался в 100 раз вблизи периастра, узла восхождения и при переходе 0 градусов долготы.
Варианты расчета___________dAlfaP__________dAlfaU________dBetta__________dEks________________dR
ОТО Ландау___________1522189,5+/-16,1___4,6+/-2,0____-0,01+/-0,00______0,0+/-0,0___________0,0+/-0,0
ОТО Ландау___VY=100__1522190,5+/-19,7___3,2+/-2,1____-0,00+/-0,00______-0,1+/-0,1__________0,0+/-0,0
Гербер1_______________1522169,2+/-19,2___2,0+/-6,7____-0,00+/-0,00______403,3+/-0,1_________1,6+/-0,0
Гербер1______VY=100___1522172,6+/-22,3___2,6+/-3,6____-0,00+/-0,00______402,5+/-0,1_________1,6+/-0,0
Л-В полный____________-182326,8+/-79,6___4,7+/-30,4___-0,09+/-0,08____349474,3+/-164,8_____1356,7+/-0,6
Л-В полный___VY=100__-180993,2+/-68,8____8,5+/-26,1___-0,17+/-0,10____293913,9+/-105,6_____1138,3+/-0,5
Юдин k1=0 k2=0__________0,4+/-16,9_____-15,9+/-14,8____0,04+/-0,01___1787409,1+/-33697,7___31811,7+/-261,0
Юдин k1=0 k2=0 VY=100__21,2+/-6,9_______-77,0+/-7,4___-32,3+/-0,8____1691640,5+/-45463,3___31064,9+/-363,6
Юдин k1=1 k2=1_________-21,7+/-4,5_______-5,2+/-8,3____0,03+/-0,00___4747236,3+/-188749,9__48419,1+/-132,9
(см. продолжение)