/[projects]/dao/DelphiScanner/Components/tpsystools_4.04/source/StUranus.pas
ViewVC logotype

Contents of /dao/DelphiScanner/Components/tpsystools_4.04/source/StUranus.pas

Parent Directory Parent Directory | Revision Log Revision Log

Revision 2671 - (show annotations) (download)
Tue Aug 25 18:15:15 2015 UTC (8 years, 9 months ago) by torben
File size: 29586 byte(s) 
Added tpsystools component
1 // Upgraded to Delphi 2009: Sebastian Zierer
2
3 (* ***** BEGIN LICENSE BLOCK *****
4 * Version: MPL 1.1
5 *
6 * The contents of this file are subject to the Mozilla Public License Version
7 * 1.1 (the "License"); you may not use this file except in compliance with
8 * the License. You may obtain a copy of the License at
9 * http://www.mozilla.org/MPL/
10 *
11 * Software distributed under the License is distributed on an "AS IS" basis,
12 * WITHOUT WARRANTY OF ANY KIND, either express or implied. See the License
13 * for the specific language governing rights and limitations under the
14 * License.
15 *
16 * The Original Code is TurboPower SysTools
17 *
18 * The Initial Developer of the Original Code is
19 * TurboPower Software
20 *
21 * Portions created by the Initial Developer are Copyright (C) 1996-2002
22 * the Initial Developer. All Rights Reserved.
23 *
24 * Contributor(s):
25 *
26 * ***** END LICENSE BLOCK ***** *)
27
28 {*********************************************************}
29 {* SysTools: StUranus.pas 4.04 *}
30 {*********************************************************}
31 {* SysTools: Astronomical Routines (for Uranus) *}
32 {*********************************************************}
33
34 {$I StDefine.inc}
35
36 unit StUranus;
37
38 interface
39
40 uses
41 StAstroP;
42
43 function ComputeUranus(JD : Double) : TStEclipticalCord;
44
45 implementation
46
47 function GetLongitude(Tau, Tau2, Tau3, Tau4, Tau5 : Double) : Double;
48 var
49 L0, L1,
50 L2, L3,
51 L4, L5 : Double;
52 begin
53 L0 := 5.48129294300 * cos(0.00000000000 + 0.00000000000 * Tau)
54 + 0.09260408252 * cos(0.89106421530 + 74.78159856700 * Tau)
55 + 0.01504247826 * cos(3.62719262190 + 1.48447270830 * Tau)
56 + 0.00365981718 * cos(1.89962189070 + 73.29712585900 * Tau)
57 + 0.00272328132 * cos(3.35823710520 + 149.56319713000 * Tau)
58 + 0.00070328499 * cos(5.39254431990 + 63.73589830300 * Tau)
59 + 0.00068892609 * cos(6.09292489050 + 76.26607127600 * Tau)
60 + 0.00061998592 * cos(2.26952040470 + 2.96894541660 * Tau)
61 + 0.00061950714 * cos(2.85098907570 + 11.04570026400 * Tau)
62 + 0.00026468869 * cos(3.14152087890 + 71.81265315100 * Tau)
63 + 0.00025710505 * cos(6.11379842940 + 454.90936653000 * Tau)
64 + 0.00021078897 * cos(4.36059465140 + 148.07872443000 * Tau)
65 + 0.00017818665 * cos(1.74436982540 + 36.64856292900 * Tau)
66 + 0.00014613471 * cos(4.73732047980 + 3.93215326310 * Tau)
67 + 0.00011162535 * cos(5.82681993690 + 224.34479570000 * Tau)
68 + 0.00010997934 * cos(0.48865493179 + 138.51749687000 * Tau)
69 + 0.00009527487 * cos(2.95516893090 + 35.16409022100 * Tau)
70 + 0.00007545543 * cos(5.23626440670 + 109.94568879000 * Tau)
71 + 0.00004220170 * cos(3.23328535510 + 70.84944530400 * Tau)
72 + 0.00004051850 * cos(2.27754158720 + 151.04766984000 * Tau)
73 + 0.00003490352 * cos(5.48305567290 + 146.59425172000 * Tau)
74 + 0.00003354607 * cos(1.06549008890 + 4.45341812490 * Tau)
75 + 0.00003144093 * cos(4.75199307600 + 77.75054398400 * Tau)
76 + 0.00002926671 * cos(4.62903695490 + 9.56122755560 * Tau)
77 + 0.00002922410 * cos(5.35236743380 + 85.82729883100 * Tau)
78 + 0.00002272790 * cos(4.36600802760 + 70.32818044200 * Tau)
79 + 0.00002148599 * cos(0.60745800902 + 38.13303563800 * Tau)
80 + 0.00002051209 * cos(1.51773563460 + 0.11187458460 * Tau)
81 + 0.00001991726 * cos(4.92437290830 + 277.03499374000 * Tau)
82 + 0.00001666910 * cos(3.62744580850 + 380.12776796000 * Tau)
83 + 0.00001533223 * cos(2.58593414270 + 52.69019803900 * Tau)
84 + 0.00001376208 * cos(2.04281409050 + 65.22037101200 * Tau)
85 + 0.00001372100 * cos(4.19641615560 + 111.43016150000 * Tau)
86 + 0.00001284183 * cos(3.11346336880 + 202.25339517000 * Tau)
87 + 0.00001281641 * cos(0.54269869505 + 222.86032299000 * Tau)
88 + 0.00001244342 * cos(0.91612680579 + 2.44768055480 * Tau)
89 + 0.00001220998 * cos(0.19901396193 + 108.46121608000 * Tau)
90 + 0.00001150993 * cos(4.17898207050 + 33.67961751300 * Tau)
91 + 0.00001150416 * cos(0.93344454002 + 3.18139373770 * Tau)
92 + 0.00001090461 * cos(1.77501638910 + 12.53017297200 * Tau)
93 + 0.00001072008 * cos(0.23564502877 + 62.25142559500 * Tau)
94 + 0.00000946195 * cos(1.19249463070 + 127.47179661000 * Tau)
95 + 0.00000707875 * cos(5.18285226580 + 213.29909544000 * Tau)
96 + 0.00000653401 * cos(0.96586909116 + 78.71375183000 * Tau)
97 + 0.00000627562 * cos(0.18210181975 + 984.60033162000 * Tau)
98 + 0.00000606827 * cos(5.43209728950 + 529.69096509000 * Tau)
99 + 0.00000559370 * cos(3.35776737700 + 0.52126486180 * Tau)
100 + 0.00000524495 * cos(2.01276707000 + 299.12639427000 * Tau)
101 + 0.00000483219 * cos(2.10553990150 + 0.96320784650 * Tau)
102 + 0.00000471288 * cos(1.40664336450 + 184.72728736000 * Tau)
103 + 0.00000467211 * cos(0.41484068933 + 145.10977901000 * Tau)
104 + 0.00000433532 * cos(5.52142978260 + 183.24281465000 * Tau)
105 + 0.00000404891 * cos(5.98689011390 + 8.07675484730 * Tau)
106 + 0.00000398996 * cos(0.33810765436 + 415.55249061000 * Tau)
107 + 0.00000395614 * cos(5.87039580950 + 351.81659231000 * Tau)
108 + 0.00000378609 * cos(2.34975805010 + 56.62235130300 * Tau)
109 + 0.00000309885 * cos(5.83301304670 + 145.63104387000 * Tau)
110 + 0.00000300379 * cos(5.64353974150 + 22.09140052800 * Tau)
111 + 0.00000294172 * cos(5.83916826230 + 39.61750834600 * Tau)
112 + 0.00000251792 * cos(1.63696775580 + 221.37585029000 * Tau)
113 + 0.00000249229 * cos(4.74617120580 + 225.82926841000 * Tau)
114 + 0.00000239334 * cos(2.35045874710 + 137.03302416000 * Tau)
115 + 0.00000224097 * cos(0.51574863468 + 84.34282612300 * Tau)
116 + 0.00000222588 * cos(2.84309380330 + 0.26063243090 * Tau)
117 + 0.00000219621 * cos(1.92212987980 + 67.66805156700 * Tau)
118 + 0.00000216549 * cos(6.14211862700 + 5.93789083320 * Tau)
119 + 0.00000216480 * cos(4.77847481360 + 340.77089205000 * Tau)
120 + 0.00000207828 * cos(5.58020570040 + 68.84370773400 * Tau)
121 + 0.00000201963 * cos(1.29693040860 + 0.04818410980 * Tau)
122 + 0.00000199146 * cos(0.95634155010 + 152.53214255000 * Tau)
123 + 0.00000193652 * cos(1.88800122610 + 456.39383924000 * Tau)
124 + 0.00000192998 * cos(0.91616058506 + 453.42489382000 * Tau)
125 + 0.00000187474 * cos(1.31924326250 + 0.16005869440 * Tau)
126 + 0.00000181934 * cos(3.53624029240 + 79.23501669200 * Tau)
127 + 0.00000173145 * cos(1.53860728050 + 160.60889740000 * Tau)
128 + 0.00000171968 * cos(5.67952685530 + 219.89137758000 * Tau)
129 + 0.00000170300 * cos(3.67717520690 + 5.41662597140 * Tau)
130 + 0.00000168648 * cos(5.87874000880 + 18.15924726500 * Tau)
131 + 0.00000164588 * cos(1.42379714840 + 106.97674337000 * Tau)
132 + 0.00000162792 * cos(3.05029377670 + 112.91463421000 * Tau)
133 + 0.00000158028 * cos(0.73811997211 + 54.17467074800 * Tau)
134 + 0.00000146653 * cos(1.26300172260 + 59.80374504000 * Tau)
135 + 0.00000143058 * cos(1.29995487560 + 35.42472265200 * Tau)
136 + 0.00000139453 * cos(5.38597723400 + 32.19514480500 * Tau)
137 + 0.00000138585 * cos(4.25994786670 + 909.81873305000 * Tau)
138 + 0.00000123840 * cos(1.37359990340 + 7.11354700080 * Tau)
139 + 0.00000110163 * cos(2.02685778980 + 554.06998748000 * Tau)
140 + 0.00000109376 * cos(5.70581833290 + 77.96299230500 * Tau)
141 + 0.00000104414 * cos(5.02820888810 + 0.75075952540 * Tau)
142 + 0.00000103562 * cos(1.45770270250 + 24.37902238800 * Tau)
143 + 0.00000103277 * cos(0.68095301267 + 14.97785352700 * Tau);
144
145 L1 := 75.02543121600 * cos(0.00000000000 + 0.00000000000 * Tau)
146 + 0.00154458244 * cos(5.24201658070 + 74.78159856700 * Tau)
147 + 0.00024456413 * cos(1.71255705310 + 1.48447270830 * Tau)
148 + 0.00009257828 * cos(0.42844639064 + 11.04570026400 * Tau)
149 + 0.00008265977 * cos(1.50220035110 + 63.73589830300 * Tau)
150 + 0.00007841715 * cos(1.31983607250 + 149.56319713000 * Tau)
151 + 0.00003899105 * cos(0.46483574024 + 3.93215326310 * Tau)
152 + 0.00002283777 * cos(4.17367534000 + 76.26607127600 * Tau)
153 + 0.00001926600 * cos(0.53013080152 + 2.96894541660 * Tau)
154 + 0.00001232727 * cos(1.58634458240 + 70.84944530400 * Tau)
155 + 0.00000791206 * cos(5.43641224140 + 3.18139373770 * Tau)
156 + 0.00000766954 * cos(1.99555409580 + 73.29712585900 * Tau)
157 + 0.00000481671 * cos(2.98401996910 + 85.82729883100 * Tau)
158 + 0.00000449798 * cos(4.13826237510 + 138.51749687000 * Tau)
159 + 0.00000445600 * cos(3.72300400330 + 224.34479570000 * Tau)
160 + 0.00000426554 * cos(4.73126059390 + 71.81265315100 * Tau)
161 + 0.00000353752 * cos(2.58324496890 + 148.07872443000 * Tau)
162 + 0.00000347735 * cos(2.45372261290 + 9.56122755560 * Tau)
163 + 0.00000317084 * cos(5.57855232070 + 52.69019803900 * Tau)
164 + 0.00000205585 * cos(2.36263144250 + 2.44768055480 * Tau)
165 + 0.00000189068 * cos(4.20242881380 + 56.62235130300 * Tau)
166 + 0.00000183762 * cos(0.28371004654 + 151.04766984000 * Tau)
167 + 0.00000179920 * cos(5.68367730920 + 12.53017297200 * Tau)
168 + 0.00000171084 * cos(3.00060075290 + 78.71375183000 * Tau)
169 + 0.00000158029 * cos(2.90931969500 + 0.96320784650 * Tau)
170 + 0.00000154670 * cos(5.59083925610 + 4.45341812490 * Tau)
171 + 0.00000153515 * cos(4.65186885940 + 35.16409022100 * Tau)
172 + 0.00000151984 * cos(2.94217326890 + 77.75054398400 * Tau)
173 + 0.00000143464 * cos(2.59049246730 + 62.25142559500 * Tau)
174 + 0.00000121452 * cos(4.14839204920 + 127.47179661000 * Tau)
175 + 0.00000115546 * cos(3.73224603790 + 65.22037101200 * Tau)
176 + 0.00000102022 * cos(4.18754517990 + 145.63104387000 * Tau)
177 + 0.00000101718 * cos(6.03385875010 + 0.11187458460 * Tau)
178 + 0.00000088202 * cos(3.99035787990 + 18.15924726500 * Tau)
179 + 0.00000087549 * cos(6.15520787580 + 202.25339517000 * Tau)
180 + 0.00000080530 * cos(2.64124743930 + 22.09140052800 * Tau)
181 + 0.00000072047 * cos(6.04545933580 + 70.32818044200 * Tau)
182 + 0.00000068570 * cos(4.05071895260 + 77.96299230500 * Tau)
183 + 0.00000059173 * cos(3.70413919080 + 67.66805156700 * Tau)
184 + 0.00000047267 * cos(3.54312460520 + 351.81659231000 * Tau)
185 + 0.00000044339 * cos(5.90865821910 + 7.11354700080 * Tau)
186 + 0.00000042534 * cos(5.72357370900 + 5.41662597140 * Tau)
187 + 0.00000038544 * cos(4.91519003850 + 222.86032299000 * Tau)
188 + 0.00000036116 * cos(5.89964278800 + 33.67961751300 * Tau)
189 + 0.00000035605 * cos(3.29197259180 + 8.07675484730 * Tau)
190 + 0.00000035524 * cos(3.32784616140 + 71.60020483000 * Tau)
191 + 0.00000034996 * cos(5.08034112150 + 38.13303563800 * Tau)
192 + 0.00000031454 * cos(5.62015632300 + 984.60033162000 * Tau)
193 + 0.00000030811 * cos(5.49591403860 + 59.80374504000 * Tau)
194 + 0.00000030608 * cos(5.46414592600 + 160.60889740000 * Tau)
195 + 0.00000029866 * cos(1.65980844670 + 447.79581953000 * Tau)
196 + 0.00000029206 * cos(1.14722640420 + 462.02291353000 * Tau)
197 + 0.00000028947 * cos(4.51867390410 + 84.34282612300 * Tau)
198 + 0.00000026627 * cos(5.54127301040 + 131.40394987000 * Tau)
199 + 0.00000026605 * cos(6.14640604130 + 299.12639427000 * Tau)
200 + 0.00000025753 * cos(4.99362028420 + 137.03302416000 * Tau)
201 + 0.00000025373 * cos(5.73584678600 + 380.12776796000 * Tau);
202
203 L2 := 0.00053033277 * cos(0.00000000000 + 0.00000000000 * Tau)
204 + 0.00002357636 * cos(2.26014661700 + 74.78159856700 * Tau)
205 + 0.00000769129 * cos(4.52561041820 + 11.04570026400 * Tau)
206 + 0.00000551533 * cos(3.25814281020 + 63.73589830300 * Tau)
207 + 0.00000541532 * cos(2.27573907420 + 3.93215326310 * Tau)
208 + 0.00000529473 * cos(4.92348433830 + 1.48447270830 * Tau)
209 + 0.00000257521 * cos(3.69059216860 + 3.18139373770 * Tau)
210 + 0.00000238835 * cos(5.85806638400 + 149.56319713000 * Tau)
211 + 0.00000181904 * cos(6.21763603410 + 70.84944530400 * Tau)
212 + 0.00000053504 * cos(1.44225240950 + 76.26607127600 * Tau)
213 + 0.00000049401 * cos(6.03101301720 + 56.62235130300 * Tau)
214 + 0.00000044753 * cos(3.90904910520 + 2.44768055480 * Tau)
215 + 0.00000044530 * cos(0.81152639478 + 85.82729883100 * Tau)
216 + 0.00000038222 * cos(1.78467827780 + 52.69019803900 * Tau)
217 + 0.00000037403 * cos(4.46228598030 + 2.96894541660 * Tau)
218 + 0.00000033029 * cos(0.86388149962 + 9.56122755560 * Tau)
219 + 0.00000029423 * cos(5.09818697710 + 73.29712585900 * Tau)
220 + 0.00000024292 * cos(2.10702559050 + 18.15924726500 * Tau)
221 + 0.00000022491 * cos(5.99320728690 + 138.51749687000 * Tau)
222 + 0.00000022135 * cos(4.81730808580 + 78.71375183000 * Tau)
223 + 0.00000021392 * cos(2.39880709310 + 77.96299230500 * Tau)
224 + 0.00000020578 * cos(2.16918786540 + 224.34479570000 * Tau)
225 + 0.00000017226 * cos(2.53537183200 + 145.63104387000 * Tau)
226 + 0.00000016777 * cos(3.46631344090 + 12.53017297200 * Tau)
227 + 0.00000012012 * cos(0.01941361902 + 22.09140052800 * Tau)
228 + 0.00000011010 * cos(0.08496274370 + 127.47179661000 * Tau)
229 + 0.00000010476 * cos(5.16453084070 + 71.60020483000 * Tau)
230 + 0.00000010466 * cos(4.45556032590 + 62.25142559500 * Tau)
231 + 0.00000008668 * cos(4.25550086980 + 7.11354700080 * Tau)
232 + 0.00000008387 * cos(5.50115930050 + 67.66805156700 * Tau)
233 + 0.00000007160 * cos(1.24903906390 + 5.41662597140 * Tau)
234 + 0.00000006109 * cos(3.36320161280 + 447.79581953000 * Tau)
235 + 0.00000006087 * cos(5.44611674380 + 65.22037101200 * Tau)
236 + 0.00000006013 * cos(4.51836836350 + 151.04766984000 * Tau)
237 + 0.00000006003 * cos(5.72500086740 + 462.02291353000 * Tau);
238
239 L3 := 0.00000120936 * cos(0.02418789918 + 74.78159856700 * Tau)
240 + 0.00000068064 * cos(4.12084267730 + 3.93215326310 * Tau)
241 + 0.00000052828 * cos(2.38964061260 + 11.04570026400 * Tau)
242 + 0.00000045806 * cos(0.00000000000 + 0.00000000000 * Tau)
243 + 0.00000045300 * cos(2.04423798410 + 3.18139373770 * Tau)
244 + 0.00000043754 * cos(2.95965039730 + 1.48447270830 * Tau)
245 + 0.00000024969 * cos(4.88741307920 + 63.73589830300 * Tau)
246 + 0.00000021061 * cos(4.54511486860 + 70.84944530400 * Tau)
247 + 0.00000019897 * cos(2.31320314140 + 149.56319713000 * Tau)
248 + 0.00000008901 * cos(1.57548871760 + 56.62235130300 * Tau)
249 + 0.00000004271 * cos(0.22777319552 + 18.15924726500 * Tau)
250 + 0.00000003613 * cos(5.39244611310 + 76.26607127600 * Tau)
251 + 0.00000003572 * cos(0.95052448578 + 77.96299230500 * Tau)
252 + 0.00000003488 * cos(4.97622811780 + 85.82729883100 * Tau)
253 + 0.00000003479 * cos(4.12969359980 + 52.69019803900 * Tau)
254 + 0.00000002696 * cos(0.37287796344 + 78.71375183000 * Tau)
255 + 0.00000002328 * cos(0.85770961794 + 145.63104387000 * Tau)
256 + 0.00000002156 * cos(5.65647821520 + 9.56122755560 * Tau);
257
258 L4 := 0.00000113855 * cos(3.14159265360 + 0.00000000000 * Tau)
259 + 0.00000005599 * cos(4.57882424420 + 74.78159856700 * Tau)
260 + 0.00000003203 * cos(0.34623003207 + 11.04570026400 * Tau)
261 + 0.00000001217 * cos(3.42199121830 + 56.62235130300 * Tau);
262
263 L5 := 0.00000000000;
264 Result := (L0 + L1*Tau + L2*Tau2 + L3*Tau3 + L4*Tau4 + L5*Tau5);
265 end;
266
267 {---------------------------------------------------------------------------}
268
269 function GetLatitude(Tau, Tau2, Tau3, Tau4, Tau5 : Double) : Double;
270 var
271 B0, B1,
272 B2, B3,
273 B4, B5 : Double;
274 begin
275 B0 := 0.01346277639 * cos(2.61877810550 + 74.78159856700 * Tau)
276 + 0.00062341405 * cos(5.08111175860 + 149.56319713000 * Tau)
277 + 0.00061601203 * cos(3.14159265360 + 0.00000000000 * Tau)
278 + 0.00009963744 * cos(1.61603876360 + 76.26607127600 * Tau)
279 + 0.00009926151 * cos(0.57630387917 + 73.29712585900 * Tau)
280 + 0.00003259455 * cos(1.26119385960 + 224.34479570000 * Tau)
281 + 0.00002972318 * cos(2.24367035540 + 1.48447270830 * Tau)
282 + 0.00002010257 * cos(6.05550401090 + 148.07872443000 * Tau)
283 + 0.00001522172 * cos(0.27960386377 + 63.73589830300 * Tau)
284 + 0.00000924055 * cos(4.03822927850 + 151.04766984000 * Tau)
285 + 0.00000760624 * cos(6.14000431920 + 71.81265315100 * Tau)
286 + 0.00000522309 * cos(3.32085194770 + 138.51749687000 * Tau)
287 + 0.00000462630 * cos(0.74256727574 + 85.82729883100 * Tau)
288 + 0.00000436843 * cos(3.38082524320 + 529.69096509000 * Tau)
289 + 0.00000434625 * cos(0.34065281858 + 77.75054398400 * Tau)
290 + 0.00000430668 * cos(3.55445034850 + 213.29909544000 * Tau)
291 + 0.00000420265 * cos(5.21279984790 + 11.04570026400 * Tau)
292 + 0.00000244698 * cos(0.78795150326 + 2.96894541660 * Tau)
293 + 0.00000232649 * cos(2.25716421380 + 222.86032299000 * Tau)
294 + 0.00000215838 * cos(1.59121704940 + 38.13303563800 * Tau)
295 + 0.00000179935 * cos(3.72487952670 + 299.12639427000 * Tau)
296 + 0.00000174895 * cos(1.23550262210 + 146.59425172000 * Tau)
297 + 0.00000173667 * cos(1.93654269130 + 380.12776796000 * Tau)
298 + 0.00000160368 * cos(5.33635436460 + 111.43016150000 * Tau)
299 + 0.00000144064 * cos(5.96239326410 + 35.16409022100 * Tau)
300 + 0.00000116363 * cos(5.73877190010 + 70.84944530400 * Tau)
301 + 0.00000106441 * cos(0.94103112994 + 70.32818044200 * Tau)
302 + 0.00000102049 * cos(2.61876256510 + 78.71375183000 * Tau);
303
304 B1 := 0.00206366162 * cos(4.12394311410 + 74.78159856700 * Tau)
305 + 0.00008563230 * cos(0.33819986165 + 149.56319713000 * Tau)
306 + 0.00001725703 * cos(2.12193159900 + 73.29712585900 * Tau)
307 + 0.00001374449 * cos(0.00000000000 + 0.00000000000 * Tau)
308 + 0.00001368860 * cos(3.06861722050 + 76.26607127600 * Tau)
309 + 0.00000450639 * cos(3.77656180980 + 1.48447270830 * Tau)
310 + 0.00000399847 * cos(2.84767037790 + 224.34479570000 * Tau)
311 + 0.00000307214 * cos(1.25456766740 + 148.07872443000 * Tau)
312 + 0.00000154336 * cos(3.78575467750 + 63.73589830300 * Tau)
313 + 0.00000112432 * cos(5.57299891500 + 151.04766984000 * Tau)
314 + 0.00000110888 * cos(5.32888676460 + 138.51749687000 * Tau)
315 + 0.00000083493 * cos(3.59152795560 + 71.81265315100 * Tau)
316 + 0.00000055573 * cos(3.40135416350 + 85.82729883100 * Tau)
317 + 0.00000053690 * cos(1.70455769940 + 77.75054398400 * Tau)
318 + 0.00000041912 * cos(1.21476607430 + 11.04570026400 * Tau)
319 + 0.00000041377 * cos(4.45476669140 + 78.71375183000 * Tau)
320 + 0.00000031959 * cos(3.77446207750 + 222.86032299000 * Tau)
321 + 0.00000030297 * cos(2.56371683640 + 2.96894541660 * Tau)
322 + 0.00000026977 * cos(5.33695500290 + 213.29909544000 * Tau)
323 + 0.00000026222 * cos(0.41620628369 + 380.12776796000 * Tau);
324
325 B2 := 0.00009211656 * cos(5.80044305790 + 74.78159856700 * Tau)
326 + 0.00000556926 * cos(0.00000000000 + 0.00000000000 * Tau)
327 + 0.00000286265 * cos(2.17729776350 + 149.56319713000 * Tau)
328 + 0.00000094969 * cos(3.84237569810 + 73.29712585900 * Tau)
329 + 0.00000045419 * cos(4.87822046060 + 76.26607127600 * Tau)
330 + 0.00000020107 * cos(5.46264485370 + 1.48447270830 * Tau)
331 + 0.00000014793 * cos(0.87983715652 + 138.51749687000 * Tau)
332 + 0.00000014261 * cos(2.84517742690 + 148.07872443000 * Tau)
333 + 0.00000013963 * cos(5.07234043990 + 63.73589830300 * Tau)
334 + 0.00000010122 * cos(5.00290894860 + 224.34479570000 * Tau)
335 + 0.00000008299 * cos(6.26655615200 + 78.71375183000 * Tau);
336
337 B3 := 0.00000267832 * cos(1.25097888290 + 74.78159856700 * Tau)
338 + 0.00000011048 * cos(3.14159265360 + 0.00000000000 * Tau)
339 + 0.00000006154 * cos(4.00663614490 + 149.56319713000 * Tau)
340 + 0.00000003361 * cos(5.77804694940 + 73.29712585900 * Tau);
341
342 B4 := 0.00000005719 * cos(2.85499529310 + 74.78159856700 * Tau);
343
344 B5 := 0.00000000000;
345 Result := (B0 + B1*Tau + B2*Tau2 + B3*Tau3 + B4*Tau4 + B5*Tau5);
346 end;
347
348 {---------------------------------------------------------------------------}
349
350 function GetRadiusVector(Tau, Tau2, Tau3, Tau4, Tau5 : Double) : Double;
351 var
352 R0, R1,
353 R2, R3,
354 R4, R5 : Double;
355 begin
356 R0 := 19.21264847900 * cos(0.00000000000 + 0.00000000000 * Tau)
357 + 0.88784984055 * cos(5.60377526990 + 74.78159856700 * Tau)
358 + 0.03440835545 * cos(0.32836098991 + 73.29712585900 * Tau)
359 + 0.02055653495 * cos(1.78295170030 + 149.56319713000 * Tau)
360 + 0.00649321851 * cos(4.52247298120 + 76.26607127600 * Tau)
361 + 0.00602248144 * cos(3.86003820460 + 63.73589830300 * Tau)
362 + 0.00496404171 * cos(1.40139934720 + 454.90936653000 * Tau)
363 + 0.00338525522 * cos(1.58002682950 + 138.51749687000 * Tau)
364 + 0.00243508222 * cos(1.57086595070 + 71.81265315100 * Tau)
365 + 0.00190521915 * cos(1.99809364500 + 1.48447270830 * Tau)
366 + 0.00161858251 * cos(2.79137863470 + 148.07872443000 * Tau)
367 + 0.00143705902 * cos(1.38368574480 + 11.04570026400 * Tau)
368 + 0.00093192359 * cos(0.17437193645 + 36.64856292900 * Tau)
369 + 0.00089805842 * cos(3.66105366330 + 109.94568879000 * Tau)
370 + 0.00071424265 * cos(4.24509327400 + 224.34479570000 * Tau)
371 + 0.00046677322 * cos(1.39976563940 + 35.16409022100 * Tau)
372 + 0.00039025681 * cos(3.36234710690 + 277.03499374000 * Tau)
373 + 0.00039009624 * cos(1.66971128870 + 70.84944530400 * Tau)
374 + 0.00036755160 * cos(3.88648934740 + 146.59425172000 * Tau)
375 + 0.00030348875 * cos(0.70100446346 + 151.04766984000 * Tau)
376 + 0.00029156264 * cos(3.18056174560 + 77.75054398400 * Tau)
377 + 0.00025785805 * cos(3.78537741500 + 85.82729883100 * Tau)
378 + 0.00025620360 * cos(5.25656292800 + 380.12776796000 * Tau)
379 + 0.00022637152 * cos(0.72519137745 + 529.69096509000 * Tau)
380 + 0.00020473163 * cos(2.79639811630 + 70.32818044200 * Tau)
381 + 0.00020471584 * cos(1.55588961500 + 202.25339517000 * Tau)
382 + 0.00017900561 * cos(0.55455488605 + 2.96894541660 * Tau)
383 + 0.00015502809 * cos(5.35405037600 + 38.13303563800 * Tau)
384 + 0.00014701566 * cos(4.90434406650 + 108.46121608000 * Tau)
385 + 0.00012896507 * cos(2.62154018240 + 111.43016150000 * Tau)
386 + 0.00012328151 * cos(5.96039150920 + 127.47179661000 * Tau)
387 + 0.00011959355 * cos(1.75044072170 + 984.60033162000 * Tau)
388 + 0.00011852996 * cos(0.99342814582 + 52.69019803900 * Tau)
389 + 0.00011696085 * cos(3.29825599110 + 3.93215326310 * Tau)
390 + 0.00011494701 * cos(0.43774027872 + 65.22037101200 * Tau)
391 + 0.00010792699 * cos(1.42104858470 + 213.29909544000 * Tau)
392 + 0.00009111446 * cos(4.99638600050 + 62.25142559500 * Tau)
393 + 0.00008420550 * cos(5.25350716620 + 222.86032299000 * Tau)
394 + 0.00008402147 * cos(5.03877516490 + 415.55249061000 * Tau)
395 + 0.00007449125 * cos(0.79491905956 + 351.81659231000 * Tau)
396 + 0.00007329454 * cos(3.97277527840 + 183.24281465000 * Tau)
397 + 0.00006046370 * cos(5.67960948360 + 78.71375183000 * Tau)
398 + 0.00005524133 * cos(3.11499484160 + 9.56122755560 * Tau)
399 + 0.00005444878 * cos(5.10575635360 + 145.10977901000 * Tau)
400 + 0.00005238103 * cos(2.62960141800 + 33.67961751300 * Tau)
401 + 0.00004079167 * cos(3.22064788670 + 340.77089205000 * Tau)
402 + 0.00003919476 * cos(4.25015288870 + 39.61750834600 * Tau)
403 + 0.00003801606 * cos(6.10985558500 + 184.72728736000 * Tau)
404 + 0.00003781219 * cos(3.45840272870 + 456.39383924000 * Tau)
405 + 0.00003686787 * cos(2.48718116540 + 453.42489382000 * Tau)
406 + 0.00003101743 * cos(4.14031063900 + 219.89137758000 * Tau)
407 + 0.00002962641 * cos(0.82977991995 + 56.62235130300 * Tau)
408 + 0.00002942239 * cos(0.42393808854 + 299.12639427000 * Tau)
409 + 0.00002940492 * cos(2.14637460320 + 137.03302416000 * Tau)
410 + 0.00002937799 * cos(3.67657450930 + 140.00196958000 * Tau)
411 + 0.00002865128 * cos(0.30996903761 + 12.53017297200 * Tau)
412 + 0.00002538032 * cos(4.85457831990 + 131.40394987000 * Tau)
413 + 0.00002363550 * cos(0.44253328372 + 554.06998748000 * Tau)
414 + 0.00002182572 * cos(2.94040431640 + 305.34616939000 * Tau);
415
416 R1 := 0.01479896370 * cos(3.67205705320 + 74.78159856700 * Tau)
417 + 0.00071212085 * cos(6.22601006670 + 63.73589830300 * Tau)
418 + 0.00068626972 * cos(6.13411265050 + 149.56319713000 * Tau)
419 + 0.00024059649 * cos(3.14159265360 + 0.00000000000 * Tau)
420 + 0.00021468152 * cos(2.60176704270 + 76.26607127600 * Tau)
421 + 0.00020857262 * cos(5.24625494220 + 11.04570026400 * Tau)
422 + 0.00011405346 * cos(0.01848461561 + 70.84944530400 * Tau)
423 + 0.00007496775 * cos(0.42360033283 + 73.29712585900 * Tau)
424 + 0.00004243800 * cos(1.41692350370 + 85.82729883100 * Tau)
425 + 0.00003926694 * cos(3.15513991320 + 71.81265315100 * Tau)
426 + 0.00003578446 * cos(2.31160668310 + 224.34479570000 * Tau)
427 + 0.00003505936 * cos(2.58354048850 + 138.51749687000 * Tau)
428 + 0.00003228835 * cos(5.25499602900 + 3.93215326310 * Tau)
429 + 0.00003060010 * cos(0.15321893225 + 1.48447270830 * Tau)
430 + 0.00002564251 * cos(0.98076846352 + 148.07872443000 * Tau)
431 + 0.00002429445 * cos(3.99440122470 + 52.69019803900 * Tau)
432 + 0.00001644719 * cos(2.65349313120 + 127.47179661000 * Tau)
433 + 0.00001583766 * cos(1.43045619200 + 78.71375183000 * Tau)
434 + 0.00001508028 * cos(5.05996325430 + 151.04766984000 * Tau)
435 + 0.00001489525 * cos(2.67559167320 + 56.62235130300 * Tau)
436 + 0.00001413112 * cos(4.57461892060 + 202.25339517000 * Tau)
437 + 0.00001403237 * cos(1.36985349740 + 77.75054398400 * Tau)
438 + 0.00001228220 * cos(1.04703640150 + 62.25142559500 * Tau)
439 + 0.00001032731 * cos(0.26459059027 + 131.40394987000 * Tau)
440 + 0.00000992085 * cos(2.17168865910 + 65.22037101200 * Tau)
441 + 0.00000861867 * cos(5.05530802220 + 351.81659231000 * Tau)
442 + 0.00000744445 * cos(3.07640148940 + 35.16409022100 * Tau)
443 + 0.00000687470 * cos(2.49912565670 + 77.96299230500 * Tau)
444 + 0.00000646851 * cos(4.47290422910 + 70.32818044200 * Tau)
445 + 0.00000623602 * cos(0.86253073820 + 9.56122755560 * Tau)
446 + 0.00000604362 * cos(0.90717667985 + 984.60033162000 * Tau)
447 + 0.00000574710 * cos(3.23070708460 + 447.79581953000 * Tau)
448 + 0.00000561839 * cos(2.71778158980 + 462.02291353000 * Tau)
449 + 0.00000530364 * cos(5.91655309050 + 213.29909544000 * Tau)
450 + 0.00000527794 * cos(5.15136007080 + 2.96894541660 * Tau);
451
452 R2 := 0.00022439904 * cos(0.69953118760 + 74.78159856700 * Tau)
453 + 0.00004727037 * cos(1.69901641490 + 63.73589830300 * Tau)
454 + 0.00001681903 * cos(4.64833551730 + 70.84944530400 * Tau)
455 + 0.00001649559 * cos(3.09660078980 + 11.04570026400 * Tau)
456 + 0.00001433755 * cos(3.52119917950 + 149.56319713000 * Tau)
457 + 0.00000770188 * cos(0.00000000000 + 0.00000000000 * Tau)
458 + 0.00000500429 * cos(6.17229032220 + 76.26607127600 * Tau)
459 + 0.00000461009 * cos(0.76676632849 + 3.93215326310 * Tau)
460 + 0.00000390371 * cos(4.49605283500 + 56.62235130300 * Tau)
461 + 0.00000389945 * cos(5.52673426380 + 85.82729883100 * Tau)
462 + 0.00000292097 * cos(0.20389012095 + 52.69019803900 * Tau)
463 + 0.00000286579 * cos(3.53357683270 + 73.29712585900 * Tau)
464 + 0.00000272898 * cos(3.84707823650 + 138.51749687000 * Tau)
465 + 0.00000219674 * cos(1.96418942890 + 131.40394987000 * Tau)
466 + 0.00000215788 * cos(0.84812474187 + 77.96299230500 * Tau)
467 + 0.00000205449 * cos(3.24758017120 + 78.71375183000 * Tau)
468 + 0.00000148554 * cos(4.89840863840 + 127.47179661000 * Tau)
469 + 0.00000128834 * cos(2.08146849520 + 3.18139373770 * Tau);
470
471 R3 := 0.00001164382 * cos(4.73453291600 + 74.78159856700 * Tau)
472 + 0.00000212367 * cos(3.34255735000 + 63.73589830300 * Tau)
473 + 0.00000196408 * cos(2.98004616320 + 70.84944530400 * Tau)
474 + 0.00000104527 * cos(0.95807937648 + 11.04570026400 * Tau)
475 + 0.00000072540 * cos(0.99701907912 + 149.56319713000 * Tau)
476 + 0.00000071681 * cos(0.02528455665 + 56.62235130300 * Tau)
477 + 0.00000054875 * cos(2.59436811270 + 3.93215326310 * Tau)
478 + 0.00000036377 * cos(5.65035573020 + 77.96299230500 * Tau)
479 + 0.00000034029 * cos(3.81553325640 + 76.26607127600 * Tau)
480 + 0.00000032081 * cos(3.59825177840 + 131.40394987000 * Tau);
481
482 R4 := 0.00000052996 * cos(3.00838033090 + 74.78159856700 * Tau)
483 + 0.00000009887 * cos(1.91399083600 + 56.62235130300 * Tau);
484
485 R5 := 0.00000000000;
486 Result := (R0 + R1*Tau + R2*Tau2 + R3*Tau3 + R4*Tau4 + R5*Tau5);
487 end;
488
489 {---------------------------------------------------------------------------}
490
491 function ComputeUranus(JD : Double) : TStEclipticalCord;
492 var
493 Tau,
494 Tau2,
495 Tau3,
496 Tau4,
497 Tau5 : Double;
498 begin
499 Tau := (JD - 2451545.0) / 365250.0;
500 Tau2 := sqr(Tau);
501 Tau3 := Tau * Tau2;
502 Tau4 := sqr(Tau2);
503 Tau5 := Tau2 * Tau3;
504
505 Result.L0 := GetLongitude(Tau, Tau2, Tau3, Tau4, Tau5);
506 Result.B0 := GetLatitude(Tau, Tau2, Tau3, Tau4, Tau5);
507 Result.R0 := GetRadiusVector(Tau, Tau2, Tau3, Tau4, Tau5);
508 end;
509
510
511 end.
  ViewVC Help
Powered by ViewVC 1.1.20