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

Uses Math,KonstPLN;

Const ED5=9; EG5=108.6; ER5=0.9833; // початкові координати Землі 9 січня 1990р.

Procedure Eagth1(var EG,ER,ER0,ESD,EH,ER1,EGD,ED:real;var DR:integer);

implementation

Procedure Eagth1;

begin

ED:=ED5; EG:=EG5; ER:=ER5;

ER0:=Evo*86400/AO*N0; {формула 4} //шлях проходження за N3 дні

ESD:=(PI*sqr(Eao)*sqrt(1-sqr(Eeo))/Etd)*N0;//формула 3 площа еліпса орбіти за N3дні

Repeat

ED:=ED+N0; //2 закон Кеплера

EH:=2*ESD/ER; //радіус-вектор через N3 днів

ER1:=sqrt(sqr(EH)+sqr(ER-sqrt(abs(sqr(ER0)-sqr(EH)))));//формула 5

EGD:=arctan((1/ER1*EH)/sqrt(1-sqr(1/ER1*EH)))*180/PI; //формула 6

EG:=EG+EGD; //кут проходу за N3 дні

while EG>=360 do

begin

EG:=EG-360

end;

ER1:=(Eao*(1-sqr(Eeo)))/(1+Eeo*cos((EG-Eap)*PI/180)); //формула 2

ER:=ER1;

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

end;

end.

unit Mars;