a 6378137,000 a -3465340,282 ^2 0,00669437999014 Ya 2638267,760 Z_a 4644094,269 Вычислим эллипсоидальные координаты (B, L, H) точки наблюдений: =√(Xp^2+Yp^2+Zp^2;) =arcsin(Zp/R); =(e^2 α)/2R; Для вычисления широты (B) воспользуемся методом итераций и

a 6378137,000 \nX_a -3465340,282\ne^2 0,00669437999014 Y_a 2638267,760
Z_a 4644094,269

Вычислим эллипсоидальные координаты (B, L, H) точки наблюдений:\nR=√(X_p^2+Y_p^2+Z_p^2;) \nC=arcsin(Z_p/R); \nP=(e^2 α)/2R;
Для вычисления широты (B) воспользуемся методом итераций и следующими формулами:\nS_1=0; \nB_1=C+S_1; \nS_2=arcsin((P*Sin(2B_1 ))/√(1-e^2 Sin^2 B_1 )); \nE=|S_2-S_1 |;
Для вычисления используется следующий алгоритм:\nS_1=0;B_1=C+S_0;\nS_2=arcsin((P*Sin(2B_1 ))/√(1-e^2 〖Sin〗^2 B_1 ));〖 B〗_2=C+S_2; E=|S_2-S_1 |;\nS_3=arcsin((P*Sin(2B_2 ))/√(1-e^2 〖Sin〗^2 B_2 ));〖 B〗_3=C+S_3; E=|S_3-S_2 |;
…\nS_n=arcsin((P*Sin(2B_(n-1) ))/√(1-e^2 〖Sin〗^2 B_(n-1) ));〖 B〗_n=C+S_n; E=|S_n-S_(n-1) |;
Итерации следует продолжать, пока не выполнится условие E_n<1*10^(-7).
После того, как условие выполнено, используется следующая формула:\nB=B_n;