Тригонометричні ефемериди планет Сонячної системи

AR1:=(Aao*(1-sqr(Aeo)))/(1+Aeo*cos((AG-Aap)*PI/180));//формула 2

AR:=AR1;

x:=x+1 //кількість циклів обчислення

Until RM<=AD; //обчислення до моменту спостереження

end;//Mars1

Procedure Mars2;

begin

AAD1:=arctan((ER*sin((AG-EG)*PI/180))/(AR-ER*cos((AG-EG)*PI/180)))*180/PI;//зміщення прямого

ARS:=sin((AG-Avy)*PI/180)*sin(Aie*PI/180); //піднесення, формула 7

AZE:=AR*ARS;

AZP:=sqrt(sqr(ER*sin((EG-AG)*PI/180))+sqr(AR*sqrt(1-sqr(ARS))-ER*cos((EG-AG)*PI/180)));

ABD:=arctan(AZE/AZP)*180/PI; //зміщення схилення, формула 8

AV:=sqrt(sqr(AZE)+sqr(AZP));//лінійна відстань між Землею і Марсом, формула 9

AYD:=2*arctan(Are/(AV*AO))*180/PI;//кутовий діаметр Марса, формула 10

AYDS:=AYD*3600;

AAD0:=arctan(sqrt(1-sqr((sqr(AV)+sqr(AR)-sqr(ER))/(2*AV*AR)))/((sqr(AV)+sqr(AR)-sqr(ER))/

(2*AV*AR)))*180/PI;

AAD:=sqrt(abs(sqr(AAD0)-sqr(ABD)));

If AAD1<0 then

begin

AAD:=-AAD;

end;

AA:=(AG+AAD)/15; //пряме піднесення Марса

If AA>=24 then

begin