hp41programs

JSUN

Heliocentric Coordinates of the Giant Planets for the HP-41:  2025/04/01->2026/04/20


Overview
 

-This routine calculates the heliocentric coordinates of Jupiter, Saturn, Uranus & Neptune.
-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 2025/10/10 0h TT   ( unit = 192 days )    ( -1 <= R00 <= +1 )

                                        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" or "J2"  ( cf "Julian & Gregorian Calendars for the HP41" )
 

 01 LBL "JSUN"
 02 HR
 03 24
 04 /
 05 X<>Y
 06 XEQ "J1"
 07 +
 08 9414
 09 -
 10 192
 11 /
 12 ENTER
 13 ENTER
 14 STO 00
 15 2
 16 *
 17 CHS
 18 11703
 19 +
 20 *
 21 8527
 22 +
 23 STO 10          
 24 CLX
 25 3
 26 *
 27 236
 28 -		  29 *
 30 13299
 31 -
 32 STO 11
 33 CLX
 34 13
 35 *
 36 CHS
 37 3
 38 -
 39 *
 40 34
 41 +
 42 *
 43 58
 44 +
 45 *
 46 4957
 47 -
 48 *
 49 29886934
 50 +
 51 STO 12          
 52 CLX
 53 34
 54 *
 55 21895
 56 +		  57 *
 58 589557
 59 +
 60 STO 07
 61 CLX
 62 2
 63 *
 64 285
 65 +
 66 *
 67 2025
 68 -
 69 STO 08          
 70 CLX
 71 14
 72 *
 73 CHS
 74 8
 75 +
 76 *
 77 34
 78 +
 79 *
 80 375
 81 -
 82 *
 83 33613
 84 -		  85 *
 86 19504737
 87 +
 88 STO 09
 89 CLX
 90 3
 91 *
 92 392
 93 +
 94 *
 95 64219
 96 +
 97 *
 98 3591328
 99 +
100 STO 04          
101 CLX
102 4
103 *
104 134
105 +
106 *
107 1166
108 -
109 *
110 22574
111 -
112 STO 05		 113 CLX
114 12
115 *
116 CHS
117 2
118 -
119 *
120 141
121 +
122 *
123 147
124 -
125 *
126 59222
127 -
128 *
129 9544159
130 +
131 STO 06          
132 CLX
133 14
134 *
135 23
136 +
137 *
138 2167
139 -
140 *		 141 160679
142 +
143 *
144 1024247
145 +
146 STO 01          
147 CLX
148 47
149 *
150 CHS
151 62
152 -
153 *
154 3651
155 +
156 *
157 384
158 +
159 4
160 10^X
161 ST/ 01
162 ST/ 04
163 ST/ 05
164 ST/ 07
165 ST/ 08
166 ST/ 10
167 ST/ 11
168 /
 169 STO 02
170 CLX
171 PI
172 *
173 CHS
174 26
175 +
176 *
177 913
178 -
179 *
180 324
181 +
182 *
183 69975
184 +
185 *
186 5181657
187 +
188 6
189 10^X
190 ST/ 06
191 ST/ 09
192 ST/ 12
193 /
194 STO 03          
195 RCL 02
196 RCL 01
197 END

    ( 355 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 2025/07/07 at 16h41  TT

   2025.0707   ENTER^
       16.41       XEQ "JSUN"   >>>>  LJup =   94°4801         = R01
                                               RDN  BJup =   -0°1419          = R02
                                              RDN   RJup =  5.147475 AU  = R03

-And in registers R04 thru R12:

       R04 = LSat =   355°9880          R07 = LUra =   57°8811           R10 = LNep =    0°2778
       R05 = bSat =   -2°1969             R08 = bUra  =  -0°2165            R11 = bNep =  -1°3182
       R06 = RSat =  9.573195            R09 = RUra =  19.521153        R12 = RNep =  29.889379
 

Note:

-These polynomial approximations are fitted to JPL DE441 ephemerides.

 
 

References:

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