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

Annotation of /dao/DelphiScanner/Components/tpsystools_4.04/source/StJup.pas

Parent Directory Parent Directory | Revision Log Revision Log

Revision 2671 - (hide annotations) (download)
Tue Aug 25 18:15:15 2015 UTC (8 years, 10 months ago) by torben
File size: 38453 byte(s) 
Added tpsystools component
1 torben 2671 // 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: StJup.pas 4.04 *}
30     {*********************************************************}
31     {* SysTools: Astronomical Routines (for Jupiter) *}
32     {*********************************************************}
33    
34     {$I StDefine.inc}
35    
36     unit StJup;
37    
38     interface
39    
40     uses
41     StAstroP;
42    
43     function ComputeJupiter(JD : Double) : TStEclipticalCord;
44    
45    
46     implementation
47    
48     {--------------------------------------------------------------------------}
49    
50     function GetLongitude(Tau, Tau2, Tau3, Tau4, Tau5 : Double) : Double;
51     var
52     L0, L1,
53     L2, L3,
54     L4, L5 : Double;
55     begin
56     L0 := 0.59954691495 * cos(0.00000000000 + 0.00000000000 * Tau)
57     + 0.09695898711 * cos(5.06191793110 + 529.69096509000 * Tau)
58     + 0.00573610145 * cos(1.44406205980 + 7.11354700080 * Tau)
59     + 0.00306389180 * cos(5.41734729980 + 1059.38193020000 * Tau)
60     + 0.00097178280 * cos(4.14264708820 + 632.78373931000 * Tau)
61     + 0.00072903096 * cos(3.64042909250 + 522.57741809000 * Tau)
62     + 0.00064263986 * cos(3.41145185200 + 103.09277422000 * Tau)
63     + 0.00039806051 * cos(2.29376744860 + 419.48464388000 * Tau)
64     + 0.00038857780 * cos(1.27231724860 + 316.39186966000 * Tau)
65     + 0.00027964622 * cos(1.78454589490 + 536.80451210000 * Tau)
66     + 0.00013589738 * cos(5.77481031590 + 1589.07289530000 * Tau)
67     + 0.00008768686 * cos(3.63000324420 + 949.17560897000 * Tau)
68     + 0.00008246362 * cos(3.58227961650 + 206.18554844000 * Tau)
69     + 0.00007368057 * cos(5.08101125610 + 735.87651353000 * Tau)
70     + 0.00006263171 * cos(0.02497643742 + 213.29909544000 * Tau)
71     + 0.00006114050 * cos(4.51319531670 + 1162.47470440000 * Tau)
72     + 0.00005305457 * cos(4.18625053490 + 1052.26838320000 * Tau)
73     + 0.00005305283 * cos(1.30671236850 + 14.22709400200 * Tau)
74     + 0.00004905419 * cos(1.32084631680 + 110.20632122000 * Tau)
75     + 0.00004647249 * cos(4.69958109500 + 3.93215326310 * Tau)
76     + 0.00003045009 * cos(4.31675960320 + 426.59819088000 * Tau)
77     + 0.00002610001 * cos(1.56667594850 + 846.08283475000 * Tau)
78     + 0.00002028191 * cos(1.06376547380 + 3.18139373770 * Tau)
79     + 0.00001920959 * cos(0.97168928755 + 639.89728631000 * Tau)
80     + 0.00001764768 * cos(2.14148077770 + 1066.49547720000 * Tau)
81     + 0.00001722983 * cos(3.88036008870 + 1265.56747860000 * Tau)
82     + 0.00001633217 * cos(3.58201089760 + 515.46387109000 * Tau)
83     + 0.00001431997 * cos(4.29683690270 + 625.67019231000 * Tau)
84     + 0.00000973278 * cos(4.09764957060 + 95.97922721800 * Tau)
85     + 0.00000884439 * cos(2.43701426120 + 412.37109687000 * Tau)
86     + 0.00000732875 * cos(6.08534113240 + 838.96928775000 * Tau)
87     + 0.00000731072 * cos(3.80591233960 + 1581.95934830000 * Tau)
88     + 0.00000709190 * cos(1.29272573660 + 742.99006053000 * Tau)
89     + 0.00000691928 * cos(6.13368222940 + 2118.76386040000 * Tau)
90     + 0.00000614464 * cos(4.10853496760 + 1478.86657410000 * Tau)
91     + 0.00000581902 * cos(4.53967717550 + 309.27832266000 * Tau)
92     + 0.00000495224 * cos(3.75567461380 + 323.50541666000 * Tau)
93     + 0.00000440854 * cos(2.95818460940 + 454.90936653000 * Tau)
94     + 0.00000417266 * cos(1.03554430160 + 2.44768055480 * Tau)
95     + 0.00000389864 * cos(4.89716105850 + 1692.16566950000 * Tau)
96     + 0.00000375657 * cos(4.70299124830 + 1368.66025280000 * Tau)
97     + 0.00000341006 * cos(5.71452525780 + 533.62311836000 * Tau)
98     + 0.00000330458 * cos(4.74049819490 + 0.04818410980 * Tau)
99     + 0.00000261540 * cos(1.87652461030 + 0.96320784650 * Tau)
100     + 0.00000261009 * cos(0.82047246448 + 380.12776796000 * Tau)
101     + 0.00000256568 * cos(3.72410724160 + 199.07200144000 * Tau)
102     + 0.00000244170 * cos(5.22020878900 + 728.76296653000 * Tau)
103     + 0.00000235141 * cos(1.22693908120 + 909.81873305000 * Tau)
104     + 0.00000220382 * cos(1.65115016000 + 543.91805910000 * Tau)
105     + 0.00000207327 * cos(1.85461666590 + 525.75881183000 * Tau)
106     + 0.00000201996 * cos(1.80684574190 + 1375.77379980000 * Tau)
107     + 0.00000197046 * cos(5.29252149020 + 1155.36115740000 * Tau)
108     + 0.00000175191 * cos(3.72966554760 + 942.06206197000 * Tau)
109     + 0.00000175184 * cos(3.22634903430 + 1898.35121790000 * Tau)
110     + 0.00000174809 * cos(5.90973505280 + 956.28915597000 * Tau)
111     + 0.00000157909 * cos(4.36483921770 + 1795.25844370000 * Tau)
112     + 0.00000150502 * cos(3.90625022620 + 74.78159856700 * Tau)
113     + 0.00000149368 * cos(4.37745104270 + 1685.05212250000 * Tau)
114     + 0.00000141445 * cos(3.13568357860 + 491.55792946000 * Tau)
115     + 0.00000137871 * cos(1.31797920780 + 1169.58825140000 * Tau)
116     + 0.00000130531 * cos(4.16867945490 + 1045.15483620000 * Tau)
117     + 0.00000117495 * cos(2.50022140890 + 1596.18644230000 * Tau)
118     + 0.00000116757 * cos(3.38920921040 + 0.52126486180 * Tau)
119     + 0.00000105895 * cos(4.55439798240 + 526.50957136000 * Tau);
120    
121     L1 := 529.93480758000 * cos(0.00000000000 + 0.00000000000 * Tau)
122     + 0.00489741194 * cos(4.22066689930 + 529.69096509000 * Tau)
123     + 0.00228918538 * cos(6.02647464020 + 7.11354700080 * Tau)
124     + 0.00027655380 * cos(4.57265956820 + 1059.38193020000 * Tau)
125     + 0.00020720943 * cos(5.45938936290 + 522.57741809000 * Tau)
126     + 0.00012105732 * cos(0.16985765041 + 536.80451210000 * Tau)
127     + 0.00006068051 * cos(4.42419502010 + 103.09277422000 * Tau)
128     + 0.00005433924 * cos(3.98478382570 + 419.48464388000 * Tau)
129     + 0.00004237795 * cos(5.89009351270 + 14.22709400200 * Tau)
130     + 0.00002211854 * cos(5.26771446620 + 206.18554844000 * Tau)
131     + 0.00001745919 * cos(4.92669378490 + 1589.07289530000 * Tau)
132     + 0.00001295769 * cos(5.55132765090 + 3.18139373770 * Tau)
133     + 0.00001173129 * cos(5.85647304350 + 1052.26838320000 * Tau)
134     + 0.00001163411 * cos(0.51450895328 + 3.93215326310 * Tau)
135     + 0.00001098735 * cos(5.30704981590 + 515.46387109000 * Tau)
136     + 0.00001007216 * cos(0.46478398551 + 735.87651353000 * Tau)
137     + 0.00001003574 * cos(3.15040301820 + 426.59819088000 * Tau)
138     + 0.00000847678 * cos(5.75805850450 + 110.20632122000 * Tau)
139     + 0.00000827329 * cos(4.80312015730 + 213.29909544000 * Tau)
140     + 0.00000816397 * cos(0.58643054886 + 1066.49547720000 * Tau)
141     + 0.00000725447 * cos(5.51827471470 + 639.89728631000 * Tau)
142     + 0.00000567845 * cos(5.98867049450 + 625.67019231000 * Tau)
143     + 0.00000474181 * cos(4.13245269170 + 412.37109687000 * Tau)
144     + 0.00000412930 * cos(5.73652891260 + 95.97922721800 * Tau)
145     + 0.00000345249 * cos(4.24159565410 + 632.78373931000 * Tau)
146     + 0.00000335817 * cos(3.73248749050 + 1162.47470440000 * Tau)
147     + 0.00000234340 * cos(4.03469970330 + 949.17560897000 * Tau)
148     + 0.00000234066 * cos(6.24302226650 + 309.27832266000 * Tau)
149     + 0.00000198525 * cos(1.50458442830 + 838.96928775000 * Tau)
150     + 0.00000194784 * cos(2.21879010910 + 323.50541666000 * Tau)
151     + 0.00000186899 * cos(6.08620565910 + 742.99006053000 * Tau)
152     + 0.00000183938 * cos(6.27963588820 + 543.91805910000 * Tau)
153     + 0.00000171380 * cos(5.41655983840 + 199.07200144000 * Tau)
154     + 0.00000130771 * cos(0.62643377351 + 728.76296653000 * Tau)
155     + 0.00000115393 * cos(0.68019050174 + 846.08283475000 * Tau)
156     + 0.00000115047 * cos(5.28641699140 + 2118.76386040000 * Tau)
157     + 0.00000107575 * cos(4.49282760120 + 956.28915597000 * Tau)
158     + 0.00000079686 * cos(5.82412400270 + 1045.15483620000 * Tau)
159     + 0.00000071643 * cos(5.34162650320 + 942.06206197000 * Tau)
160     + 0.00000069618 * cos(5.97263450280 + 532.87235883000 * Tau)
161     + 0.00000066824 * cos(5.73365126530 + 21.34064100200 * Tau)
162     + 0.00000065635 * cos(0.12924191430 + 526.50957136000 * Tau)
163     + 0.00000064850 * cos(6.08803490290 + 1581.95934830000 * Tau)
164     + 0.00000058509 * cos(0.58626971028 + 1155.36115740000 * Tau)
165     + 0.00000057939 * cos(0.99453087342 + 1596.18644230000 * Tau)
166     + 0.00000057368 * cos(5.96851304800 + 1169.58825140000 * Tau)
167     + 0.00000056600 * cos(1.41198438840 + 533.62311836000 * Tau)
168     + 0.00000054935 * cos(5.42806383720 + 10.29494073800 * Tau)
169     + 0.00000052309 * cos(5.72661448390 + 117.31986822000 * Tau)
170     + 0.00000052016 * cos(0.22981299129 + 1368.66025280000 * Tau)
171     + 0.00000050418 * cos(6.08075147810 + 525.75881183000 * Tau)
172     + 0.00000047418 * cos(3.62611843240 + 1478.86657410000 * Tau)
173     + 0.00000046678 * cos(0.51144073175 + 1265.56747860000 * Tau)
174     + 0.00000039888 * cos(4.16158013600 + 1692.16566950000 * Tau)
175     + 0.00000033558 * cos(0.09913904872 + 302.16477566000 * Tau)
176     + 0.00000032827 * cos(5.03596689460 + 220.41264244000 * Tau)
177     + 0.00000032449 * cos(5.37492530700 + 508.35032409000 * Tau)
178     + 0.00000029483 * cos(5.42208897100 + 1272.68102560000 * Tau)
179     + 0.00000029379 * cos(3.35927241530 + 4.66586644600 * Tau)
180     + 0.00000029307 * cos(0.75907909735 + 88.86568021700 * Tau)
181     + 0.00000025195 * cos(1.60723063390 + 831.85574075000 * Tau);
182    
183     L2 := 0.00047233598 * cos(4.32148323550 + 7.11354700080 * Tau)
184     + 0.00038965550 * cos(0.00000000000 + 0.00000000000 * Tau)
185     + 0.00030629053 * cos(2.93021440220 + 529.69096509000 * Tau)
186     + 0.00003189317 * cos(1.05504615600 + 522.57741809000 * Tau)
187     + 0.00002729292 * cos(4.84545481350 + 536.80451210000 * Tau)
188     + 0.00002723358 * cos(3.41411526640 + 1059.38193020000 * Tau)
189     + 0.00001721069 * cos(4.18734385160 + 14.22709400200 * Tau)
190     + 0.00000383258 * cos(5.76790714390 + 419.48464388000 * Tau)
191     + 0.00000377524 * cos(0.76048964872 + 515.46387109000 * Tau)
192     + 0.00000367498 * cos(6.05509120410 + 103.09277422000 * Tau)
193     + 0.00000337386 * cos(3.78644384240 + 3.18139373770 * Tau)
194     + 0.00000308200 * cos(0.69356654052 + 206.18554844000 * Tau)
195     + 0.00000218408 * cos(3.81389191350 + 1589.07289530000 * Tau)
196     + 0.00000198883 * cos(5.33996443440 + 1066.49547720000 * Tau)
197     + 0.00000197445 * cos(2.48356402050 + 3.93215326310 * Tau)
198     + 0.00000155862 * cos(1.40642426470 + 1052.26838320000 * Tau)
199     + 0.00000146230 * cos(3.81373196840 + 639.89728631000 * Tau)
200     + 0.00000141932 * cos(1.63435169020 + 426.59819088000 * Tau)
201     + 0.00000129570 * cos(5.83738872530 + 412.37109687000 * Tau)
202     + 0.00000117327 * cos(1.41435462590 + 625.67019231000 * Tau)
203     + 0.00000096733 * cos(4.03383427890 + 110.20632122000 * Tau)
204     + 0.00000090823 * cos(1.10630629040 + 95.97922721800 * Tau)
205     + 0.00000087292 * cos(2.52235174820 + 632.78373931000 * Tau)
206     + 0.00000078769 * cos(4.63726131330 + 543.91805910000 * Tau)
207     + 0.00000072392 * cos(2.21716670030 + 735.87651353000 * Tau)
208     + 0.00000058475 * cos(0.83216317444 + 199.07200144000 * Tau)
209     + 0.00000056910 * cos(3.12292059850 + 213.29909544000 * Tau)
210     + 0.00000048622 * cos(1.67283791620 + 309.27832266000 * Tau)
211     + 0.00000040150 * cos(4.02485444740 + 21.34064100200 * Tau)
212     + 0.00000039784 * cos(0.62416945827 + 323.50541666000 * Tau)
213     + 0.00000035718 * cos(2.32581247000 + 728.76296653000 * Tau)
214     + 0.00000029255 * cos(3.60838327800 + 10.29494073800 * Tau)
215     + 0.00000027814 * cos(3.23992013740 + 838.96928775000 * Tau)
216     + 0.00000025993 * cos(4.50118298290 + 742.99006053000 * Tau)
217     + 0.00000025620 * cos(2.51240623860 + 1162.47470440000 * Tau)
218     + 0.00000025194 * cos(1.21868110690 + 1045.15483620000 * Tau)
219     + 0.00000023591 * cos(3.00532139310 + 956.28915597000 * Tau)
220     + 0.00000019458 * cos(4.29028644670 + 532.87235883000 * Tau)
221     + 0.00000017660 * cos(0.80953941560 + 508.35032409000 * Tau)
222     + 0.00000017058 * cos(4.20001977720 + 2118.76386040000 * Tau)
223     + 0.00000017040 * cos(1.83402146640 + 526.50957136000 * Tau)
224     + 0.00000015355 * cos(5.81037986940 + 1596.18644230000 * Tau)
225     + 0.00000015292 * cos(0.68174165476 + 942.06206197000 * Tau)
226     + 0.00000014661 * cos(3.99989622590 + 117.31986822000 * Tau)
227     + 0.00000013920 * cos(5.95169568480 + 316.39186966000 * Tau)
228     + 0.00000013639 * cos(1.80336677960 + 302.16477566000 * Tau)
229     + 0.00000013230 * cos(2.51856643600 + 88.86568021700 * Tau)
230     + 0.00000012756 * cos(4.36856232410 + 1169.58825140000 * Tau)
231     + 0.00000010986 * cos(4.43586634640 + 525.75881183000 * Tau)
232     + 0.00000009681 * cos(1.71563161050 + 1581.95934830000 * Tau)
233     + 0.00000009437 * cos(2.17684563460 + 1155.36115740000 * Tau)
234     + 0.00000008812 * cos(3.29452783340 + 220.41264244000 * Tau)
235     + 0.00000008690 * cos(3.31924493610 + 831.85574075000 * Tau)
236     + 0.00000007823 * cos(5.75672228350 + 846.08283475000 * Tau)
237     + 0.00000007549 * cos(2.70955516780 + 533.62311836000 * Tau)
238     + 0.00000006685 * cos(2.17560093280 + 1265.56747860000 * Tau)
239     + 0.00000006285 * cos(0.49939863541 + 949.17560897000 * Tau);
240    
241     L3 := 0.00006501665 * cos(2.59862880480 + 7.11354700080 * Tau)
242     + 0.00001356524 * cos(1.34635886410 + 529.69096509000 * Tau)
243     + 0.00000470716 * cos(2.47503977880 + 14.22709400200 * Tau)
244     + 0.00000416960 * cos(3.24451243210 + 536.80451210000 * Tau)
245     + 0.00000352851 * cos(2.97360159000 + 522.57741809000 * Tau)
246     + 0.00000154880 * cos(2.07565585820 + 1059.38193020000 * Tau)
247     + 0.00000086771 * cos(2.51431584320 + 515.46387109000 * Tau)
248     + 0.00000044378 * cos(0.00000000000 + 0.00000000000 * Tau)
249     + 0.00000033538 * cos(3.82633794500 + 1066.49547720000 * Tau)
250     + 0.00000028457 * cos(2.44754756060 + 206.18554844000 * Tau)
251     + 0.00000023737 * cos(1.27667172310 + 412.37109687000 * Tau)
252     + 0.00000022644 * cos(2.98231326770 + 543.91805910000 * Tau)
253     + 0.00000019798 * cos(2.10099934010 + 639.89728631000 * Tau)
254     + 0.00000019740 * cos(1.40255938970 + 419.48464388000 * Tau)
255     + 0.00000018768 * cos(1.59368403500 + 103.09277422000 * Tau)
256     + 0.00000017033 * cos(2.30214681200 + 21.34064100200 * Tau)
257     + 0.00000016774 * cos(2.59821460670 + 1589.07289530000 * Tau)
258     + 0.00000016214 * cos(3.14521117300 + 625.67019231000 * Tau)
259     + 0.00000016055 * cos(3.36030126300 + 1052.26838320000 * Tau)
260     + 0.00000013392 * cos(2.75973892200 + 95.97922721800 * Tau)
261     + 0.00000013234 * cos(2.53862244340 + 199.07200144000 * Tau)
262     + 0.00000012611 * cos(6.26578110400 + 426.59819088000 * Tau)
263     + 0.00000008701 * cos(1.76334960740 + 10.29494073800 * Tau)
264     + 0.00000008637 * cos(2.26563256290 + 110.20632122000 * Tau)
265     + 0.00000006725 * cos(3.42566433320 + 309.27832266000 * Tau)
266     + 0.00000006527 * cos(4.03869562910 + 728.76296653000 * Tau)
267     + 0.00000005675 * cos(2.52096417680 + 508.35032409000 * Tau)
268     + 0.00000005399 * cos(2.91184687110 + 1045.15483620000 * Tau)
269     + 0.00000005368 * cos(5.25196153540 + 323.50541666000 * Tau)
270     + 0.00000003996 * cos(4.30290261180 + 88.86568021700 * Tau)
271     + 0.00000003857 * cos(3.52381361550 + 302.16477566000 * Tau)
272     + 0.00000003774 * cos(4.09125315150 + 735.87651353000 * Tau)
273     + 0.00000003269 * cos(1.43175991270 + 956.28915597000 * Tau)
274     + 0.00000002783 * cos(4.35817507670 + 1596.18644230000 * Tau)
275     + 0.00000002661 * cos(1.25276590760 + 213.29909544000 * Tau)
276     + 0.00000002656 * cos(5.01505839850 + 838.96928775000 * Tau)
277     + 0.00000002553 * cos(2.23785673280 + 117.31986822000 * Tau)
278     + 0.00000002371 * cos(2.89662409240 + 742.99006053000 * Tau)
279     + 0.00000002279 * cos(2.35581871230 + 942.06206197000 * Tau);
280    
281     L4 := 0.00000669483 * cos(0.85282421090 + 7.11354700080 * Tau)
282     + 0.00000114019 * cos(3.14159265360 + 0.00000000000 * Tau)
283     + 0.00000099961 * cos(0.74258947751 + 14.22709400200 * Tau)
284     + 0.00000050024 * cos(1.65346208250 + 536.80451210000 * Tau)
285     + 0.00000043585 * cos(5.82026386620 + 529.69096509000 * Tau)
286     + 0.00000031813 * cos(4.85829986650 + 522.57741809000 * Tau)
287     + 0.00000014742 * cos(4.29061635780 + 515.46387109000 * Tau)
288     + 0.00000008899 * cos(0.71478520741 + 1059.38193020000 * Tau)
289     + 0.00000004957 * cos(1.29502259430 + 543.91805910000 * Tau)
290     + 0.00000004484 * cos(2.31715516630 + 1066.49547720000 * Tau)
291     + 0.00000004251 * cos(0.48326797501 + 21.34064100200 * Tau)
292     + 0.00000003100 * cos(3.00245542680 + 412.37109687000 * Tau)
293     + 0.00000002055 * cos(0.39858940218 + 639.89728631000 * Tau)
294     + 0.00000001902 * cos(4.25925620270 + 199.07200144000 * Tau)
295     + 0.00000001762 * cos(4.90536207310 + 625.67019231000 * Tau)
296     + 0.00000001695 * cos(4.26147580800 + 206.18554844000 * Tau)
297     + 0.00000001375 * cos(5.25546955670 + 1052.26838320000 * Tau)
298     + 0.00000001203 * cos(4.71614633840 + 95.97922721800 * Tau)
299     + 0.00000001086 * cos(1.28604571170 + 1589.07289530000 * Tau);
300    
301     L5 := 0.00000049577 * cos(5.25658966180 + 7.11354700080 * Tau)
302     + 0.00000015761 * cos(5.25126837480 + 14.22709400200 * Tau)
303     + 0.00000004343 * cos(0.01461869263 + 536.80451210000 * Tau)
304     + 0.00000001526 * cos(1.09739911440 + 522.57741809000 * Tau)
305     + 0.00000000845 * cos(3.14159265360 + 0.00000000000 * Tau);
306     Result := (L0 + L1*Tau + L2*Tau2 + L3*Tau3 + L4*Tau4 + L5*Tau5);
307     end;
308    
309     {---------------------------------------------------------------------------}
310    
311     function GetLatitude(Tau, Tau2, Tau3, Tau4, Tau5 : Double) : Double;
312     var
313     B0, B1,
314     B2, B3,
315     B4, B5 : Double;
316     begin
317     B0 := 0.02268615703 * cos(3.55852606720 + 529.69096509000 * Tau)
318     + 0.00110090358 * cos(0.00000000000 + 0.00000000000 * Tau)
319     + 0.00109971634 * cos(3.90809347390 + 1059.38193020000 * Tau)
320     + 0.00008101427 * cos(3.60509573370 + 522.57741809000 * Tau)
321     + 0.00006437782 * cos(0.30627121409 + 536.80451210000 * Tau)
322     + 0.00006043996 * cos(4.25883108790 + 1589.07289530000 * Tau)
323     + 0.00001106880 * cos(2.98534421930 + 1162.47470440000 * Tau)
324     + 0.00000944328 * cos(1.67522288400 + 426.59819088000 * Tau)
325     + 0.00000941651 * cos(2.93619072400 + 1052.26838320000 * Tau)
326     + 0.00000894088 * cos(1.75447429920 + 7.11354700080 * Tau)
327     + 0.00000835861 * cos(5.17881973230 + 103.09277422000 * Tau)
328     + 0.00000767280 * cos(2.15473594060 + 632.78373931000 * Tau)
329     + 0.00000684220 * cos(3.67808770100 + 213.29909544000 * Tau)
330     + 0.00000629223 * cos(0.64343282328 + 1066.49547720000 * Tau)
331     + 0.00000558524 * cos(0.01354830508 + 846.08283475000 * Tau)
332     + 0.00000531670 * cos(2.70305954350 + 110.20632122000 * Tau)
333     + 0.00000464449 * cos(1.17337249190 + 949.17560897000 * Tau)
334     + 0.00000431072 * cos(2.60825000490 + 419.48464388000 * Tau)
335     + 0.00000351433 * cos(4.61062990710 + 2118.76386040000 * Tau)
336     + 0.00000132160 * cos(4.77816990670 + 742.99006053000 * Tau)
337     + 0.00000123148 * cos(3.34968181380 + 1692.16566950000 * Tau)
338     + 0.00000116379 * cos(1.38688232030 + 323.50541666000 * Tau)
339     + 0.00000115038 * cos(5.04892295440 + 316.39186966000 * Tau)
340     + 0.00000103762 * cos(3.70103838110 + 515.46387109000 * Tau)
341     + 0.00000103402 * cos(2.31878999570 + 1478.86657410000 * Tau)
342     + 0.00000102420 * cos(3.15293785440 + 1581.95934830000 * Tau);
343    
344     B1 := 0.00177351787 * cos(5.70166488490 + 529.69096509000 * Tau)
345     + 0.00003230171 * cos(5.77941619340 + 1059.38193020000 * Tau)
346     + 0.00003081364 * cos(5.47464296530 + 522.57741809000 * Tau)
347     + 0.00002211914 * cos(4.73477480210 + 536.80451210000 * Tau)
348     + 0.00001694232 * cos(3.14159265360 + 0.00000000000 * Tau)
349     + 0.00000346445 * cos(4.74595174110 + 1052.26838320000 * Tau)
350     + 0.00000234264 * cos(5.18856099930 + 1066.49547720000 * Tau)
351     + 0.00000196154 * cos(6.18554286640 + 7.11354700080 * Tau)
352     + 0.00000150468 * cos(3.92721226090 + 1589.07289530000 * Tau)
353     + 0.00000114128 * cos(3.43897271830 + 632.78373931000 * Tau)
354     + 0.00000096667 * cos(2.91426304090 + 949.17560897000 * Tau)
355     + 0.00000081671 * cos(5.07666097500 + 1162.47470440000 * Tau)
356     + 0.00000076599 * cos(2.50522188660 + 103.09277422000 * Tau)
357     + 0.00000076572 * cos(0.61288981445 + 419.48464388000 * Tau)
358     + 0.00000073875 * cos(5.49958292150 + 515.46387109000 * Tau)
359     + 0.00000060544 * cos(5.44740084360 + 213.29909544000 * Tau)
360     + 0.00000049915 * cos(3.94799616570 + 735.87651353000 * Tau)
361     + 0.00000046032 * cos(0.53850360901 + 110.20632122000 * Tau)
362     + 0.00000045123 * cos(1.89516645240 + 846.08283475000 * Tau)
363     + 0.00000036561 * cos(4.69828392840 + 543.91805910000 * Tau)
364     + 0.00000036019 * cos(6.10952578760 + 316.39186966000 * Tau)
365     + 0.00000031975 * cos(4.92452714630 + 1581.95934830000 * Tau);
366    
367     B2 := 0.00008094051 * cos(1.46322843660 + 529.69096509000 * Tau)
368     + 0.00000813244 * cos(3.14159265360 + 0.00000000000 * Tau)
369     + 0.00000742415 * cos(0.95691639003 + 522.57741809000 * Tau)
370     + 0.00000398951 * cos(2.89888666450 + 536.80451210000 * Tau)
371     + 0.00000342226 * cos(1.44683789730 + 1059.38193020000 * Tau)
372     + 0.00000073948 * cos(0.40724675866 + 1052.26838320000 * Tau)
373     + 0.00000046151 * cos(3.48036895770 + 1066.49547720000 * Tau)
374     + 0.00000029717 * cos(1.92504171330 + 1589.07289530000 * Tau)
375     + 0.00000029314 * cos(0.99088831805 + 515.46387109000 * Tau)
376     + 0.00000022753 * cos(4.27124052440 + 7.11354700080 * Tau)
377     + 0.00000013916 * cos(2.92242387340 + 543.91805910000 * Tau)
378     + 0.00000012067 * cos(5.22168932480 + 632.78373931000 * Tau)
379     + 0.00000010703 * cos(4.88024222480 + 949.17560897000 * Tau)
380     + 0.00000006078 * cos(6.21089108430 + 1045.15483620000 * Tau);
381    
382     B3 := 0.00000251624 * cos(3.38087923080 + 529.69096509000 * Tau)
383     + 0.00000121738 * cos(2.73311837200 + 522.57741809000 * Tau)
384     + 0.00000048694 * cos(1.03689996680 + 536.80451210000 * Tau)
385     + 0.00000010988 * cos(2.31463561350 + 1052.26838320000 * Tau)
386     + 0.00000008067 * cos(2.76729757620 + 515.46387109000 * Tau)
387     + 0.00000007287 * cos(4.25268318970 + 1059.38193020000 * Tau)
388     + 0.00000006205 * cos(1.78115827370 + 1066.49547720000 * Tau)
389     + 0.00000003627 * cos(1.13028917220 + 543.91805910000 * Tau)
390     + 0.00000002798 * cos(3.14159265360 + 0.00000000000 * Tau);
391    
392     B4 := 0.00000015050 * cos(4.52956999640 + 522.57741809000 * Tau)
393     + 0.00000005370 * cos(4.47427159140 + 529.69096509000 * Tau)
394     + 0.00000004456 * cos(5.43908581050 + 536.80451210000 * Tau)
395     + 0.00000003422 * cos(0.00000000000 + 0.00000000000 * Tau)
396     + 0.00000001833 * cos(4.51807036230 + 515.46387109000 * Tau)
397     + 0.00000001322 * cos(4.20117611580 + 1052.26838320000 * Tau);
398    
399     B5 := 0.00000001445 * cos(0.09198554072 + 522.57741809000 * Tau);
400     Result := (B0 + B1*Tau + B2*Tau2 + B3*Tau3 + B4*Tau4 + B5*Tau5);
401     end;
402    
403     {---------------------------------------------------------------------------}
404    
405     function GetRadiusVector(Tau, Tau2, Tau3, Tau4, Tau5 : Double) : Double;
406     var
407     R0, R1,
408     R2, R3,
409     R4, R5 : Double;
410     begin
411     R0 := 5.20887429470 * cos(0.00000000000 + 0.00000000000 * Tau)
412     + 0.25209327020 * cos(3.49108640020 + 529.69096509000 * Tau)
413     + 0.00610599902 * cos(3.84115365600 + 1059.38193020000 * Tau)
414     + 0.00282029465 * cos(2.57419879930 + 632.78373931000 * Tau)
415     + 0.00187647391 * cos(2.07590380080 + 522.57741809000 * Tau)
416     + 0.00086792941 * cos(0.71001090609 + 419.48464388000 * Tau)
417     + 0.00072062869 * cos(0.21465694745 + 536.80451210000 * Tau)
418     + 0.00065517227 * cos(5.97995850840 + 316.39186966000 * Tau)
419     + 0.00030135275 * cos(2.16132058450 + 949.17560897000 * Tau)
420     + 0.00029134620 * cos(1.67759243710 + 103.09277422000 * Tau)
421     + 0.00023947340 * cos(0.27457854894 + 7.11354700080 * Tau)
422     + 0.00023453209 * cos(3.54023147300 + 735.87651353000 * Tau)
423     + 0.00022283710 * cos(4.19362773550 + 1589.07289530000 * Tau)
424     + 0.00013032600 * cos(2.96043055740 + 1162.47470440000 * Tau)
425     + 0.00012749004 * cos(2.71550102860 + 1052.26838320000 * Tau)
426     + 0.00009703346 * cos(1.90669572400 + 206.18554844000 * Tau)
427     + 0.00009161431 * cos(4.41352618940 + 213.29909544000 * Tau)
428     + 0.00007894539 * cos(2.47907551400 + 426.59819088000 * Tau)
429     + 0.00007057978 * cos(2.18184753110 + 1265.56747860000 * Tau)
430     + 0.00006137755 * cos(6.26417542510 + 846.08283475000 * Tau)
431     + 0.00005477093 * cos(5.65729325170 + 639.89728631000 * Tau)
432     + 0.00004170012 * cos(2.01605033910 + 515.46387109000 * Tau)
433     + 0.00004136890 * cos(2.72219979680 + 625.67019231000 * Tau)
434     + 0.00003502519 * cos(0.56531297394 + 1066.49547720000 * Tau)
435     + 0.00002616955 * cos(2.00993967130 + 1581.95934830000 * Tau)
436     + 0.00002499966 * cos(4.55182055940 + 838.96928775000 * Tau)
437     + 0.00002127644 * cos(6.12751461750 + 742.99006053000 * Tau)
438     + 0.00001911876 * cos(0.85621927419 + 412.37109687000 * Tau)
439     + 0.00001610549 * cos(3.08867789270 + 1368.66025280000 * Tau)
440     + 0.00001479484 * cos(2.68026191370 + 1478.86657410000 * Tau)
441     + 0.00001230708 * cos(1.89042979700 + 323.50541666000 * Tau)
442     + 0.00001216810 * cos(1.80171561020 + 110.20632122000 * Tau)
443     + 0.00001014959 * cos(1.38673237670 + 454.90936653000 * Tau)
444     + 0.00000998579 * cos(2.87208940110 + 309.27832266000 * Tau)
445     + 0.00000961072 * cos(4.54876989810 + 2118.76386040000 * Tau)
446     + 0.00000885708 * cos(4.14785948470 + 533.62311836000 * Tau)
447     + 0.00000821465 * cos(1.59342534400 + 1898.35121790000 * Tau)
448     + 0.00000812036 * cos(5.94091899140 + 909.81873305000 * Tau)
449     + 0.00000776700 * cos(3.67696954690 + 728.76296653000 * Tau)
450     + 0.00000727162 * cos(3.98824686400 + 1155.36115740000 * Tau)
451     + 0.00000655289 * cos(2.79065604220 + 1685.05212250000 * Tau)
452     + 0.00000653981 * cos(3.38150775270 + 1692.16566950000 * Tau)
453     + 0.00000620798 * cos(4.82284338960 + 956.28915597000 * Tau)
454     + 0.00000614784 * cos(2.27624915600 + 942.06206197000 * Tau)
455     + 0.00000562120 * cos(0.08095987241 + 543.91805910000 * Tau)
456     + 0.00000542221 * cos(0.28360266386 + 525.75881183000 * Tau);
457    
458     R1 := 0.01271801596 * cos(2.64937511120 + 529.69096509000 * Tau)
459     + 0.00061661771 * cos(3.00076251020 + 1059.38193020000 * Tau)
460     + 0.00053443592 * cos(3.89717644230 + 522.57741809000 * Tau)
461     + 0.00041390257 * cos(0.00000000000 + 0.00000000000 * Tau)
462     + 0.00031185167 * cos(4.88276663530 + 536.80451210000 * Tau)
463     + 0.00011847190 * cos(2.41329588180 + 419.48464388000 * Tau)
464     + 0.00009166360 * cos(4.75979408590 + 7.11354700080 * Tau)
465     + 0.00003403605 * cos(3.34688538000 + 1589.07289530000 * Tau)
466     + 0.00003203446 * cos(5.21083285480 + 735.87651353000 * Tau)
467     + 0.00003175763 * cos(2.79297987070 + 103.09277422000 * Tau)
468     + 0.00002806064 * cos(3.74223693580 + 515.46387109000 * Tau)
469     + 0.00002676575 * cos(4.33052878700 + 1052.26838320000 * Tau)
470     + 0.00002600003 * cos(3.63435101620 + 206.18554844000 * Tau)
471     + 0.00002412207 * cos(1.46947308300 + 426.59819088000 * Tau)
472     + 0.00002100507 * cos(3.92762682310 + 639.89728631000 * Tau)
473     + 0.00001646182 * cos(5.30953510950 + 1066.49547720000 * Tau)
474     + 0.00001641257 * cos(4.41628669820 + 625.67019231000 * Tau)
475     + 0.00001049866 * cos(3.16113622960 + 213.29909544000 * Tau)
476     + 0.00001024802 * cos(2.55432643020 + 412.37109687000 * Tau)
477     + 0.00000806404 * cos(2.67750801380 + 632.78373931000 * Tau)
478     + 0.00000740996 * cos(2.17094630560 + 1162.47470440000 * Tau)
479     + 0.00000676928 * cos(6.24953479790 + 838.96928775000 * Tau)
480     + 0.00000567076 * cos(4.57655414710 + 742.99006053000 * Tau)
481     + 0.00000484689 * cos(2.46882793190 + 949.17560897000 * Tau)
482     + 0.00000468895 * cos(4.70973463480 + 543.91805910000 * Tau)
483     + 0.00000444683 * cos(0.40281181402 + 323.50541666000 * Tau)
484     + 0.00000415894 * cos(5.36836018210 + 728.76296653000 * Tau)
485     + 0.00000401738 * cos(4.60528841540 + 309.27832266000 * Tau)
486     + 0.00000347378 * cos(4.68148808720 + 14.22709400200 * Tau)
487     + 0.00000337555 * cos(3.16781951120 + 956.28915597000 * Tau)
488     + 0.00000260753 * cos(5.34290306100 + 846.08283475000 * Tau)
489     + 0.00000246603 * cos(3.92313823540 + 942.06206197000 * Tau)
490     + 0.00000220084 * cos(4.84210964960 + 1368.66025280000 * Tau)
491     + 0.00000203217 * cos(5.59995425430 + 1155.36115740000 * Tau)
492     + 0.00000200190 * cos(4.43888814440 + 1045.15483620000 * Tau)
493     + 0.00000197134 * cos(3.70551461390 + 2118.76386040000 * Tau)
494     + 0.00000196005 * cos(3.75877587140 + 199.07200144000 * Tau)
495     + 0.00000183504 * cos(4.26526769700 + 95.97922721800 * Tau)
496     + 0.00000180134 * cos(4.40165491160 + 532.87235883000 * Tau)
497     + 0.00000170225 * cos(4.84647488870 + 526.50957136000 * Tau)
498     + 0.00000146335 * cos(6.12958365530 + 533.62311836000 * Tau)
499     + 0.00000133483 * cos(1.32245735860 + 110.20632122000 * Tau)
500     + 0.00000132076 * cos(4.51187950810 + 525.75881183000 * Tau);
501    
502     R2 := 0.00079644833 * cos(1.35865896600 + 529.69096509000 * Tau)
503     + 0.00008251618 * cos(5.77773935440 + 522.57741809000 * Tau)
504     + 0.00007029864 * cos(3.27476965830 + 536.80451210000 * Tau)
505     + 0.00005314006 * cos(1.83835109710 + 1059.38193020000 * Tau)
506     + 0.00001860833 * cos(2.97682139370 + 7.11354700080 * Tau)
507     + 0.00000964466 * cos(5.48031822020 + 515.46387109000 * Tau)
508     + 0.00000836267 * cos(4.19889881720 + 419.48464388000 * Tau)
509     + 0.00000497920 * cos(3.14159265360 + 0.00000000000 * Tau)
510     + 0.00000426570 * cos(2.22753101800 + 639.89728631000 * Tau)
511     + 0.00000406453 * cos(3.78250730350 + 1066.49547720000 * Tau)
512     + 0.00000377316 * cos(2.24248352870 + 1589.07289530000 * Tau)
513     + 0.00000362943 * cos(5.36761847270 + 206.18554844000 * Tau)
514     + 0.00000342048 * cos(6.09922969320 + 1052.26838320000 * Tau)
515     + 0.00000339043 * cos(6.12690864040 + 625.67019231000 * Tau)
516     + 0.00000332578 * cos(0.00328961161 + 426.59819088000 * Tau)
517     + 0.00000279920 * cos(4.26162555830 + 412.37109687000 * Tau)
518     + 0.00000257290 * cos(0.96295364983 + 632.78373931000 * Tau)
519     + 0.00000229777 * cos(0.70530766213 + 735.87651353000 * Tau)
520     + 0.00000200783 * cos(3.06850623370 + 543.91805910000 * Tau)
521     + 0.00000199807 * cos(4.42884165320 + 103.09277422000 * Tau)
522     + 0.00000138606 * cos(2.93235671610 + 14.22709400200 * Tau)
523     + 0.00000113535 * cos(0.78713911289 + 728.76296653000 * Tau)
524     + 0.00000094565 * cos(1.70498041070 + 838.96928775000 * Tau)
525     + 0.00000086025 * cos(5.14434751990 + 323.50541666000 * Tau)
526     + 0.00000083469 * cos(0.05834873484 + 309.27832266000 * Tau)
527     + 0.00000080328 * cos(2.98122361800 + 742.99006053000 * Tau)
528     + 0.00000075198 * cos(1.60495195910 + 956.28915597000 * Tau)
529     + 0.00000070451 * cos(1.50988357480 + 213.29909544000 * Tau)
530     + 0.00000066572 * cos(5.47307178080 + 199.07200144000 * Tau)
531     + 0.00000061649 * cos(6.10137889850 + 1045.15483620000 * Tau)
532     + 0.00000056203 * cos(0.95534810533 + 1162.47470440000 * Tau)
533     + 0.00000051904 * cos(5.58435625610 + 942.06206197000 * Tau)
534     + 0.00000050057 * cos(2.72063162320 + 532.87235883000 * Tau)
535     + 0.00000044548 * cos(5.52445621410 + 508.35032409000 * Tau)
536     + 0.00000044282 * cos(0.27118152557 + 526.50957136000 * Tau)
537     + 0.00000039833 * cos(5.94566506230 + 95.97922721800 * Tau);
538    
539     R3 := 0.00003519257 * cos(6.05800633850 + 529.69096509000 * Tau)
540     + 0.00001073239 * cos(1.67321345760 + 536.80451210000 * Tau)
541     + 0.00000915666 * cos(1.41329676120 + 522.57741809000 * Tau)
542     + 0.00000341593 * cos(0.52296542656 + 1059.38193020000 * Tau)
543     + 0.00000254893 * cos(1.19625473530 + 7.11354700080 * Tau)
544     + 0.00000221512 * cos(0.95225226237 + 515.46387109000 * Tau)
545     + 0.00000089729 * cos(3.14159265360 + 0.00000000000 * Tau)
546     + 0.00000069078 * cos(2.26885282310 + 1066.49547720000 * Tau)
547     + 0.00000057827 * cos(1.41389745340 + 543.91805910000 * Tau)
548     + 0.00000057653 * cos(0.52580117593 + 639.89728631000 * Tau)
549     + 0.00000051079 * cos(5.98016364680 + 412.37109687000 * Tau)
550     + 0.00000046935 * cos(1.57864237960 + 625.67019231000 * Tau)
551     + 0.00000042824 * cos(6.11689609100 + 419.48464388000 * Tau)
552     + 0.00000037477 * cos(1.18262762330 + 14.22709400200 * Tau)
553     + 0.00000033816 * cos(1.66671706950 + 1052.26838320000 * Tau)
554     + 0.00000033531 * cos(0.84784977903 + 206.18554844000 * Tau)
555     + 0.00000031195 * cos(1.04290245900 + 1589.07289530000 * Tau)
556     + 0.00000030023 * cos(4.63236245030 + 426.59819088000 * Tau)
557     + 0.00000020804 * cos(2.50071243810 + 728.76296653000 * Tau)
558     + 0.00000015023 * cos(0.89136998434 + 199.07200144000 * Tau)
559     + 0.00000014466 * cos(0.96040197071 + 508.35032409000 * Tau)
560     + 0.00000012969 * cos(1.50233788550 + 1045.15483620000 * Tau)
561     + 0.00000012319 * cos(2.60952614500 + 735.87651353000 * Tau)
562     + 0.00000011654 * cos(3.55513510120 + 323.50541666000 * Tau)
563     + 0.00000011160 * cos(1.79041437550 + 309.27832266000 * Tau)
564     + 0.00000010554 * cos(6.27845112680 + 956.28915597000 * Tau)
565     + 0.00000009812 * cos(6.26016859520 + 103.09277422000 * Tau)
566     + 0.00000009301 * cos(3.45126812480 + 838.96928775000 * Tau);
567    
568     R4 := 0.00000128628 * cos(0.08419309557 + 536.80451210000 * Tau)
569     + 0.00000113458 * cos(4.24858855780 + 529.69096509000 * Tau)
570     + 0.00000082650 * cos(3.29754909410 + 522.57741809000 * Tau)
571     + 0.00000037883 * cos(2.73326611140 + 515.46387109000 * Tau)
572     + 0.00000026694 * cos(5.69142588560 + 7.11354700080 * Tau)
573     + 0.00000017650 * cos(5.40012536920 + 1059.38193020000 * Tau)
574     + 0.00000012612 * cos(6.01560416060 + 543.91805910000 * Tau)
575     + 0.00000009287 * cos(0.76813946494 + 1066.49547720000 * Tau)
576     + 0.00000008107 * cos(5.68228065710 + 14.22709400200 * Tau)
577     + 0.00000006978 * cos(1.42751292060 + 412.37109687000 * Tau)
578     + 0.00000006271 * cos(5.12286932530 + 639.89728631000 * Tau)
579     + 0.00000005377 * cos(3.33501947270 + 625.67019231000 * Tau)
580     + 0.00000002911 * cos(3.40334805050 + 1052.26838320000 * Tau)
581     + 0.00000002593 * cos(4.16090412980 + 728.76296653000 * Tau)
582     + 0.00000002562 * cos(2.89802035070 + 426.59819088000 * Tau);
583    
584     R5 := 0.00000011188 * cos(4.75249399940 + 536.80451210000 * Tau)
585     + 0.00000004255 * cos(5.91516229170 + 522.57741809000 * Tau)
586     + 0.00000002079 * cos(5.56781555860 + 515.46387109000 * Tau)
587     + 0.00000001908 * cos(4.29659647290 + 543.91805910000 * Tau)
588     + 0.00000001875 * cos(3.69357495840 + 7.11354700080 * Tau)
589     + 0.00000001612 * cos(4.13222808530 + 1059.38193020000 * Tau)
590     + 0.00000001590 * cos(5.49312796170 + 1066.49547720000 * Tau);
591    
592     Result := (R0 + R1*Tau + R2*Tau2 + R3*Tau3 + R4*Tau4 + R5*Tau5);
593     end;
594    
595     {---------------------------------------------------------------------------}
596    
597     function ComputeJupiter(JD : Double) : TStEclipticalCord;
598     var
599     Tau,
600     Tau2,
601     Tau3,
602     Tau4,
603     Tau5 : Double;
604     begin
605     Tau := (JD - 2451545.0) / 365250.0;
606     Tau2 := sqr(Tau);
607     Tau3 := Tau * Tau2;
608     Tau4 := sqr(Tau2);
609     Tau5 := Tau2 * Tau3;
610    
611     Result.L0 := GetLongitude(Tau, Tau2, Tau3, Tau4, Tau5);
612     Result.B0 := GetLatitude(Tau, Tau2, Tau3, Tau4, Tau5);
613     Result.R0 := GetRadiusVector(Tau, Tau2, Tau3, Tau4, Tau5);
614     end;
615    
616     end.
  ViewVC Help
Powered by ViewVC 1.1.20