hp41programs

JSUN

2021 - Heliocentric Coordinates of the Giant Planets for the HP-41


Overview
 

-This routine calculates the heliocentric coordinates of Jupiter, Saturn, Uranus & Neptune from 2020/12/31 to 2022/01/19
-The ecliptic longitudes L & latitudes B are referred to the mean ecliptic & equinox of the date.
 

Program Listing
 

-All the angles are expressed in decimal degrees, all the distances R Sun-Planet in Astronomical Units
 
 

Data Registers:              R00 = time since 2020/12/31   ( unit = (1/384) day )

                                        R01 = LJup   R04 = LSat   R07 = LUra   R10 = LNep
                                        R02 = BJup   R05 = BSat   R08 = BUra   R11 = BNep
                                        R03 = RJup   R06 = RSat   R09 = RUra   R12 = RNep
Flags: /
Subroutine:  "J1"   ( cf "Julian & Gregorian Calendars for the HP41" )
 
 
 

 01 LBL "JSUN"
 02 HR
 03 24
 04 /
 05 X<>Y
 06 XEQ "J1"
 07 +
 08 7862
 09 -
 10 192
 11 /
 12 ENTER
 13 ENTER
 14 STO 00
 15 11
 16 *
 17 214
 18 +
 19 *
 20 59071
 21 +
 22 *
 23 3095005
 24 +
 25 STO 04
 26 CLX
 27 4
 28 *
 29 27
 30 +
 31 *
 32 2465
 33 -
 34 *
 35 6712
 36 -
 37 STO 05
 38 CLX
 39 14
 40 *
 41 3
 42 +
 43 *
 44 48
 45 +
 46 *
 47 2595
 48 -
 49 *
 50 35613
 51 -
 52 *
 53 9954212
 54 +
 55 STO 06
 56 CLX
 57 34
 58 *
 59 21335
 60 +
 61 *
 62 414831
 63 +
 64 STO 07
 65 CLX
 66 3
 67 *
 68 241
 69 +
 70 *
 71 4166
 72 -
 73 STO 08
 74 CLX
 75 7
 76 *
 77 13
 78 +
 79 *
 80 3
 81 +
 82 *
 83 437
 84 +
 85 *
 86 26033
 87 +
 88 *
 89 CHS
 90 19746697
 91 +
 92 STO 09
 93 CLX
 94 4
 95 *
 96 11702
 97 +
 98 *
 99 3513871
100 +
101 STO 10
102 RDN
103 ST+ X
104 277
105 -
106 *
107 11218
108 -
109 STO 11
110 CLX
111 7
112 *
113 8
114 -
115 *
116 26
117 -
118 *
119 127
120 -
121 *
122 3499
123 -
124 *
125 29923557
126 +
127 STO 12
128 CLX
129 19
130 *
131 CHS
132 130
133 -
134 *
135 1847
136 +
137 *
138 170166
139 +
140 *
141 3236758
142 +
143 STO 01
144 CLX
145 52
146 *
147 359
148 +
149 *
150 2830
151 -
152 *
153 8882
154 -
155  E4
156 ST/ 01
157 ST/ 04
158 ST/ 05
159 ST/ 07
160 ST/ 08
161 ST/ 10
162 ST/ 11
163 /
164 STO 02
165 CLX
166 45
167 *
168 CHS
169 821
170 +
171 *
172 6676
173 +
174 *
175 54325
176 -
177 *
178 5035918
179 +
180  E6
181 ST/ 06
182 ST/ 09
183 ST/ 12
184 /
185 STO 03
186 RCL 02
187 RCL 01
188 END

 
     ( 354 bytes / SIZE 013 )
 
 

           STACK             INPUTS         OUTPUTS
               Z                  /             RJup 
               Y       YYYY.MNDD             BJup
               X        HH.MNSS (TT)             LJup 


Example:   Calculate the heliocentric ecliptic coordinates of Jupiter Saturn Uranus & Neptune on 2021/12/25 at 16h41  TT

   2021.1225   ENTER^
       16.41       XEQ "JSUN"   >>>>  LJup =  338°6694        = R01
                                               RDN  BJup =   -1°1045          = R02
                                              RDN   RJup =  4.994083 AU  = R03

-And in registers R04 thru R12:

       R04 = LSat =   314.6769          R07 = LUra =   43.3491           R10 = LNep = 352.4095
       R05 = bSat =   -0.8842             R08 = bUra  =  -0.3953            R11 = bNep =  -1.1458
       R06 = RSat =  9.921168           R09 = RUra =  19.723613        R12 = RNep =  29.920386
 

Note:

-These polynomial approximations are fitted to JPL DE431 ephemerides.

 
 

References:

[1]  ftp://ssd.jpl.nasa.gov/pub/eph/planets/ascii/
[2]  Aldo Vitagliano SOLEX  http://www.solexorb.it/