Вычисление элементов орбиты ИСЗ по координатам и составляющим скорости

Вычисляем постоянные площадей
с1=y·z-y∙z
с1=1389427.298 ·-18123.244--5308.053∙(893020.655)=-20440728562.186 м2/с
с2=z·x-z∙x
с2=(893020.655)·-3233.033--18123.244∙-787090.954=-17151807787.330 м2/с
с3=x·y-x∙y
с3=-787090.954·-5308.053--3233.033∙1389427.298 =8669986929.620м2/с
c=с12+с22+с32
c=-20440728562.1862+-17151807787.3302+8669986929.6202=28056667084.426
2. Вычисляем фокальный параметр
p=c2μ
p=28056667084.42623.9860044∙1014=1974851.227 м
3. Вычисляем угол наклона плоскости орбиты к плоскости экватора
i=arccosс3c
i=arccos8669986929.62028056667084.426=72°00΄00˝
4. Вычисляем долготу восходящего узла Ω
sin Ω=c1c∙sini
sin Ω=-20440728562.18628056667084.426 ∙sin72°00΄00˝=-0.766044443
cos Ω=--17151807787.330 28056667084.426 ∙sin72°00΄00˝=0.642787609
Ω=130º00΄00˝
5 . Вычисляем квадрат скорости спутника и геоцентрическое расстояние до спутника
v2=x2+y2+z2
v2=(-3233.033 )2+-5308.0532+-18123.2442= 367079929.665 м2/с2
r=x2+y2+z2
r=(-787090.954)2+1389427.298 2+(893020.655)2=1829619.161 м
6.Вычисляем постоянную энергии
h=v2-2∙μr
h=367079929.665-2∙3.9860044∙10141829619.161=-68639643.453 м2/с2
7. Вычисляем большую полуось орбиты спутника
a=-μh
a=-3.9860044∙1014-68639643.453=5807146.132 м
8. Вычисляем составляющие вектора Лапласа
D=x∙x+y∙y+z∙z
D=-21014894728.687
D=μr+h
D=3.9860044∙10141829619.161+-68639643.453=149220143.106
f1=D∙x-D∙x
f1= -185391693907996.000
f2=D∙y-D∙y
f2=95782345517287.400
f3=D∙z-D∙z
f3=-247601401311347.000
f=f12+f22+f32
f=-1853916939079962+(95782345517287.4)2+(-247601401311347)2=323806719841037.000
9