tpsystools_4.04/source/StMerc.pas

// Upgraded to Delphi 2009: Sebastian Zierer


(* ***** BEGIN LICENSE BLOCK *****
 * Version: MPL 1.1
 *
 * The contents of this file are subject to the Mozilla Public License Version
 * 1.1 (the "License"); you may not use this file except in compliance with
 * the License. You may obtain a copy of the License at
 * http://www.mozilla.org/MPL/
 *
 * Software distributed under the License is distributed on an "AS IS" basis,
 * WITHOUT WARRANTY OF ANY KIND, either express or implied. See the License
 * for the specific language governing rights and limitations under the
 * License.
 *
 * The Original Code is TurboPower SysTools
 *
 * The Initial Developer of the Original Code is
 * TurboPower Software
 *
 * Portions created by the Initial Developer are Copyright (C) 1996-2002
 * the Initial Developer. All Rights Reserved.
 *
 * Contributor(s):
 *
 * ***** END LICENSE BLOCK ***** *)

{*********************************************************}
{* SysTools: StMerc.pas 4.04                             *}
{*********************************************************}
{* SysTools: Astronomical Routines (for Mercury)         *}
{*********************************************************}

{$I StDefine.inc}

unit StMerc;

interface

uses
  StAstroP;

function ComputeMercury(JD : Double) : TStEclipticalCord;


implementation

function GetLongitude(Tau, Tau2, Tau3, Tau4, Tau5 : Double) : Double;
var
  L0, L1,
  L2, L3,
  L4, L5   : Double;
begin
  L0 := 4.40250710140 * cos(0.00000000000 +      0.00000000000 * Tau)
      + 0.40989414976 * cos(1.48302034190 +  26087.90314200000 * Tau)
      + 0.05046294199 * cos(4.47785489540 +  52175.80628300000 * Tau)
      + 0.00855346843 * cos(1.16520322350 +  78263.70942500000 * Tau)
      + 0.00165590362 * cos(4.11969163180 + 104351.61257000000 * Tau)
      + 0.00034561897 * cos(0.77930765817 + 130439.51571000000 * Tau)
      + 0.00007583476 * cos(3.71348400510 + 156527.41885000000 * Tau)
      + 0.00003559740 * cos(1.51202669420 +   1109.37855210000 * Tau)
      + 0.00001803463 * cos(4.10333178410 +   5661.33204920000 * Tau)
      + 0.00001726012 * cos(0.35832239908 + 182615.32199000000 * Tau)
      + 0.00001589923 * cos(2.99510417810 +  25028.52121100000 * Tau)
      + 0.00001364682 * cos(4.59918318740 +  27197.28169400000 * Tau)
      + 0.00001017332 * cos(0.88031439040 +  31749.23519100000 * Tau)
      + 0.00000714182 * cos(1.54144865260 +  24978.52458900000 * Tau)
      + 0.00000643759 * cos(5.30266110790 +  21535.94964400000 * Tau)
      + 0.00000451137 * cos(6.04989275290 +  51116.42435300000 * Tau)
      + 0.00000404200 * cos(3.28228847030 + 208703.22513000000 * Tau)
      + 0.00000352441 * cos(5.24156297100 +  20426.57109200000 * Tau)
      + 0.00000345212 * cos(2.79211901540 +  15874.61759500000 * Tau)
      + 0.00000343313 * cos(5.76531885340 +    955.59974161000 * Tau)
      + 0.00000339214 * cos(5.86327765000 +  25558.21217600000 * Tau)
      + 0.00000325335 * cos(1.33674334780 +  53285.18483500000 * Tau)
      + 0.00000272947 * cos(2.49451163980 +    529.69096509000 * Tau)
      + 0.00000264336 * cos(3.91705094010 +  57837.13833200000 * Tau)
      + 0.00000259587 * cos(0.98732428184 +   4551.95349710000 * Tau)
      + 0.00000238793 * cos(0.11343953378 +   1059.38193020000 * Tau)
      + 0.00000234830 * cos(0.26672118900 +  11322.66409800000 * Tau)
      + 0.00000216645 * cos(0.65987207348 +  13521.75144200000 * Tau)
      + 0.00000208995 * cos(2.09178234010 +  47623.85278600000 * Tau)
      + 0.00000183359 * cos(2.62878670780 +  27043.50288300000 * Tau)
      + 0.00000181629 * cos(2.43413502470 +  25661.30495100000 * Tau)
      + 0.00000175965 * cos(4.53636829860 +  51066.42773100000 * Tau)
      + 0.00000172643 * cos(2.45200164170 +  24498.83024600000 * Tau)
      + 0.00000142316 * cos(3.36003948840 +  37410.56724000000 * Tau)
      + 0.00000137942 * cos(0.29098447849 +  10213.28554600000 * Tau)
      + 0.00000125219 * cos(3.72079804430 +  39609.65458300000 * Tau)
      + 0.00000118233 * cos(2.78149786370 +  77204.32749400000 * Tau)
      + 0.00000106422 * cos(4.20572116250 +  19804.82729200000 * Tau);

  L1 := 26088.1470620 * cos(0.00000000000 +      0.00000000000 * Tau)
      + 0.01126007832 * cos(6.21703971000 +  26087.90314200000 * Tau)
      + 0.00303471395 * cos(3.05565472360 +  52175.80628300000 * Tau)
      + 0.00080538452 * cos(6.10454743370 +  78263.70942500000 * Tau)
      + 0.00021245035 * cos(2.83531934450 + 104351.61257000000 * Tau)
      + 0.00005592094 * cos(5.82675673330 + 130439.51571000000 * Tau)
      + 0.00001472233 * cos(2.51845458400 + 156527.41885000000 * Tau)
      + 0.00000388318 * cos(5.48039225890 + 182615.32199000000 * Tau)
      + 0.00000352244 * cos(3.05238094400 +   1109.37855210000 * Tau)
      + 0.00000102743 * cos(2.14879173780 + 208703.22513000000 * Tau)
      + 0.00000093540 * cos(6.11791163930 +  27197.28169400000 * Tau)
      + 0.00000090579 * cos(0.00045481669 +  24978.52458900000 * Tau)
      + 0.00000051941 * cos(5.62107554050 +   5661.33204920000 * Tau)
      + 0.00000044370 * cos(4.57348500460 +  25028.52121100000 * Tau)
      + 0.00000028070 * cos(3.04195430990 +  51066.42773100000 * Tau)
      + 0.00000027295 * cos(5.09210138840 + 234791.12827000000 * Tau);

  L2 := 0.00053049845 * cos(0.00000000000 +      0.00000000000 * Tau)
      + 0.00016903658 * cos(4.69072300650 +  26087.90314200000 * Tau)
      + 0.00007396711 * cos(1.34735624670 +  52175.80628300000 * Tau)
      + 0.00003018297 * cos(4.45643539700 +  78263.70942500000 * Tau)
      + 0.00001107419 * cos(1.26226537550 + 104351.61257000000 * Tau)
      + 0.00000378173 * cos(4.31998055900 + 130439.51571000000 * Tau)
      + 0.00000122998 * cos(1.06868541050 + 156527.41885000000 * Tau)
      + 0.00000038663 * cos(4.08011610180 + 182615.32199000000 * Tau)
      + 0.00000014898 * cos(4.63343085810 +   1109.37855210000 * Tau)
      + 0.00000011861 * cos(0.79187646439 + 208703.22513000000 * Tau);

  L3 := 0.00000188077 * cos(0.03466830117 +   52175.80628300000 * Tau)
      + 0.00000142152 * cos(3.12505452600 +   26087.90314200000 * Tau)
      + 0.00000096877 * cos(3.00378171920 +   78263.70942500000 * Tau)
      + 0.00000043669 * cos(6.01867965830 +  104351.61257000000 * Tau)
      + 0.00000035395 * cos(0.00000000000 +       0.00000000000 * Tau)
      + 0.00000018045 * cos(2.77538373990 +  130439.51571000000 * Tau)
      + 0.00000006971 * cos(5.81808665740 +  156527.41885000000 * Tau)
      + 0.00000002556 * cos(2.57014364450 +  182615.32199000000 * Tau);

  L4 := 0.00000114078 * cos(3.14159265360 +       0.00000000000 * Tau)
      + 0.00000003247 * cos(2.02848007620 +   26087.90314200000 * Tau)
      + 0.00000001914 * cos(1.41731803760 +   78263.70942500000 * Tau)
      + 0.00000001727 * cos(4.50137643800 +   52175.80628300000 * Tau)
      + 0.00000001237 * cos(4.49970181060 +  104351.61257000000 * Tau)
      + 0.00000000645 * cos(1.26591776990 +  130439.51571000000 * Tau);

  L5 := 0.00000000877 * cos(3.14159265360 +  0.00000000000 * Tau);
  Result := (L0 + L1*Tau + L2*Tau2 + L3*Tau3 + L4*Tau4 + L5*Tau5);
end;

{-------------------------------------------------------------------------}

function GetLatitude(Tau, Tau2, Tau3, Tau4, Tau5 : Double) : Double;
var
  B0, B1,
  B2, B3,
  B4, B5   : Double;
begin

  B0 := 0.11737528962 * cos(1.98357498770 +  26087.90314200000 * Tau)
      + 0.02388076996 * cos(5.03738959690 +  52175.80628300000 * Tau)
      + 0.01222839532 * cos(3.14159265360 +      0.00000000000 * Tau)
      + 0.00543251810 * cos(1.79644363960 +  78263.70942500000 * Tau)
      + 0.00129778770 * cos(4.83232503960 + 104351.61257000000 * Tau)
      + 0.00031866927 * cos(1.58088495670 + 130439.51571000000 * Tau)
      + 0.00007963301 * cos(4.60972126350 + 156527.41885000000 * Tau)
      + 0.00002014189 * cos(1.35324164690 + 182615.32199000000 * Tau)
      + 0.00000513953 * cos(4.37835409310 + 208703.22513000000 * Tau)
      + 0.00000208584 * cos(2.02020294150 +  24978.52458900000 * Tau)
      + 0.00000207674 * cos(4.91772564070 +  27197.28169400000 * Tau)
      + 0.00000132013 * cos(1.11908492280 + 234791.12827000000 * Tau)
      + 0.00000121395 * cos(1.81271752060 +  53285.18483500000 * Tau)
      + 0.00000100454 * cos(5.65684734210 +  20426.57109200000 * Tau);

  B1 := 0.00429151362 * cos(3.50169780390 +  26087.90314200000 * Tau)
      + 0.00146233668 * cos(3.14159265360 +      0.00000000000 * Tau)
      + 0.00022675295 * cos(0.01515366880 +  52175.80628300000 * Tau)
      + 0.00010894981 * cos(0.48540174006 +  78263.70942500000 * Tau)
      + 0.00006353462 * cos(3.42943919980 + 104351.61257000000 * Tau)
      + 0.00002495743 * cos(0.16051210665 + 130439.51571000000 * Tau)
      + 0.00000859585 * cos(3.18452433650 + 156527.41885000000 * Tau)
      + 0.00000277503 * cos(6.21020774180 + 182615.32199000000 * Tau)
      + 0.00000086233 * cos(2.95244391820 + 208703.22513000000 * Tau)
      + 0.00000027696 * cos(0.29068938889 +  27197.28169400000 * Tau)
      + 0.00000026133 * cos(5.97708962690 + 234791.12827000000 * Tau);

  B2 := 0.00011830934 * cos(4.79065585780 +  26087.90314200000 * Tau)
      + 0.00001913516 * cos(0.00000000000 +      0.00000000000 * Tau)
      + 0.00001044801 * cos(1.21216540540 +  52175.80628300000 * Tau)
      + 0.00000266213 * cos(4.43418336530 +  78263.70942500000 * Tau)
      + 0.00000170280 * cos(1.62255638710 + 104351.61257000000 * Tau)
      + 0.00000096300 * cos(4.80023692020 + 130439.51571000000 * Tau)
      + 0.00000044692 * cos(1.60758267770 + 156527.41885000000 * Tau)
      + 0.00000018316 * cos(4.66904655380 + 182615.32199000000 * Tau)
      + 0.00000006927 * cos(1.43404888930 + 208703.22513000000 * Tau);

  B3 := 0.00000235423 * cos(0.35387524604 +  26087.90314200000 * Tau)
      + 0.00000160537 * cos(0.00000000000 +      0.00000000000 * Tau)
      + 0.00000018904 * cos(4.36275460260 +  52175.80628300000 * Tau)
      + 0.00000006376 * cos(2.50715381440 +  78263.70942500000 * Tau)
      + 0.00000004580 * cos(6.14257817570 + 104351.61257000000 * Tau)
      + 0.00000003061 * cos(3.12497552680 + 130439.51571000000 * Tau)
      + 0.00000001732 * cos(6.26642412060 + 156527.41885000000 * Tau);

  B4 := 0.00000004276 * cos(1.74579932120 +  26087.90314200000 * Tau)
      + 0.00000001023 * cos(3.14159265360 +      0.00000000000 * Tau);

  B5 := 0.00000000000 * cos(0.00000000000 +      0.00000000000 * Tau);
  Result := (B0 + B1*Tau + B2*Tau2 + B3*Tau3 + B4*Tau4 + B5*Tau5);
end;

{-------------------------------------------------------------------------}

function GetRadiusVector(Tau, Tau2, Tau3, Tau4, Tau5 : Double) : Double;
var
  R0, R1,
  R2, R3,
  R4, R5  : Double;
begin
  R0 := 0.39528271652 * cos(0.00000000000 +      0.00000000000 * Tau)
      + 0.07834131817 * cos(6.19233722600 +  26087.90314200000 * Tau)
      + 0.00795525557 * cos(2.95989690100 +  52175.80628300000 * Tau)
      + 0.00121281763 * cos(6.01064153810 +  78263.70942500000 * Tau)
      + 0.00021921969 * cos(2.77820093970 + 104351.61257000000 * Tau)
      + 0.00004354065 * cos(5.82894543260 + 130439.51571000000 * Tau)
      + 0.00000918228 * cos(2.59650562600 + 156527.41885000000 * Tau)
      + 0.00000289955 * cos(1.42441936950 +  25028.52121100000 * Tau)
      + 0.00000260033 * cos(3.02817753480 +  27197.28169400000 * Tau)
      + 0.00000201855 * cos(5.64725040350 + 182615.32199000000 * Tau)
      + 0.00000201499 * cos(5.59227724200 +  31749.23519100000 * Tau)
      + 0.00000141980 * cos(6.25264202640 +  24978.52458900000 * Tau)
      + 0.00000100144 * cos(3.73435608690 +  21535.94964400000 * Tau);

  R1 := 0.00217347739 * cos(4.65617158660 +  26087.90314200000 * Tau)
      + 0.00044141826 * cos(1.42385543980 +  52175.80628300000 * Tau)
      + 0.00010094479 * cos(4.47466326320 +  78263.70942500000 * Tau)
      + 0.00002432804 * cos(1.24226083430 + 104351.61257000000 * Tau)
      + 0.00001624367 * cos(0.00000000000 +      0.00000000000 * Tau)
      + 0.00000603996 * cos(4.29303116560 + 130439.51571000000 * Tau)
      + 0.00000152851 * cos(1.06060779810 + 156527.41885000000 * Tau)
      + 0.00000039202 * cos(4.11136751420 + 182615.32199000000 * Tau);

  R2 := 0.00003117867 * cos(3.08231840300 +  26087.90314200000 * Tau)
      + 0.00001245396 * cos(6.15183317420 +  52175.80628300000 * Tau)
      + 0.00000424822 * cos(2.92583352960 +  78263.70942500000 * Tau)
      + 0.00000136130 * cos(5.97983925840 + 104351.61257000000 * Tau)
      + 0.00000042175 * cos(2.74936980630 + 130439.51571000000 * Tau)
      + 0.00000021759 * cos(3.14159265360 +      0.00000000000 * Tau)
      + 0.00000012793 * cos(5.80143162210 + 156527.41885000000 * Tau);

  R3 := 0.00000032676 * cos(1.67971635360 +  26087.90314200000 * Tau)
      + 0.00000024166 * cos(4.63403169000 +  52175.80628300000 * Tau)
      + 0.00000012133 * cos(1.38983781540 +  78263.70942500000 * Tau)
      + 0.00000005140 * cos(4.43915386930 + 104351.61257000000 * Tau)
      + 0.00000001981 * cos(1.20733880270 + 130439.51571000000 * Tau);

  R4 := 0.00000000000;

  R5 := 0.00000000000;
  Result := (R0 + R1*Tau + R2*Tau2 + R3*Tau3 + R4*Tau4 + R5*Tau5);
end;

{---------------------------------------------------------------------------}

function ComputeMercury(JD : Double) : TStEclipticalCord;
var
  Tau,
  Tau2,
  Tau3,
  Tau4,
  Tau5      : Double;
begin
  Tau  := (JD - 2451545.0) / 365250.0;
  Tau2 := sqr(Tau);
  Tau3 := Tau * Tau2;
  Tau4 := sqr(Tau2);
  Tau5 := Tau2 * Tau3;

  Result.L0 := GetLongitude(Tau, Tau2, Tau3, Tau4, Tau5);
  Result.B0 := GetLatitude(Tau, Tau2, Tau3, Tau4, Tau5);
  Result.R0 := GetRadiusVector(Tau, Tau2, Tau3, Tau4, Tau5);
end;


end.
1	torben	2671	// Upgraded to Delphi 2009: Sebastian Zierer
2
3
4			(* *** BEGIN LICENSE BLOCK ***
5			* Version: MPL 1.1
6			*
7			* The contents of this file are subject to the Mozilla Public License Version
8			* 1.1 (the "License"); you may not use this file except in compliance with
9			* the License. You may obtain a copy of the License at
10			* http://www.mozilla.org/MPL/
11			*
12			* Software distributed under the License is distributed on an "AS IS" basis,
13			* WITHOUT WARRANTY OF ANY KIND, either express or implied. See the License
14			* for the specific language governing rights and limitations under the
15			* License.
16			*
17			* The Original Code is TurboPower SysTools
18			*
19			* The Initial Developer of the Original Code is
20			* TurboPower Software
21			*
22			* Portions created by the Initial Developer are Copyright (C) 1996-2002
23			* the Initial Developer. All Rights Reserved.
24			*
25			* Contributor(s):
26			*
27			* *** END LICENSE BLOCK *** *)
28
29			{*********************************************************}
30			{* SysTools: StMerc.pas 4.04 *}
31			{*********************************************************}
32			{* SysTools: Astronomical Routines (for Mercury) *}
33			{*********************************************************}
34
35			{$I StDefine.inc}
36
37			unit StMerc;
38
39			interface
40
41			uses
42			StAstroP;
43
44			function ComputeMercury(JD : Double) : TStEclipticalCord;
45
46
47			implementation
48
49			function GetLongitude(Tau, Tau2, Tau3, Tau4, Tau5 : Double) : Double;
50			var
51			L0, L1,
52			L2, L3,
53			L4, L5 : Double;
54			begin
55			L0 := 4.40250710140 * cos(0.00000000000 + 0.00000000000 * Tau)
56			+ 0.40989414976 * cos(1.48302034190 + 26087.90314200000 * Tau)
57			+ 0.05046294199 * cos(4.47785489540 + 52175.80628300000 * Tau)
58			+ 0.00855346843 * cos(1.16520322350 + 78263.70942500000 * Tau)
59			+ 0.00165590362 * cos(4.11969163180 + 104351.61257000000 * Tau)
60			+ 0.00034561897 * cos(0.77930765817 + 130439.51571000000 * Tau)
61			+ 0.00007583476 * cos(3.71348400510 + 156527.41885000000 * Tau)
62			+ 0.00003559740 * cos(1.51202669420 + 1109.37855210000 * Tau)
63			+ 0.00001803463 * cos(4.10333178410 + 5661.33204920000 * Tau)
64			+ 0.00001726012 * cos(0.35832239908 + 182615.32199000000 * Tau)
65			+ 0.00001589923 * cos(2.99510417810 + 25028.52121100000 * Tau)
66			+ 0.00001364682 * cos(4.59918318740 + 27197.28169400000 * Tau)
67			+ 0.00001017332 * cos(0.88031439040 + 31749.23519100000 * Tau)
68			+ 0.00000714182 * cos(1.54144865260 + 24978.52458900000 * Tau)
69			+ 0.00000643759 * cos(5.30266110790 + 21535.94964400000 * Tau)
70			+ 0.00000451137 * cos(6.04989275290 + 51116.42435300000 * Tau)
71			+ 0.00000404200 * cos(3.28228847030 + 208703.22513000000 * Tau)
72			+ 0.00000352441 * cos(5.24156297100 + 20426.57109200000 * Tau)
73			+ 0.00000345212 * cos(2.79211901540 + 15874.61759500000 * Tau)
74			+ 0.00000343313 * cos(5.76531885340 + 955.59974161000 * Tau)
75			+ 0.00000339214 * cos(5.86327765000 + 25558.21217600000 * Tau)
76			+ 0.00000325335 * cos(1.33674334780 + 53285.18483500000 * Tau)
77			+ 0.00000272947 * cos(2.49451163980 + 529.69096509000 * Tau)
78			+ 0.00000264336 * cos(3.91705094010 + 57837.13833200000 * Tau)
79			+ 0.00000259587 * cos(0.98732428184 + 4551.95349710000 * Tau)
80			+ 0.00000238793 * cos(0.11343953378 + 1059.38193020000 * Tau)
81			+ 0.00000234830 * cos(0.26672118900 + 11322.66409800000 * Tau)
82			+ 0.00000216645 * cos(0.65987207348 + 13521.75144200000 * Tau)
83			+ 0.00000208995 * cos(2.09178234010 + 47623.85278600000 * Tau)
84			+ 0.00000183359 * cos(2.62878670780 + 27043.50288300000 * Tau)
85			+ 0.00000181629 * cos(2.43413502470 + 25661.30495100000 * Tau)
86			+ 0.00000175965 * cos(4.53636829860 + 51066.42773100000 * Tau)
87			+ 0.00000172643 * cos(2.45200164170 + 24498.83024600000 * Tau)
88			+ 0.00000142316 * cos(3.36003948840 + 37410.56724000000 * Tau)
89			+ 0.00000137942 * cos(0.29098447849 + 10213.28554600000 * Tau)
90			+ 0.00000125219 * cos(3.72079804430 + 39609.65458300000 * Tau)
91			+ 0.00000118233 * cos(2.78149786370 + 77204.32749400000 * Tau)
92			+ 0.00000106422 * cos(4.20572116250 + 19804.82729200000 * Tau);
93
94			L1 := 26088.1470620 * cos(0.00000000000 + 0.00000000000 * Tau)
95			+ 0.01126007832 * cos(6.21703971000 + 26087.90314200000 * Tau)
96			+ 0.00303471395 * cos(3.05565472360 + 52175.80628300000 * Tau)
97			+ 0.00080538452 * cos(6.10454743370 + 78263.70942500000 * Tau)
98			+ 0.00021245035 * cos(2.83531934450 + 104351.61257000000 * Tau)
99			+ 0.00005592094 * cos(5.82675673330 + 130439.51571000000 * Tau)
100			+ 0.00001472233 * cos(2.51845458400 + 156527.41885000000 * Tau)
101			+ 0.00000388318 * cos(5.48039225890 + 182615.32199000000 * Tau)
102			+ 0.00000352244 * cos(3.05238094400 + 1109.37855210000 * Tau)
103			+ 0.00000102743 * cos(2.14879173780 + 208703.22513000000 * Tau)
104			+ 0.00000093540 * cos(6.11791163930 + 27197.28169400000 * Tau)
105			+ 0.00000090579 * cos(0.00045481669 + 24978.52458900000 * Tau)
106			+ 0.00000051941 * cos(5.62107554050 + 5661.33204920000 * Tau)
107			+ 0.00000044370 * cos(4.57348500460 + 25028.52121100000 * Tau)
108			+ 0.00000028070 * cos(3.04195430990 + 51066.42773100000 * Tau)
109			+ 0.00000027295 * cos(5.09210138840 + 234791.12827000000 * Tau);
110
111			L2 := 0.00053049845 * cos(0.00000000000 + 0.00000000000 * Tau)
112			+ 0.00016903658 * cos(4.69072300650 + 26087.90314200000 * Tau)
113			+ 0.00007396711 * cos(1.34735624670 + 52175.80628300000 * Tau)
114			+ 0.00003018297 * cos(4.45643539700 + 78263.70942500000 * Tau)
115			+ 0.00001107419 * cos(1.26226537550 + 104351.61257000000 * Tau)
116			+ 0.00000378173 * cos(4.31998055900 + 130439.51571000000 * Tau)
117			+ 0.00000122998 * cos(1.06868541050 + 156527.41885000000 * Tau)
118			+ 0.00000038663 * cos(4.08011610180 + 182615.32199000000 * Tau)
119			+ 0.00000014898 * cos(4.63343085810 + 1109.37855210000 * Tau)
120			+ 0.00000011861 * cos(0.79187646439 + 208703.22513000000 * Tau);
121
122			L3 := 0.00000188077 * cos(0.03466830117 + 52175.80628300000 * Tau)
123			+ 0.00000142152 * cos(3.12505452600 + 26087.90314200000 * Tau)
124			+ 0.00000096877 * cos(3.00378171920 + 78263.70942500000 * Tau)
125			+ 0.00000043669 * cos(6.01867965830 + 104351.61257000000 * Tau)
126			+ 0.00000035395 * cos(0.00000000000 + 0.00000000000 * Tau)
127			+ 0.00000018045 * cos(2.77538373990 + 130439.51571000000 * Tau)
128			+ 0.00000006971 * cos(5.81808665740 + 156527.41885000000 * Tau)
129			+ 0.00000002556 * cos(2.57014364450 + 182615.32199000000 * Tau);
130
131			L4 := 0.00000114078 * cos(3.14159265360 + 0.00000000000 * Tau)
132			+ 0.00000003247 * cos(2.02848007620 + 26087.90314200000 * Tau)
133			+ 0.00000001914 * cos(1.41731803760 + 78263.70942500000 * Tau)
134			+ 0.00000001727 * cos(4.50137643800 + 52175.80628300000 * Tau)
135			+ 0.00000001237 * cos(4.49970181060 + 104351.61257000000 * Tau)
136			+ 0.00000000645 * cos(1.26591776990 + 130439.51571000000 * Tau);
137
138			L5 := 0.00000000877 * cos(3.14159265360 + 0.00000000000 * Tau);
139			Result := (L0 + L1Tau + L2Tau2 + L3Tau3 + L4Tau4 + L5*Tau5);
140			end;
141
142			{-------------------------------------------------------------------------}
143
144			function GetLatitude(Tau, Tau2, Tau3, Tau4, Tau5 : Double) : Double;
145			var
146			B0, B1,
147			B2, B3,
148			B4, B5 : Double;
149			begin
150
151			B0 := 0.11737528962 * cos(1.98357498770 + 26087.90314200000 * Tau)
152			+ 0.02388076996 * cos(5.03738959690 + 52175.80628300000 * Tau)
153			+ 0.01222839532 * cos(3.14159265360 + 0.00000000000 * Tau)
154			+ 0.00543251810 * cos(1.79644363960 + 78263.70942500000 * Tau)
155			+ 0.00129778770 * cos(4.83232503960 + 104351.61257000000 * Tau)
156			+ 0.00031866927 * cos(1.58088495670 + 130439.51571000000 * Tau)
157			+ 0.00007963301 * cos(4.60972126350 + 156527.41885000000 * Tau)
158			+ 0.00002014189 * cos(1.35324164690 + 182615.32199000000 * Tau)
159			+ 0.00000513953 * cos(4.37835409310 + 208703.22513000000 * Tau)
160			+ 0.00000208584 * cos(2.02020294150 + 24978.52458900000 * Tau)
161			+ 0.00000207674 * cos(4.91772564070 + 27197.28169400000 * Tau)
162			+ 0.00000132013 * cos(1.11908492280 + 234791.12827000000 * Tau)
163			+ 0.00000121395 * cos(1.81271752060 + 53285.18483500000 * Tau)
164			+ 0.00000100454 * cos(5.65684734210 + 20426.57109200000 * Tau);
165
166			B1 := 0.00429151362 * cos(3.50169780390 + 26087.90314200000 * Tau)
167			+ 0.00146233668 * cos(3.14159265360 + 0.00000000000 * Tau)
168			+ 0.00022675295 * cos(0.01515366880 + 52175.80628300000 * Tau)
169			+ 0.00010894981 * cos(0.48540174006 + 78263.70942500000 * Tau)
170			+ 0.00006353462 * cos(3.42943919980 + 104351.61257000000 * Tau)
171			+ 0.00002495743 * cos(0.16051210665 + 130439.51571000000 * Tau)
172			+ 0.00000859585 * cos(3.18452433650 + 156527.41885000000 * Tau)
173			+ 0.00000277503 * cos(6.21020774180 + 182615.32199000000 * Tau)
174			+ 0.00000086233 * cos(2.95244391820 + 208703.22513000000 * Tau)
175			+ 0.00000027696 * cos(0.29068938889 + 27197.28169400000 * Tau)
176			+ 0.00000026133 * cos(5.97708962690 + 234791.12827000000 * Tau);
177
178			B2 := 0.00011830934 * cos(4.79065585780 + 26087.90314200000 * Tau)
179			+ 0.00001913516 * cos(0.00000000000 + 0.00000000000 * Tau)
180			+ 0.00001044801 * cos(1.21216540540 + 52175.80628300000 * Tau)
181			+ 0.00000266213 * cos(4.43418336530 + 78263.70942500000 * Tau)
182			+ 0.00000170280 * cos(1.62255638710 + 104351.61257000000 * Tau)
183			+ 0.00000096300 * cos(4.80023692020 + 130439.51571000000 * Tau)
184			+ 0.00000044692 * cos(1.60758267770 + 156527.41885000000 * Tau)
185			+ 0.00000018316 * cos(4.66904655380 + 182615.32199000000 * Tau)
186			+ 0.00000006927 * cos(1.43404888930 + 208703.22513000000 * Tau);
187
188			B3 := 0.00000235423 * cos(0.35387524604 + 26087.90314200000 * Tau)
189			+ 0.00000160537 * cos(0.00000000000 + 0.00000000000 * Tau)
190			+ 0.00000018904 * cos(4.36275460260 + 52175.80628300000 * Tau)
191			+ 0.00000006376 * cos(2.50715381440 + 78263.70942500000 * Tau)
192			+ 0.00000004580 * cos(6.14257817570 + 104351.61257000000 * Tau)
193			+ 0.00000003061 * cos(3.12497552680 + 130439.51571000000 * Tau)
194			+ 0.00000001732 * cos(6.26642412060 + 156527.41885000000 * Tau);
195
196			B4 := 0.00000004276 * cos(1.74579932120 + 26087.90314200000 * Tau)
197			+ 0.00000001023 * cos(3.14159265360 + 0.00000000000 * Tau);
198
199			B5 := 0.00000000000 * cos(0.00000000000 + 0.00000000000 * Tau);
200			Result := (B0 + B1Tau + B2Tau2 + B3Tau3 + B4Tau4 + B5*Tau5);
201			end;
202
203			{-------------------------------------------------------------------------}
204
205			function GetRadiusVector(Tau, Tau2, Tau3, Tau4, Tau5 : Double) : Double;
206			var
207			R0, R1,
208			R2, R3,
209			R4, R5 : Double;
210			begin
211			R0 := 0.39528271652 * cos(0.00000000000 + 0.00000000000 * Tau)
212			+ 0.07834131817 * cos(6.19233722600 + 26087.90314200000 * Tau)
213			+ 0.00795525557 * cos(2.95989690100 + 52175.80628300000 * Tau)
214			+ 0.00121281763 * cos(6.01064153810 + 78263.70942500000 * Tau)
215			+ 0.00021921969 * cos(2.77820093970 + 104351.61257000000 * Tau)
216			+ 0.00004354065 * cos(5.82894543260 + 130439.51571000000 * Tau)
217			+ 0.00000918228 * cos(2.59650562600 + 156527.41885000000 * Tau)
218			+ 0.00000289955 * cos(1.42441936950 + 25028.52121100000 * Tau)
219			+ 0.00000260033 * cos(3.02817753480 + 27197.28169400000 * Tau)
220			+ 0.00000201855 * cos(5.64725040350 + 182615.32199000000 * Tau)
221			+ 0.00000201499 * cos(5.59227724200 + 31749.23519100000 * Tau)
222			+ 0.00000141980 * cos(6.25264202640 + 24978.52458900000 * Tau)
223			+ 0.00000100144 * cos(3.73435608690 + 21535.94964400000 * Tau);
224
225			R1 := 0.00217347739 * cos(4.65617158660 + 26087.90314200000 * Tau)
226			+ 0.00044141826 * cos(1.42385543980 + 52175.80628300000 * Tau)
227			+ 0.00010094479 * cos(4.47466326320 + 78263.70942500000 * Tau)
228			+ 0.00002432804 * cos(1.24226083430 + 104351.61257000000 * Tau)
229			+ 0.00001624367 * cos(0.00000000000 + 0.00000000000 * Tau)
230			+ 0.00000603996 * cos(4.29303116560 + 130439.51571000000 * Tau)
231			+ 0.00000152851 * cos(1.06060779810 + 156527.41885000000 * Tau)
232			+ 0.00000039202 * cos(4.11136751420 + 182615.32199000000 * Tau);
233
234			R2 := 0.00003117867 * cos(3.08231840300 + 26087.90314200000 * Tau)
235			+ 0.00001245396 * cos(6.15183317420 + 52175.80628300000 * Tau)
236			+ 0.00000424822 * cos(2.92583352960 + 78263.70942500000 * Tau)
237			+ 0.00000136130 * cos(5.97983925840 + 104351.61257000000 * Tau)
238			+ 0.00000042175 * cos(2.74936980630 + 130439.51571000000 * Tau)
239			+ 0.00000021759 * cos(3.14159265360 + 0.00000000000 * Tau)
240			+ 0.00000012793 * cos(5.80143162210 + 156527.41885000000 * Tau);
241
242			R3 := 0.00000032676 * cos(1.67971635360 + 26087.90314200000 * Tau)
243			+ 0.00000024166 * cos(4.63403169000 + 52175.80628300000 * Tau)
244			+ 0.00000012133 * cos(1.38983781540 + 78263.70942500000 * Tau)
245			+ 0.00000005140 * cos(4.43915386930 + 104351.61257000000 * Tau)
246			+ 0.00000001981 * cos(1.20733880270 + 130439.51571000000 * Tau);
247
248			R4 := 0.00000000000;
249
250			R5 := 0.00000000000;
251			Result := (R0 + R1Tau + R2Tau2 + R3Tau3 + R4Tau4 + R5*Tau5);
252			end;
253
254			{---------------------------------------------------------------------------}
255
256			function ComputeMercury(JD : Double) : TStEclipticalCord;
257			var
258			Tau,
259			Tau2,
260			Tau3,
261			Tau4,
262			Tau5 : Double;
263			begin
264			Tau := (JD - 2451545.0) / 365250.0;
265			Tau2 := sqr(Tau);
266			Tau3 := Tau * Tau2;
267			Tau4 := sqr(Tau2);
268			Tau5 := Tau2 * Tau3;
269
270			Result.L0 := GetLongitude(Tau, Tau2, Tau3, Tau4, Tau5);
271			Result.B0 := GetLatitude(Tau, Tau2, Tau3, Tau4, Tau5);
272			Result.R0 := GetRadiusVector(Tau, Tau2, Tau3, Tau4, Tau5);
273			end;
274
275
276			end.