Вычисление постоянных интегрирования по заданным координатам и компонентам скорости спутника

Для вычислений была использована программа Excel. Результаты всех вычислений представлены в таблице ниже.
№ п/п
формулы
вычисления
Для t1=1530 сек
Для t2=1620 сек
1 Постоянные интегралов площадей c1=ynzn-znyn
26156995582,952 26156995582,952
c2=znxn-xnzn
-45305240045,133 -45305240045,133
c3=xnyn-ynxn
913149411,676 913149411,676
2 Постоянная интеграла энергии rn=xn2+yn2+zn2
7616,871 7616,871
rn=xn2+yn2+zn2
6869222,629 6869222,629
h=rn2-2μrn
-58048827,981 -58048827,981
3 Постоянные интегралов Лапласа f1=-μxnrn+yn∙c3-zn∙c2
201685283765,386 201685283765,489
f2=-μynrn+zn∙c1-xn∙c3
123130483893,914 123130483893,991
f3=-μznrn+xn∙c2-yn∙c1
331791354000,542 331791354000,367
Контроль:
- постоянные интегралов площадей:
c1, 1=c1, 2=c1=26156995582,952
c2, 1=c2, 2=c2=-45305240045,133
c3, 1=c3, 2=c3=913149411,676
- постоянная интеграла энергии:
h1=h2=-58048827,981
- постоянные интегралов Лапласа:
f1, 1≅f1,2≅f1, среднее=201685283765,44
f2, 1≅f2,2≅f2, среднее=123130483893,95
f3, 1≅f3,2≅f3, среднее=331791354000,45