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

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;