Heliocentric Coordinates of the Giant Planets for the HP-41: 2024/04/01->2025/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 2024/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 9049 09 - 10 192 11 / 12 ENTER 13 ENTER 14 STO 00 15 2 16 * 17 CHS 18 11709 19 + 20 * 21 3586273 22 + 23 STO 10 24 CLX 25 3 26 * 27 247 |
28
- 29 * 30 12839 31 - 32 STO 11 33 CLX 34 17 35 * 36 11 37 + 38 * 39 5 40 - 41 * 42 42 43 - 44 * 45 5018 46 - 47 * 48 29896466 49 + 50 STO 12 51 CLX 52 34 53 * 54 21765 |
55
+ 56 * 57 548058 58 + 59 STO 07 60 CLX 61 2 62 * 63 277 64 + 65 * 66 2559 67 - 68 STO 08 69 CLX 70 17 71 * 72 2 73 - 74 * 75 16 76 - 77 * 78 455 79 - 80 * 81 32042 |
82
- 83 * 84 19567206 85 + 86 STO 09 87 CLX 88 5 89 * 90 368 91 + 92 * 93 62768 94 + 95 * 96 3470640 97 + 98 STO 04 99 CLX 100 5 101 * 102 109 103 + 104 * 105 1630 106 - 107 * 108 19901 |
109
- 110 STO 05 111 CLX 112 14 113 * 114 13 115 + 116 * 117 96 118 + 119 * 120 831 121 - 122 * 123 57375 124 - 125 * 126 9655387 127 + 128 STO 06 129 CLX 130 17 131 * 132 101 133 - 134 * 135 1962 |
136
- 137 * 138 168762 139 + 140 * 141 710988 142 + 143 STO 01 144 CLX 145 56 146 * 147 CHS 148 240 149 + 150 * 151 3335 152 + 153 * 154 6441 155 - 156 E4 157 ST/ 01 158 ST/ 04 159 ST/ 05 160 ST/ 07 161 ST/ 08 162 ST/ 10 163 ST/ 11 |
164
/ 165 STO 02 166 CLX 167 17 168 * 169 35 170 - 171 * 172 911 173 - 174 * 175 5648 176 + 177 * 178 58568 179 + 180 * 181 5056269 182 + 183 E6 184 ST/ 06 185 ST/ 09 186 ST/ 12 187 / 188 STO 03 189 RCL 02 190 RCL 01 191 END |
( 358 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 2024/07/07 at 16h41 TT
2024.0707 ENTER^
16.41 XEQ "JSUN"
>>>> LJup = 62°7637
= R01
RDN BJup = -0°8015
= R02
RDN RJup = 5.028970 AU = R03
-And in registers R04 thru R12:
R04 = LSat =
343°9898 R07 = LUra
= 53°7376
R10 = LNep = 358°0521
R05 = bSat =
-1°9075
R08 = bUra = -0°2695
R11 = bNep = -1°2717
R06 = RSat = 9.683356
R09 = RUra =
19.582836 R12 = RNep = 29.898921
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/