[compuin]

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
  //Compute T as Julian millennia since J2000.0
  T := (At_JDE - 2451545) / 365250;
       
  //Compute heliocentric, ecliptical Longitude L in radians
  L := Mercury_L01(T) + Mercury_L02(T);

  //Compute heliocentric, ecliptical Latitude B in radians
  B := Mercury_B01(T) + Mercury_B02(T);

  //Compute heliocentric distance R in AU
  R := Mercury_R01(T) + Mercury_R02(T);   
  
  //Convert longitude and latitude from radians into degrees.
  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;