Расчет подспутниковой точки спутника на ГСО
 

Помогите рассчитать подспутниковую точку КА на ГСО с прецессией. Т.е. спутник не "висит" в одной точки, а выписывает на орбите восьмерку. Задача состоит в том, чтобы определить по известным координатам станции наблюдения и разовому съему азимутальных и угломестных координат, подспутниковую точку КА. Существует формула для расчета ПСТ, но она корректно работает только для КА без прецессии:double LaES- географическая широта double LoES - географическая долгота double LoAzimut - азимут спутника double LoUgolMesta- угол места спутника d_sh2 - подспутниковая точкаdouble d_sh2=-atan(tan(LoAzimut)*sin(LaES))*180/M_PI+LoES*180/M_PI; Как видно из формулы в ней не используются угломестные координаты.Необходимая точность +-0.5 гр.