RITE("JAHR");READLN;READ(J);
XYZ(D,MO,J);
WRITE("GROSSE HALBACHSE A");READLN;READ(A);
WRITE("BAHNEXZENTRIZITAET E");READLN;READ(F);
WRITE("DELTA T");READLN;READ(T);

Q:=A*(1-F);
TAU:=K*T;

(*++++Hilfsgroessen++++*)
GR:=F/Q/Q/Q*SQR (TAU);
CHI:=SQRT ((1-F)/Q/Q/Q)*TAU;


(*-------Startwert fuer E-------*)
E0:=0;C3:=1/6;

(*        Iterationsverfahren        *)
WHILE ABS (E-E0)>1E-7 DO BEGIN
E:=E0;
C3:=F3 (E);
H0:=1/C3/GR;H1:=H0/2;H2:=SQRT (1+4/27*H0);
Z:=EXP (LN (H1*H2+H1)/3)-EXP (LN (H1*H2-H1)/3);
E0:=CHI*Z
END;

R:=A*(1-F*COS (E));
V:=2*ARCTAN (SQRT ((1+F)/(1-F))*TAN (E/2));
WRITELN("RADIUSVEKTOR R:",R:10:6);
WRITELN("WAHRE ANOMALIE V:",V*180/PI:10:4," *");

(*=====Gauss'sche Konstanten=====*)
WRITE("LAENGE DES PERIHELS VOM AUFSTEIGENDEN KNOTEN");READLN;READ(W);
WRITE("LAENGE DES AUFSTEIGENDEN KNOTENS");READLN;READ(L);
I:=EKL;
W:=W*PI/180;L:=L*PI/180;I:=I*PI/180;
PX:=COS (W)*COS (L)-SIN (W)*SIN (L)*COS (I);
  H1:=COS (W)*SIN (L)+SIN (W)*COS (L)*COS (I);
  H2:=SIN (W)*SIN(I);
QX:=-SIN (W)*COS (L)-COS (W)*SIN (L)*COS (I);
  H3:=-SIN (W)*SIN (L)+COS (W)*COS (L)*COS (I);
  H4:=COS (W)*SIN (I);
  CE