hp41programs

Astronomical Ephemeris 4 Astronomical Ephemeris (IV) for the HP41
 

Overview
 

-The following program calculates the geocentric and heliocentric longitudes and latitudes of the Sun and the major planets, including Pluto.
-They are referred to the standard equinox J2000.
-The data are given by the JPL in reference [1] and the results are valid between  3000 BC and AD 3000
-Low accuracy is to be expected for Saturn & Uranus and good accuracy for the Sun, Mercury, Venus and Mars.
 

Data Registers
 

-Registers R00 thru R15 are used.

 R00 = T = time since 2000/01/01 0h TT is to be initialized before executing this program   -5 < T < 1         ( Unit = 1000 years = 365250 days )
 R01, R02, R03 contain the Sun's rectangular ecliptic coordinates. Execute "SUN"  first.
 R04 = geocentric longitude in decimal degrees
 R05 = geocentric latitude    ------------------
 R06 = distance to the Earth in Astronomical Units
 R07 = heliocentric longitude in decimal degrees
 R08 = heliocentric latitude   -------------------
 R09 = radius vector in Astronomical Units

 R10 = semi-major axis
 R11 = eccentricity
 R12 = inclination over the ecliptic J2000
 R13 = mean longitude
 R14 = longitude of the perihelion
 R15 = longitude of the ascending node
 

 >>>> The HP-41 must be set in DEG mode.
 

Program Listing
 

  LBL 01 = Mercury
  LBL 02 = Venus
  LBL 04 = Mars
  LBL 05 = Jupiter
  LBL 06 = Saturn
  LBL 07 = Uranus
  LBL 08 = Neptune
  LBL 09 = Pluto
 
 

  01  LBL "SUN"
  02  CLX
  03  STO 01
  04  STO 02
  05  STO 03
  06  SIGN
  07  STO 10
  08  1673
  09  RCL 00
  10  37
  11  *
  12  -
  13  STO 11
  14  1
  15  RCL 00
  16  134
  17  *
  18  +
  19  STO 12
  20  359993731
  21  RCL 00
  22  *
  23  80026
  24  -
  25  STO 13
  26  3180
  27  RCL 00
  28  *
  29  77070
  30  -
  31  STO 14
  32  174887
  33  RCL 00
  34  2412
  35  *
  36  -
  37  XEQ 14
  38  RCL 06
  39  XEQ 12
  40  STO 01
  41  RDN
  42  STO 02
  43  X<>Y
  44  STO 03
  45  RCL 06
  46  RCL 05
  47  RCL 04
  48  RTN
  49  LBL 01
  50  .3781
  51  STO 10
  52  20564
  53  RCL 00
  54  21
  55  *
  56  +
  57  STO 11
  58  7006
  59  RCL 00
  60  59
  61  *
  62  -
  63  STO 12
  64  RCL 00
  65  1494726749
  66  *
  67  250206
  68  +
  69  STO 13
  70  RCL 00
  71  1594
  72  *
  73  77458
  74  +
  75  STO 14
  76  48340
  77  RCL 00
  78  1221
  79  *
  80  -
  81  XEQ 14
  82  RTN
  83  LBL 02
  84  .7233
  85  STO 10        
  86  676
  87  RCL 00
  88  51
  89  *
  90  -
  91  STO 11
  92  3398
  93  RCL 00
  94  4
  95  *
  96  +
  97  STO 12
  98  RCL 00
  99  585178156
100  *
101  181179
102  +
103  STO 13
104  131768
105  RCL 00
106  568
107  *
108  +
109  STO 14
110  76673
111  RCL 00
112  2727
113  *
114  -
115  XEQ 14
116  RTN
117  LBL 04
118  1.5237
119  STO 10
120  9337
121  RCL 00
122  92
123  *
124  +
125  STO 11
126  1852
127  RCL 00
128  73
129  *
130  -
131  STO 12
132  RCL 00
133  191402993
134  *
135  4830
136  -
137  STO 13
138  RCL 00
139  4522
140  *
141  23917
142  -
143  STO 14
144  49713
145  RCL 00
146  2685
147  *
148  -
149  XEQ 14
150  RTN
151  LBL 05
152  5.202
153  STO 10
154  4854
155  RCL 00
156  180
157  *
158  +
159  STO 11
160  1299
161  RCL 00
162  32
163  *
164  -
165  STO 12        
166  30349037
167  RCL 00
168  12
169  *
170  -
171  RCL 00
172  *
173  34293
174  +
175  RCL 00
176  383.5
177  *
178  356
179  P-R
180  6
181  /
182  -
183  -
184  STO 13
185  14275
186  RCL 00
187  1820
188  *
189  +
190  STO 14
191  100293
192  RCL 00
193  1303
194  *
195  +
196  XEQ 14
197  RTN
198  LBL 06
199  9.541
200  STO 10
201  5551
202  RCL 00
203  320
204  *
205  -
206  STO 11
207  2494
208  RCL 00
209  45
210  *
211  +
212  STO 12
213  12221149
214  RCL 00
215  26
216  *
217  +
218  RCL 00
219  *
220  50059
221  +
222  RCL 00
223  383.5
224  *
225  873
226  P-R
227  6.5
228  /
229  -
230  +
231  STO 13
232  92861
233  RCL 00
234  5418
235  *
236  +
237  STO 14
238  113640
239  RCL 00
240  2502
241  *
242  -
243  XEQ 14
244  RTN
245  LBL 07
246  19.188
247  RCL 00        
248  500
249  /
250  -
251  STO 10
252  4686
253  RCL 00
254  16
255  *
256  -
257  STO 11
258  773
259  RCL 00
260  18
261  *
262  -
263  STO 12
264  4284951
265  RCL 00
266  58
267  *
268  +
269  RCL 00
270  *
271  45803
272  -
273  RCL 00
274  76.7
275  *
276  177
277  P-R
278  5.52
279  *
280  -
281  +
282  STO 13
283  172434
284  RCL 00
285  927
286  *
287  +
288  STO 14
289  73963
290  RCL 00
291  574
292  *
293  +
294  XEQ 14
295  RTN
296  LBL 08
297  30.07
298  STO 10
299  895
300  RCL 00
301  8
302  *
303  +
304  STO 11
305  1770
306  RCL 00
307  ST+ X
308  +
309  STO 12
310  2184652
311  RCL 00
312  41
313  *
314  -
315  RCL 00
316  *
317  55780
318  -
319  RCL 00
320  76.7
321  *
322  102
323  P-R
324  6.7
325  *
326  -
327  -
328  STO 13        
329  46682
330  RCL 00
331  101
332  *
333  +
334  STO 14
335  131786
336  RCL 00
337  61
338  *
339  -
340  XEQ 14
341  RTN
342  LBL 09
343  39.487
344  RCL 00
345  22
346  /
347  +
348  STO 10
349  24885
350  RCL 00
351  60
352  *
353  +
354  STO 11
355  17141
356  STO 12
357  1451804
358  RCL 00
359  1263
360  *
361  -
362  RCL 00
363  *
364  238963
365  +
366  STO 13
367  224097
368  RCL 00
369  97
370  *
371  -
372  STO 14
373  110302
374  RCL 00
375  81
376  *
377  -
378  GTO 14
379  LBL 12
380  X<>Y
381  RDN
382  P-R
383  R^
384  X<>Y
385  P-R
386  RTN
387  LBL 13
388  R-P
389  X<>Y
390  RDN
391  R-P
392  R^
393  X<>Y
394  RTN
395  LBL 14
396   E3
397  ST/ 12
398  ST/ 13
399  ST/ 14
400  /
401  STO 15
402   E5
403  ST/ 11
404  9
405  STO 09
406  RCL 13
407  RCL 14        
408  -
409  360
410  MOD
411  ENTER^
412  ENTER^
413  LBL 11
414  SIN
415  RCL 11
416  R-D
417  *
418  +
419  DSE 09
420  GTO 11
421  2
422  /
423  1
424  RCL 11
425  +
426  1
427  LASTX
428  -
429  /
430  SQRT
431  P-R
432  LASTX
433  /
434  R-P
435  RDN
436  ST+ X
437  1
438  R^
439  ST+ X
440  COS
441  RCL 11
442  *
443  -
444  RCL 10
445  *
446  STO 09
447  CLX
448  RCL 15
449  -
450  RCL 14
451  +
452  COS
453  RCL 12
454  LASTX
455  SIN
456  P-R
457  X<>Y
458  ASIN
459  STO 08
460  X<> Z
461  R-P
462  CLX
463  RCL 15
464  +
465  STO 07
466  RCL 09
467  XEQ 12
468  RCL 03
469  ST+ T
470  CLX
471  RCL 02
472  ST+ Z
473  CLX
474  RCL 01
475  +
476  XEQ 13
477  STO 06
478  X<> Z
479  STO 05
480  X<>Y
481  STO 04
482  END
 

 
    ( 899 bytes / SIZE 016 )
 
 

   STACK    INPUTS                 OUTPUTS
        Z         /      distance to the earth ( AU )
        Y         /       geocentric latitude ( deg )
        X         /      geocentric longitude ( deg )

 
Example:       T = -3   STO 00    i-e    -1001.1209   Gregorian =  -1001.1219  Julian      (  year -1001 =  1002 BC )
 

   •   XEQ "SUN"    >>>>    Long = -60°324              =  R04                     R01 = X =  0.487138
                               RDN      Lat  =  -0°356               =  R05         and       R02 = Y = -0.854871
                               RDN     Dist  =  0.983943 AU     =  R06                     R03 = Z =  -0.006105

   •   Mercury    R/S  or   XEQ 01   >>>>    Long = -78°210              =  R04                     R07 = -122°772
                                                       RDN      Lat  =  -0°516               =  R05         and       R08 = -0°658
                                                       RDN     Dist  =  1.265888 AU     =  R06                     R09 =  0.462176 AU

   •     Venus     R/S  or   XEQ 02   >>>>    Long = -32°829              =  R04                      R07 =  5°941
                                                       RDN      Lat  =  -1°920               =  R05         and       R08 = -3°323
                                                       RDN     Dist  =  1.4392 AU         =  R06                     R09 =  0.7266 AU

   •     Mars      R/S  or   XEQ 04   >>>>     Long = -29°972              =  R04                      R07 = -9°152
                                                       RDN      Lat  =  -1°398               =  R05         and       R08 = -1°905
                                                       RDN     Dist  =  2.157 AU           =  R06                     R09 =  1.400 AU

   •    Jupiter     R/S  or   XEQ 05   >>>>     Long =  62°420              =  R04                      R07 =  71°774
                                                       RDN      Lat  =  -0°736               =  R05         and       R08 = -0°581
                                                       RDN     Dist  =  4.493 AU          =  R06                      R09 =  5.092 AU

   •   Saturn      R/S  or   XEQ 06   >>>>    Long = 108°735              =  R04                      R07 =  109°922
                                                       RDN      Lat  =  -0°558              =  R05         and       R08 = -0°459
                                                       RDN     Dist  =  8.042 AU          =  R06                     R09 =  9.010 AU

   •  Uranus     R/S  or   XEQ 07     >>>>    Long =  52°952              =  R04                      R07 =  55°605
                                                       RDN      Lat  =  -0°260               =  R05         and       R08 = -0°237
                                                       RDN     Dist  =  19.118 AU         =  R06                     R09 =  19.528 AU

   •  Neptune   R/S  or   XEQ 08     >>>>    Long = -128°938            =  R04                      R07 =  229°330
                                                       RDN      Lat  =  1°718                 =  R05         and       R08 =  1°749
                                                       RDN     Dist  =  30.676 AU         =  R06                     R09 =  30.332 AU

   •   Pluto      R/S  or   XEQ 09     >>>>    Long = 172°634              =  R04                      R07 =  171°182
                                                       RDN      Lat  = 15°319               =  R05         and       R08 =  15°045
                                                       RDN     Dist  =  31.498 AU        =  R06                     R09 =  32.081 AU
 

Reference:

[1]    http://ssd.jpl.nasa.gov/txt/aprx_pos_planets.pdf