Hola
Esto es lo que hace LBR_For_Mercury...hay un procedure por cada planeta
Código Delphi
[-]
procedure LBR_For_Mercury(At_JDE:extended; out L: extended; out B: extended; out R: extended);
var
T: extended;
begin
T := (At_JDE - 2451545) / 365250;
L := Mercury_L01(T) + Mercury_L02(T);
B := Mercury_B01(T) + Mercury_B02(T);
R := Mercury_R01(T) + Mercury_R02(T);
L := L * 45 / Arctan(1);
If Abs(L) >= 360 Then L := L - 360 * Int(L / 360);
If L < 0 Then L := L + 360;
B := B * 45 / Arctan(1);
end;