hp41programs

Astro5 Astronomical Ephemeris (VIII) for the HP41CV
 

Overview
 

1°)  HP41CV  without Module

 a)  Program1 [1900 , 2100 ]  ( 2116 bytes )

 b)  Program2 [ 1900/03/01 , 2100/02/28 ]  ( 1958 bytes )
 c)
 Program3 [ 1900/03/01 , 2100/02/28 ]  ( 1834 bytes )
 d)  Program4 [ 1900/03/01 , 2100/02/28 ]  ( 1517 bytes )

2°)  HP41CV  with X-Functions Module
- Program ( 867 bytes + 1516 or 1512 bytes )



1°)  HP41CV  Without Module
 

      a)  Program1  [ 1900 , 2100 ]  ( 2116 bytes )


-With geocentric coordinates, nutation, aberration and time-light are taken into account ( approximately )
-Refraction is computed with Bennet formula.


Registers
 

 R00 = Time from 2000/01/01 0h TT ( unit = 100 Julian years )            ( Registers R15-R16 are to be initialized before executing "EPH" )

 R01 = X   Rectangular
 R02 = Y   coordinates of the Sun,

 R03 = Geocentric longitude ( deg )
 R04 = Geocentric latitude ( deg )
 R05 = Distance Earth-Planet ( AU )

 R06 = Right-ascension ( deg )
 R07 = Declination ( deg )

 R08 = Heliocentric longitude ( deg )
 R09 = Heliocentric latitude ( deg )
 R10 = Radius vector ( AU )

 R11 = Azimuth ( deg )  ( reckoned clockwise from South )
 R12 = height ( deg )

 R13: temp

 R14 = local true sidereal time ( deg )
                            

•  R15 = Longitude of the observer ( deg )  positive East
•  R16 = Latitude of the observer ( deg)

Flags: /
Subroutines:  /



 01 LBL "EPH"
 02 19
 03 10^X
 04 SQRT
 05 /
 06 STO 00
 07 RDN
 08 HR
 09 STO 14
 10 24
 11 /
 12 X<>Y
 13 INT
 14 ST- L
 15 LASTX
 16 2
 17 10^X
 18 STO 02
 19 *
 20 ENTER
 21 INT
 22 STO 01
 23 -
 24 RCL 02
 25 *
 26 X<> 01
 27 1
 28 E^X
 29 -
 30 12
 31 /
 32 +
 33 INT
 34 ST- 01
 35 LASTX
 36 367
 37 *
 38 INT
 39 ST+ 01
 40 SIGN
 41 %
 42 INT
 43 .75
 44 CHS
 45 ST* Z
 46 *
 47 X<>Y
 48 RCL 01
 49 +
 50 INT
 51 +
 52 INT
 53 730434
 54 -
 55 +
 56 36525
 57 /
 58 ST+ 00
 59 36000770
 60 STO 06
 61 *
 62 3
 63 10^X
 64 STO 09
 65 /
 66 RCL 02
 67 +
 68 31
 69 1/X
 70 -
 71 RCL 14
 72 15
 73 *
 74 +
 75 RCL 15
 76 +
 77 STO 14
 78 1671
 79 STO 08
 80 48
 81 +
 82 STO 07
 83 4
 84 CHS
 85 STO 13
 86 80026
 87 CHS
 88 STO 12
 89 2963
 90 +
 91 STO 10
 92 RCL 00
 93 445267
 94 *
 95 68
 96 -
 97 STO 03
 98 SIN
 99 RCL 00
100 20
101 *
102 18
103 -
104 COS
105 -
106 RCL 00
107 32964
108 *
109 22
110 -
111 COS
112 -
113 ST+ X
114 X<> 03
115 COS
116 PI
117 %
118 STO 05
119 CLST
120 STO 01
121 STO 02
122 STO 04
123 STO 11
124 XEQ 14
125 CLX
126 RCL 13
127 RCL 05
128 P-R
129 STO 01
130 RDN
131 STO 02
132 CLX
133 RCL 11
134 RTN
135 LBL 03
136 RCL 00
137 3599905
138 %
139 2.96
140 -
141 STO 11
142 STO 12
143 RCL 00
144 44526711
145 %
146 68.25
147 -
148 STO 13
149 ST+ X
150 STO 05
151 STO 10
152 ST- 11
153 STO 06
154 ST+ X
155 STO 03
156 RCL 00
157 47719887
158 %
159 128.43
160 +
161 ST- 05
162 STO 07
163 ST+ 10
164 ST+ X
165 STO 08
166 RCL 00
167 48320202
168 %
169 86.66
170 +
171 STO 09
172 RCL 07
173 +
174 SIN
175 281
176 *
177 RCL 07
178 RCL 09
179 -
180 SIN
181 278
182 *
183 +
184 RCL 06
185 RCL 09
186 -
187 SIN
188 173
189 *
190 +
191 RCL 05
192 RCL 09
193 +
194 SIN
195 55
196 *
197 +
198 RCL 05
199 RCL 09          
200 -
201 SIN
202 46
203 *
204 +
205 RCL 06
206 RCL 09
207 +
208 SIN
209 33
210 *
211 +
212 RCL 08
213 RCL 09
214 +
215 SIN
216 17
217 *
218 +
219 RCL 10
220 RCL 09
221 -
222 SIN
223 RCL 08
224 RCL 09
225 -
226 SIN
227 +
228 9
229 *
230 +
231 RCL 09
232 RCL 11
233 +
234 SIN
235 RCL 09
236 SIN
237 641
238 *
239 -
240 8
241 *
242 -
243 RCL 09
244 RCL 10
245 +
246 SIN
247 4
248 *
249 +
250 STO 04
251 RCL 07
252 SIN
253 6289
254 *
255 RCL 06
256 SIN
257 658
258 *
259 +
260 RCL 08
261 SIN
262 214
263 *
264 +
265 RCL 09
266 ST+ X
267 STO 09
268 SIN
269 114
270 *
271 -
272 RCL 06
273 RCL 08
274 -
275 SIN
276 59
277 *
278 +
279 RCL 05
280 RCL 12
281 -
282 SIN
283 57
284 *
285 +
286 RCL 10
287 SIN
288 53
289 *
290 +
291 RCL 11
292 SIN
293 46
294 *
295 -
296 RCL 12
297 RCL 07
298 -
299 SIN
300 41
301 *
302 -
303 RCL 05
304 SIN
305 98
306 *
307 RCL 07
308 RCL 09
309 +
310 SIN
311 -
312 13
313 *
314 +
315 RCL 07
316 RCL 09
317 -
318 SIN
319 RCL 03
320 RCL 07
321 -
322 SIN
323 +
324 11
325 *
326 +
327 RCL 12
328 SIN
329 37
330 *
331 RCL 13
332 SIN
333 7
334 *
335 +
336 RCL 07
337 RCL 12
338 +
339 SIN
340 6
341 *
342 +
343 RCL 06
344 RCL 09
345 -
346 SIN
347 3
348 *
349 -
350 RCL 07
351 RCL 08
352 +
353 SIN
354 ST+ X
355 -
356 RCL 13
357 RCL 07
358 -
359 SIN
360 +
361 RCL 12
362 ST- 10
363 RCL 13
364 +
365 SIN
366 -
367 5
368 *
369 -
370 RCL 05
371 ST+ X
372 SIN
373 9
374 *
375 +
376 RCL 05
377 RCL 12
378 +
379 SIN
380 8
381 *
382 -
383 RCL 06
384 RCL 12
385 +
386 SIN
387 7
388 *
389 -
390 RCL 10
391 SIN
392 RCL 06
393 RCL 08
394 +
395 SIN
396 +
397 RCL 03
398 SIN
399 +
400 RCL 05
401 RCL 08
402 -
403 SIN
404 +
405 RCL 00          
406 132
407 *
408 60
409 -
410 SIN
411 -
412 4
413 *
414 +
415 RCL 08
416 RCL 12
417 -
418 SIN
419 RCL 05
420 RCL 09
421 +
422 SIN
423 -
424 3
425 *
426 +
427 481267881
428 RCL 00
429 ST+ X
430 -
431 RCL 00
432 *
433 211728
434 +
435 +
436 3
437 10^X
438 ST/ 04
439 /
440 360
441 MOD
442 STO 03
443 2574
444 RCL 07
445 COS
446 28
447 *
448 RCL 05
449 COS
450 5
451 *
452 +
453 RCL 06
454 COS
455 4
456 *
457 +
458 5
459 *
460 -
461 6
462 10^X
463 /
464 STO 05
465 XEQ 00
466 RTN
467 LBL 01
468 387.1
469 STO 09
470 20563
471 STO 08
472 7005
473 RCL 00
474 2
475 STO 13
476 *
477 +
478 STO 11
479 1556
480 STO 07
481 77456
482 STO 10
483 250205
484 STO 12
485 149474072
486 STO 06
487 48331
488 1186
489 XEQ 13
490 RTN
491 LBL 02
492 723.33
493 STO 09
494 677
495 STO 08
496 5
497 CHS
498 STO 13
499 3395
500 RCL 00
501 +
502 STO 11
503 1402
504 STO 07
505 131564
506 STO 10
507 181179
508 STO 12
509 58519213
510 STO 06
511 76680
512 901
513 XEQ 13
514 RTN
515 LBL 04
516 XEQ 12
517 RCL 04
518 10
519 *
520 STO 05
521 21
522 -
523 COS
524 14
525 *
526 RCL 09
527 22
528 *
529 RCL 03
530 +
531 ST+ 06
532 12
533 -
534 COS
535 6
536 *
537 -
538 RCL 06
539 49
540 +
541 COS
542 7
543 *
544 +
545 RCL 08
546 RCL 11
547 +
548 RCL 12
549 +
550 22
551 +
552 COS
553 RCL 06
554 ST+ X
555 12
556 +
557 COS
558 -
559 4
560 *
561 +
562 RCL 05
563 SIN
564 RCL 00
565 *
566 PI
567 *
568 +
569 RCL 04
570 7077
571 *
572 SIN
573 ST+ X
574 -
575 STO 03
576 CLX
577 STO 04
578 STO 05
579 1523.7
580 STO 09
581 9340
582 STO 08
583 1850
584 STO 11
585 9
586 STO 13
587 -
588 STO 07
589 4829
590 CHS
591 STO 12
592 19141696
593 STO 06
594 23940
595 CHS
596 STO 10
597 49558
598 772
599 XEQ 14
600 RTN
601 LBL 05
602 XEQ 12
603 RCL 08
604 73
605 +
606 COS
607 RCL 10
608 RCL 07          
609 -
610 41
611 +
612 SIN
613 -
614 509
615 SQRT
616 *
617 RCL 08
618 3
619 *
620 28
621 +
622 COS
623 RCL 09
624 ST+ X
625 STO 12
626 25
627 +
628 COS
629 +
630 5
631 *
632 -
633 RCL 06
634 RCL 12
635 +
636 69
637 -
638 COS
639 RCL 06
640 RCL 11
641 +
642 STO 04
643 12
644 +
645 SIN
646 4
647 *
648 -
649 RCL 07
650 COS
651 +
652 RCL 06
653 ST+ X
654 RCL 09
655 +
656 79
657 +
658 COS
659 -
660 RCL 11
661 ST- 06
662 32
663 -
664 COS
665 RCL 00
666 *
667 +
668 4
669 *
670 +
671 RCL 09
672 75
673 -
674 SIN
675 37
676 *
677 +
678 RCL 06
679 61
680 -
681 SIN
682 6
683 *
684 RCL 10
685 33
686 -
687 SIN
688 8
689 *
690 +
691 RCL 11
692 7.3
693 -
694 SIN
695 47
696 *
697 -
698 7
699 *
700 +
701 STO 03
702 RCL 06
703 63
704 -
705 SIN
706 5
707 *
708 STO 04
709 CLX
710 STO 05
711 5203
712 STO 09
713 1303
714 RCL 00
715 5
716 *
717 -
718 STO 11
719 3036303
720 STO 06
721 4849
722 STO 08
723 16
724 STO 13
725 1613
726 STO 07
727 34310
728 STO 12
729 14331
730 STO 10
731 86133
732 +
733 1021
734 XEQ 14
735 RTN
736 LBL 06
737 XEQ 12
738 CLX
739 21
740 -
741 COS
742 RCL 11
743 ST+ X
744 STO 12
745 75
746 +
747 COS
748 -
749 RCL 05
750 R^
751 STO 05
752 +
753 SIN
754 -
755 8
756 *
757 RCL 08
758 30
759 -
760 COS
761 RCL 10
762 54
763 +
764 COS
765 +
766 RCL 09
767 RCL 11
768 +
769 STO 10
770 RCL 06
771 -
772 25
773 -
774 COS
775 +
776 RCL 07
777 RCL 11
778 +
779 STO 04
780 50
781 -
782 SIN
783 5
784 *
785 +
786 RCL 11
787 7.3
788 -
789 SIN
790 90
791 *
792 +
793 9
794 *
795 +
796 RCL 07
797 RCL 03
798 -
799 COS
800 RCL 08
801 RCL 13
802 -
803 RCL 05
804 -
805 79
806 -
807 COS
808 +
809 RCL 09
810 ST+ X
811 ST+ 07
812 30
813 +
814 STO 12
815 COS
816 38
817 *
818 -
819 6
820 *
821 -
822 RCL 10          
823 SIN
824 7
825 *
826 RCL 11
827 32
828 -
829 COS
830 5
831 *
832 RCL 00
833 *
834 -
835 RCL 09
836 14
837 +
838 COS
839 59
840 *
841 +
842 ST+ X
843 +
844 STO 03
845 RCL 12
846 3
847 +
848 COS
849 20
850 ST+ 04
851 ST- 08
852 *
853 RCL 07
854 28
855 -
856 COS
857 RCL 04
858 COS
859 3
860 *
861 -
862 6
863 *
864 +
865 RCL 11
866 SIN
867 4
868 *
869 +
870 STO 04
871 RCL 12
872 SIN
873 19
874 *
875 RCL 08
876 COS
877 8
878 *
879 +
880 STO 05
881 1223511
882 STO 06
883 9543
884 STO 09
885 5553
886 STO 08
887 35
888 CHS
889 STO 13
890 2489
891 RCL 00
892 4
893 *
894 -
895 STO 11
896 1964
897 STO 07
898 50061
899 STO 12
900 93057
901 STO 10
902 20609
903 +
904 877
905 XEQ 14
906 RTN
907 LBL 07
908 XEQ 12
909 5
910 *
911 CHS
912 R^
913 STO 04
914 RCL 05
915 +
916 STO 08
917 10
918 STO 10
919 +
920 SIN
921 7
922 *
923 +
924 RCL 12
925 ST- 07
926 RCL 13
927 ST+ X
928 STO 11
929 +
930 20
931 -
932 SIN
933 ST+ X
934 RCL 07
935 28
936 +
937 COS
938 +
939 6
940 *
941 +
942 RCL 13
943 RCL 10
944 +
945 SIN
946 RCL 10
947 ST- 06
948 *
949 -
950 RCL 06
951 RCL 12
952 -
953 COS
954 RCL 12
955 RCL 04
956 -
957 COS
958 -
959 15
960 *
961 +
962 RCL 12
963 RCL 03
964 -
965 51
966 -
967 COS
968 ST+ Y
969 RCL 05
970 RCL 12
971 +
972 STO 08
973 11
974 -
975 COS
976 +
977 RCL 05
978 62.2
979 -
980 SIN
981 22
982 *
983 +
984 39
985 ST+ 04
986 *
987 +
988 RCL 03
989 73
990 +
991 SIN
992 RCL 04
993 SIN
994 +
995 35
996 *
997 -
998 RCL 11
999 19
1000 +
1001 STO 11
1002 SIN
1003 209
1004 *
1005 -
1006 STO 03
1007 RCL 11
1008 14
1009 +
1010 COS
1011 RCL 08
1012 SIN
1013 -
1014 6
1015 *
1016 STO 04
1017 RCL 11        
1018 COS
1019 34
1020 *
1021 STO 05
1022 429864
1023 STO 06
1024 19192
1025 STO 09
1026 45955
1027 CHS
1028 STO 12
1029 4630
1030 STO 08
1031 3
1032 CHS
1033 STO 13
1034 773
1035 RCL 00
1036 +
1037 STO 11
1038 1486
1039 STO 07
1040 172993
1041 STO 10
1042 74005
1043 521
1044 XEQ 14
1045 RTN
1046 LBL 08
1047 XEQ 12
1048 22
1049 *
1050 R^
1051 STO 04
1052 50
1053 -
1054 COS
1055 24
1056 *
1057 +
1058 RCL 06
1059 RCL 13
1060 ST+ 04
1061 RCL 05
1062 +
1063 ST- 07
1064 -
1065 COS
1066 9
1067 *
1068 +
1069 RCL 04
1070 77
1071 -
1072 COS
1073 RCL 07
1074 16
1075 +
1076 COS
1077 +
1078 5
1079 *
1080 +
1081 RCL 13
1082 7
1083 +
1084 SIN
1085 71
1086 *
1087 +
1088 RCL 05
1089 27.8
1090 +
1091 COS
1092 584
1093 *
1094 +
1095 STO 03
1096 RCL 04
1097 69
1098 -
1099 SIN
1100 RCL 13
1101 36
1102 -
1103 COS
1104 -
1105 9
1106 *
1107 STO 04
1108 RCL 13
1109 6
1110 +
1111 COS
1112 17
1113 *
1114 CHS
1115 STO 05
1116 219883
1117 STO 06
1118 30070
1119 STO 09
1120 899
1121 STO 08
1122 CLX
1123 STO 13
1124 1770
1125 RCL 00
1126 9
1127 *
1128 -
1129 STO 11
1130 1426
1131 STO 07
1132 55654
1133 CHS
1134 STO 12
1135 48124
1136 STO 10
1137 83660
1138 +
1139 1102
1140 XEQ 14
1141 RTN
1142 LBL 09
1143 RCL 00
1144 2890
1145 *
1146 65
1147 +
1148 COS
1149 7
1150 *
1151 STO 03
1152 RCL 00
1153 STO 05
1154 ST+ 05
1155 9
1156 *
1157 4
1158 -
1159 STO 04
1160 39489
1161 STO 09
1162 24900
1163 STO 08
1164 CLX
1165 STO 13
1166 17140
1167 STO 11
1168 224075
1169 STO 10
1170 238927
1171 STO 12
1172 146567
1173 STO 06
1174 36260
1175 -
1176 1397
1177 STO 07
1178 GTO 14
1179 LBL 12
1180 102
1181 RCL 00
1182 2.06
1183 *
1184 STO 04
1185 *
1186 STO 13
1187 4
1188 LASTX
1189 *
1190 STO 05
1191 52
1192 *
1193 STO 12
1194 19
1195 RCL 04
1196 *
1197 STO 11
1198 593
1199 LASTX
1200 *
1201 STO 07
1202 LASTX
1203 880
1204 *
1205 STO 08
1206 +
1207 STO 06
1208 LASTX
1209 RCL 07
1210 -
1211 STO 09
1212 +
1213 STO 10
1214 RCL 05
1215 ST+ X
1216 RCL 04
1217 31
1218 *
1219 STO 03
1220 RCL 13
1221 RCL 05
1222 -
1223 11
1224 -
1225 COS
1226 RTN
1227 LBL 13
1228 0
1229 STO 03
1230 STO 04
1231 STO 05
1232 RDN
1233 LBL 14
1234 RCL 00        
1235 ST* 06
1236 ST* 07
1237 ST* 13
1238 *
1239 +
1240 X<> 13
1241 ST+ 08
1242 5
1243 10^X
1244 ST/ 08
1245 RCL 06
1246 RCL 12
1247 +
1248 RCL 07
1249 RCL 10
1250 +
1251 9
1252 STO 06
1253 SQRT
1254 10^X
1255 ST/ 03
1256 ST/ 04
1257 ST/ 05
1258 ST/ 09
1259 ST/ 11
1260 ST/ 13
1261 ST/ Z
1262 /
1263 STO 10
1264 -
1265 ENTER
1266 ENTER
1267 LBL 11
1268 SIN
1269 RCL 08
1270 R-D
1271 *
1272 +
1273 DSE 06
1274 GTO 11
1275 STO 07
1276 2
1277 /
1278 TAN
1279 RCL 08
1280 1
1281 ST+ Y
1282 RCL 08
1283 -
1284 /
1285 SQRT
1286 *
1287 ATAN
1288 ST+ X
1289 1
1290 RCL 07
1291 COS
1292 RCL 08
1293 *
1294 -
1295 RCL 09
1296 *
1297 RCL 05
1298 +
1299 X<> 10
1300 +
1301 RCL 13
1302 ST+ 03
1303 -
1304 COS
1305 RCL 11
1306 LASTX
1307 SIN
1308 P-R
1309 X<>Y
1310 ASIN
1311 RCL 04
1312 +
1313 STO 09
1314 X<> Z
1315 R-P
1316 CLX
1317 RCL 03
1318 +
1319 STO 08
1320 RCL 10
1321 XEQ 10
1322 RCL 02
1323 ST+ Z
1324 CLX
1325 RCL 01
1326 +
1327 XEQ 12
1328 STO 05
1329 RDN
1330 STO 03
1331 STO 13
1332 X<>Y
1333 STO 04
1334 RCL 10
1335 SQRT
1336 1/X
1337 RCL 10
1338 ST/ Y
1339 X^2
1340 1
1341 -
1342 RCL 05
1343 X^2
1344 ST+ Y
1345 ST+ X
1346 /
1347 ST* Y
1348 -
1349 RCL 05
1350 176
1351 /
1352 ST* Y
1353 +
1354 ST- 03
1355 LBL 00
1356 RCL 00
1357 1934
1358 *
1359 55
1360 +
1361 1
1362 P-R
1363 391
1364 /
1365 549.4
1366 SQRT
1367 -
1368 RCL 00
1369 77
1370 /
1371 +
1372 STO 12
1373 COS
1374 X<>Y
1375 209
1376 /
1377 *
1378 X<> 12
1379 RCL 04
1380 RCL 03
1381 LASTX
1382 -
1383 STO 03
1384 XEQ 13
1385 STO 06
1386 RCL 16
1387 RCL Z
1388 STO 07
1389 90
1390 STO 11
1391 ST- Z
1392 R^
1393 +
1394 RCL 12
1395 +
1396 RCL 14
1397 -
1398 XEQ 13
1399 ST- 11
1400 X<>Y
1401 ENTER
1402 COS
1403 RCL 05
1404 409
1405 *
1406 /
1407 -
1408 STO 12
1409 10.3
1410 RCL 12        
1411 5.11
1412 +
1413 /
1414 +
1415 TAN
1416 1/X
1417 59
1418 /
1419 +
1420 X<>Y
1421 RCL 11
1422 RTN
1423 LBL 13
1424 1
1425 XEQ 10
1426 RDN
1427 R-P
1428 X<> Z
1429 ST- Y
1430 X<> Z
1431 P-R
1432 R^
1433 XEQ 12
1434 RDN
1435 RTN
1436 LBL 10
1437 X<>Y
1438 RDN
1439 P-R
1440 R^
1441 X<>Y
1442 P-R
1443 RTN
1444 LBL 12
1445 R-P
1446 X<>Y
1447 RDN
1448 R-P
1449 R^
1450 X<>Y
1451 END
 
 
    ( 2116 bytes / SIZE 017 )
 

    STACK           INPUTS                       OUTPUTS
         Z      YYYY.MNDD             h0 = apparent altitude  ( deg )
         Y      HH.MNSS (UT)                h  = true altitude  ( deg )
         X    DT = TT - UT ( sec)                 Az = Azimuth  ( deg )
 

   YYYY.MNDD   ENTER^
       HH.MNSS     ENTER^
        TT - UT        XEQ "EPH"   returns the azimuthal coordinates of the Sun.

-Then  R/S  or  XEQ 03   gives the coordinates of the Moon
-Then  R/S  or  XEQ 01   ------------------------  Mercury
     ............................................................................
-Then  R/S  or  XEQ 09   ------------------------- Pluto.


Example:         Seattle    Longitude = -122°333   Latitude = 47°600

   -122.333  STO 15
       47.6      STO 16

   Date      =  2024/05/01
   Time      =  12h41m  UT
 TT - UT  = 69.2 seconds


THE SUN
 

   2024.0501  ENTER^
      12.41       ENTER^
       69.2        XEQ "EPH"


               Az =   244°041   = R11                                                      
  RDN      h  =    -2°383    = R12
  RDN     h0 =    -1°687                        ( often meaningless when < 0 )

and   R01 = X = 0.75315 AU
         R02 = Y = 0.66949 AU

  R03 = 41°628 = Geocentric longitude                 R06 = 39.193 = right-ascension
  R04 = 0  = Geocentric latitude                            R07 = 15.321  = declination
  R05 = 1.00770 AU = distance Earth-Sun


-And  R14 =  local true sidereal time = 8927°8538 = 287°8538 = 19h1903 = 19h11m25s


THE MOON

   R/S  or  XEQ 03     >>>>   Az =  -27.342                                
                                  RDN      h  =  15.330
                                  RDN     h0 =  15.390
                  

  R03 = 312.289 = Geocentric longitude                R06  = -43.840 = right-ascension
  R04 =  -4.678   = Geocentric latitude                   R07 =  -21.598  = declination
  R05 = 0.002496 AU = distance Earth-Moon


MERCURY


   R/S  or  XEQ 01   >>>>    Az =   267.533                         
                                RDN      h  =    3.976
                                RDN     h0 =    4.165
                              

  R03 =  17.417 = Geocentric longitude                  R06 = 16.986 = right-ascension           R08 = -100.006 = heliocentric longitude
  R04 = -2.421  = Geocentric latitude                     R07 =  4.600  = declination                 R09 = -3.661      = heliocentric latitude
  R05 = 0.70519 AU = distance Earth-Mercury                                                                R10 = 0.46660  AU = radius vector


VENUS


   R/S  or  XEQ 02    >>>>     Az =   252.974                                        
                                 RDN      h  =   -0.239
                                 RDN     h0 =     0.278

  R03 = 32.524 = Geocentric longitude                 R06 =  30.754 = right-ascension          R08 = 19.828 = heliocentric longitude
  R04 = -1.214  = Geocentric latitude                   R07 = 11.208  = declination                 R09 =  -2.851  = heliocentric latitude
  R05 = 1.7019 AU = distance Earth-Venus                                                                     R10 =  0.7251 AU = radius vector


MARS


   R/S  or  XEQ 04   >>>>    Az =  -76.892
                                RDN      h  =    10.519
                                RDN     h0 =    10.604

  R03 =   0.679 = Geocentric longitude                   R06 = 1.126 = right-ascension            R08 = -27.879 = heliocentric longitude
  R04 =  -1.264   = Geocentric latitude                   R07 = -0.890  = declination                R09 = -1.807   = heliocentric latitude
  R05 = 1.9745 AU = distance Earth-Mars                                                                        R10 = 1.3817 AU = radius vector


JUPITER


   R/S  or  XEQ 05     >>>>    Az = 232.447                    
                                  RDN      h  = -7.534
                                  RDN     h0 = -7.615   ( often meaningless when < 0 )

  R03 = 54.249 = Geocentric longitude                   R06 = 52.076  = right-ascension             R08 = 56.778 = heliocentric longitude
  R04 = -0.759   = Geocentric latitude                    R07 =  18.097  = declination                   R09 = -0.908  = heliocentric latitude
  R05 =  5.989 AU = distance Earth-Jupiter                                                                           R10 = 5.010 AU = radius vector


SATURN


   R/S  or  XEQ 06     >>>>   Az = -62.988                                       
                                  RDN      h  =  13.916
                                  RDN     h0 =  13.982
  R03 = -13.321 = Geocentric longitude                   R06 = -11.570 = right-ascension            R08 = -18.194 = heliocentric longitude
  R04 = -1.752  = Geocentric latitude                      R07 =  -6.873  = declination                   R09 = -1.851  = heliocentric latitude
  R05 = 10.248 AU = distance Earth-Saturn                                                                          R10 = 9.704 AU = radius vector


URANUS


   R/S  or  XEQ 07   >>>>   Az =   234.010                                 
                                RDN      h  =  -6.436
                                RDN     h0 =  -6.503   ( often meaningless when < 0 )

  R03 = 52.415 = Geocentric longitude                    R06 =  50.077 = right-ascension             R08 =  52.975 = heliocentric longitude
  R04 = -0.264  = Geocentric latitude                       R07 =  18.119  = declination                   R09 = -0.277  = heliocentric latitude
  R05 = 20.568 AU = distance Earth-Uranus                                                                           R10 = 19.579 AU = radius vector


NEPTUNE


   R/S  or  XEQ 08   >>>>   Az =   -75.223                                
                                 RDN      h  =  11.070
                                 RDN     h0 =  11.152

  R03 = -1.050  = Geocentric longitude                    R06 =  -0.474 = right-ascension                R08 = -2.354 = heliocentric longitude
  R04 = -1.231  = Geocentric latitude                       R07 =  -1.547   = declination                    R09 = -1.261  = heliocentric latitude
  R05 = 30.624 AU = distance Earth-Neptune                                                                          R10 = 29.891 AU = radius vector


PLUTO


   R/S  or  XEQ 09     >>>>   Az = -16.725                                   
                                  RDN      h  =  18.064
                                  RDN     h0 =  18.115

  R03 = -57.901 = Geocentric longitude                    R06 =  -54.904 = right-ascension             R08 =  -59.529 = heliocentric longitude
  R04 = -3.028  = Geocentric latitude                       R07 =  -22.641  = declination                    R09 = -3.013  = heliocentric latitude
  R05 = 34.825 AU = distance Earth-Pluto                                                                                R10 = 35.006 AU = radius vector


Notes:

-HP41 in mode DEG
-Z output = h0 = apparent altitude but it is not in a register.


      b)  Program2  [ 1900/03/01 , 2100/02/28 ] ( 1958 bytes )


-Nutation, aberration and time-light are not computed
-Refraction is computed with Bennet formula.


Registers
 

 R00 = Time from 2000/01/01 0h TT ( unit = 100 Julian years )            ( Registers R15-R16 are to be initialized before executing "EPH" )

 R01 = X   Rectangular
 R02 = Y   coordinates of the Sun,

 R03 = Geocentric longitude ( deg )
 R04 = Geocentric latitude ( deg )
 R05 = Distance Earth-Planet ( AU )

 R06 = Right-ascension ( deg )
 R07 = Declination ( deg )

 R08 = Heliocentric longitude ( deg )
 R09 = Heliocentric latitude ( deg )
 R10 = Radius vector ( AU )

 R11 = Azimuth ( deg )  ( reckoned clockwise from South )
 R12 = height ( deg )

 R13: temp

 R14 = local sidereal time ( deg )
                            

•  R15 = Longitude of the observer ( deg )  positive East
•  R16 = Latitude of the observer ( deg)

Flags: /
Subroutines:  /



 01 LBL "EPH"
 02 19
 03 10^X
 04 SQRT
 05 /
 06 STO 00
 07 RDN
 08 HR
 09 STO 14
 10 24
 11 /
 12 X<>Y
 13 INT
 14 STO 01
 15 X<> L
 16 FRC
 17 2
 18 10^X
 19 STO 02
 20 STO 03
 21 *
 22 INT
 23 LASTX
 24 FRC
 25 RCL 02
 26 *
 27 X<>Y
 28 3
 29 +
 30 6
 31 X>Y?
 32 DSE 01
 33 INT
 34 X<=Y?
 35 CHS
 36 +
 37 4
 38 +
 39 30.6
 40 *
 41 INT
 42 +
 43 RCL 01
 44 365.25
 45 ST* 02
 46 *
 47 INT
 48 +
 49 730563
 50 -
 51 +
 52 RCL 02
 53 /
 54 ST+ 00
 55 36000770
 56 STO 06
 57 *
 58 3
 59 10^X
 60 STO 09
 61 /
 62 RCL 03
 63 +
 64 31
 65 1/X
 66 -
 67 RCL 14
 68 15
 69 *
 70 +
 71 RCL 15
 72 +
 73 STO 14
 74 1671
 75 STO 08
 76 48
 77 +
 78 STO 07
 79 4
 80 CHS
 81 STO 13
 82 80026
 83 CHS
 84 STO 12
 85 2963
 86 +
 87 STO 10
 88 RCL 00
 89 445267
 90 *
 91 68
 92 -
 93 STO 03
 94 SIN
 95 RCL 00
 96 20
 97 *
 98 18
 99 -
100 COS
101 -
102 RCL 00
103 32964
104 *
105 22
106 -
107 COS
108 -
109 ST+ X
110 X<> 03
111 COS
112 PI
113 %
114 STO 05
115 CLST
116 STO 01
117 STO 02
118 STO 04
119 STO 11
120 XEQ 14
121 CLX
122 RCL 03
123 RCL 05
124 P-R
125 STO 01
126 RDN
127 STO 02
128 CLX
129 RCL 11
130 RTN
131 LBL 03
132 RCL 00
133 3599905
134 %
135 2.96
136 -
137 STO 11
138 STO 12
139 RCL 00
140 44526711
141 %
142 68.25
143 -
144 STO 13
145 ST+ X
146 STO 05
147 STO 10
148 ST- 11
149 STO 06
150 ST+ X
151 STO 03
152 RCL 00
153 47719887
154 %
155 128.43
156 +
157 ST- 05
158 STO 07
159 ST+ 10
160 ST+ X
161 STO 08
162 RCL 00
163 48320202
164 %
165 86.66
166 +
167 STO 09
168 RCL 07
169 +
170 SIN
171 281
172 *
173 RCL 07
174 RCL 09
175 -
176 SIN
177 278
178 *
179 +
180 RCL 06
181 RCL 09          
182 -
183 SIN
184 173
185 *
186 +
187 RCL 05
188 RCL 09
189 +
190 SIN
191 55
192 *
193 +
194 RCL 05
195 RCL 09
196 -
197 SIN
198 46
199 *
200 +
201 RCL 06
202 RCL 09
203 +
204 SIN
205 33
206 *
207 +
208 RCL 08
209 RCL 09
210 +
211 SIN
212 17
213 *
214 +
215 RCL 10
216 RCL 09
217 -
218 SIN
219 RCL 08
220 RCL 09
221 -
222 SIN
223 +
224 9
225 *
226 +
227 RCL 09
228 RCL 11
229 +
230 SIN
231 RCL 09
232 SIN
233 641
234 *
235 -
236 8
237 *
238 -
239 STO 04
240 RCL 07
241 SIN
242 6289
243 *
244 RCL 06
245 SIN
246 658
247 *
248 +
249 RCL 08
250 SIN
251 214
252 *
253 +
254 RCL 09
255 ST+ X
256 STO 09
257 SIN
258 114
259 *
260 -
261 RCL 06
262 RCL 08
263 -
264 SIN
265 59
266 *
267 +
268 RCL 05
269 RCL 12
270 -
271 SIN
272 57
273 *
274 +
275 RCL 10
276 SIN
277 53
278 *
279 +
280 RCL 11
281 SIN
282 46
283 *
284 -
285 RCL 12
286 RCL 07
287 -
288 SIN
289 41
290 *
291 -
292 RCL 05
293 SIN
294 98
295 *
296 RCL 07
297 RCL 09
298 +
299 SIN
300 -
301 13
302 *
303 +
304 RCL 07
305 RCL 09
306 -
307 SIN
308 RCL 03
309 RCL 07
310 -
311 SIN
312 +
313 11
314 *
315 +
316 RCL 12
317 SIN
318 37
319 *
320 RCL 13
321 SIN
322 7
323 *
324 +
325 RCL 07
326 RCL 12
327 +
328 SIN
329 6
330 *
331 +
332 RCL 06
333 RCL 09
334 -
335 SIN
336 3
337 *
338 -
339 RCL 07
340 RCL 08
341 +
342 SIN
343 ST+ X
344 -
345 RCL 13
346 RCL 07
347 -
348 SIN
349 +
350 RCL 12
351 ST- 10
352 RCL 13
353 +
354 SIN
355 -
356 5
357 *
358 -
359 RCL 05
360 ST+ X
361 SIN
362 9
363 *
364 +
365 RCL 05
366 RCL 12
367 +
368 SIN
369 8
370 *
371 -
372 RCL 06          
373 RCL 12
374 +
375 SIN
376 7
377 *
378 -
379 RCL 10
380 SIN
381 RCL 06
382 RCL 08
383 +
384 SIN
385 +
386 RCL 03
387 SIN
388 +
389 RCL 05
390 RCL 08
391 -
392 SIN
393 +
394 4
395 *
396 +
397 481267881
398 RCL 00
399 ST+ X
400 -
401 RCL 00
402 *
403 211728
404 +
405 +
406 3
407 10^X
408 ST/ 04
409 /
410 360
411 MOD
412 STO 03
413 2574
414 RCL 07
415 COS
416 28
417 *
418 RCL 05
419 COS
420 5
421 *
422 +
423 RCL 06
424 COS
425 4
426 *
427 +
428 5
429 *
430 -
431 6
432 10^X
433 /
434 STO 05
435 XEQ 00
436 RTN
437 LBL 01
438 387.1
439 STO 09
440 20563
441 STO 08
442 7005
443 RCL 00
444 2
445 STO 13
446 *
447 +
448 STO 11
449 1556
450 STO 07
451 77456
452 STO 10
453 250205
454 STO 12
455 149474072
456 STO 06
457 48331
458 1186
459 XEQ 13
460 RTN
461 LBL 02
462 723.33
463 STO 09
464 677
465 STO 08
466 5
467 CHS
468 STO 13
469 3395
470 RCL 00
471 +
472 STO 11
473 1402
474 STO 07
475 131564
476 STO 10
477 181179
478 STO 12
479 58519213
480 STO 06
481 76680
482 901
483 XEQ 13
484 RTN
485 LBL 04
486 XEQ 12
487 RCL 04
488 10
489 *
490 21
491 -
492 COS
493 14
494 *
495 RCL 09
496 22
497 *
498 RCL 03
499 +
500 ST+ 06
501 12
502 -
503 COS
504 6
505 *
506 -
507 RCL 06
508 49
509 +
510 COS
511 7
512 *
513 +
514 STO 03
515 CLX
516 STO 04
517 STO 05
518 1523.7
519 STO 09
520 9340
521 STO 08
522 1850
523 STO 11
524 9
525 STO 13
526 -
527 STO 07
528 4829
529 CHS
530 STO 12
531 19141696
532 STO 06
533 23940
534 CHS
535 STO 10
536 49558
537 772
538 XEQ 14
539 RTN
540 LBL 05
541 XEQ 12
542 RCL 08
543 73
544 +
545 COS
546 RCL 10
547 RCL 07
548 -
549 41
550 +
551 SIN
552 -
553 509
554 SQRT
555 *
556 RCL 08
557 3
558 *
559 28
560 +
561 COS
562 RCL 09          
563 ST+ X
564 STO 12
565 25
566 +
567 COS
568 +
569 5
570 *
571 -
572 RCL 06
573 RCL 12
574 +
575 69
576 -
577 COS
578 RCL 06
579 RCL 11
580 +
581 STO 04
582 12
583 +
584 SIN
585 4
586 *
587 -
588 RCL 07
589 COS
590 +
591 RCL 06
592 ST+ X
593 STO 12
594 RCL 09
595 +
596 79
597 +
598 COS
599 -
600 RCL 11
601 ST- 06
602 32
603 -
604 COS
605 RCL 00
606 *
607 +
608 4
609 *
610 +
611 RCL 09
612 75
613 -
614 SIN
615 37
616 *
617 +
618 RCL 06
619 61
620 -
621 SIN
622 6
623 *
624 RCL 10
625 33
626 -
627 SIN
628 8
629 *
630 +
631 RCL 11
632 7.3
633 -
634 SIN
635 47
636 *
637 -
638 7
639 *
640 +
641 STO 03
642 RCL 06
643 63
644 -
645 SIN
646 5
647 *
648 STO 04
649 CLX
650 STO 05
651 5203
652 STO 09
653 1303
654 RCL 00
655 5
656 *
657 -
658 STO 11
659 3036303
660 STO 06
661 4849
662 STO 08
663 16
664 STO 13
665 1613
666 STO 07
667 34310
668 STO 12
669 14331
670 STO 10
671 86133
672 +
673 1021
674 XEQ 14
675 RTN
676 LBL 06
677 XEQ 12
678 CLX
679 21
680 -
681 COS
682 RCL 11
683 ST+ X
684 STO 12
685 75
686 +
687 COS
688 -
689 RCL 05
690 R^
691 STO 05
692 +
693 SIN
694 -
695 8
696 *
697 RCL 08
698 30
699 -
700 COS
701 RCL 10
702 54
703 +
704 COS
705 +
706 RCL 09
707 RCL 11
708 +
709 STO 10
710 RCL 06
711 -
712 25
713 -
714 COS
715 +
716 RCL 07
717 RCL 11
718 +
719 STO 04
720 50
721 -
722 SIN
723 5
724 *
725 +
726 RCL 11
727 7.3
728 -
729 SIN
730 90
731 *
732 +
733 9
734 *
735 +
736 RCL 07
737 RCL 03
738 -
739 COS
740 RCL 08
741 RCL 13
742 -
743 RCL 05
744 -
745 79
746 -
747 COS
748 +
749 RCL 09          
750 ST+ X
751 ST+ 07
752 30
753 +
754 STO 12
755 COS
756 38
757 *
758 -
759 6
760 *
761 -
762 RCL 10
763 SIN
764 7
765 *
766 RCL 11
767 32
768 -
769 COS
770 5
771 *
772 RCL 00
773 *
774 -
775 RCL 09
776 14
777 +
778 COS
779 59
780 *
781 +
782 ST+ X
783 +
784 STO 03
785 RCL 12
786 3
787 +
788 COS
789 20
790 ST+ 04
791 ST- 08
792 *
793 RCL 07
794 28
795 -
796 COS
797 RCL 04
798 COS
799 3
800 *
801 -
802 6
803 *
804 +
805 RCL 11
806 SIN
807 4
808 *
809 +
810 STO 04
811 RCL 12
812 SIN
813 19
814 *
815 RCL 08
816 COS
817 8
818 *
819 +
820 STO 05
821 1223511
822 STO 06
823 9543
824 STO 09
825 5553
826 STO 08
827 35
828 CHS
829 STO 13
830 2489
831 RCL 00
832 4
833 *
834 -
835 STO 11
836 1964
837 STO 07
838 50061
839 STO 12
840 93057
841 STO 10
842 20609
843 +
844 877
845 XEQ 14
846 RTN
847 LBL 07
848 XEQ 12
849 5
850 *
851 CHS
852 R^
853 STO 04
854 RCL 05
855 +
856 STO 08
857 10
858 STO 10
859 +
860 SIN
861 7
862 *
863 +
864 RCL 12
865 ST- 07
866 RCL 13
867 ST+ X
868 STO 11
869 +
870 20
871 -
872 SIN
873 ST+ X
874 RCL 07
875 28
876 +
877 COS
878 +
879 6
880 *
881 +
882 RCL 13
883 RCL 10
884 +
885 SIN
886 RCL 10
887 ST- 06
888 *
889 -
890 RCL 06
891 RCL 12
892 -
893 COS
894 RCL 12
895 RCL 04
896 -
897 COS
898 -
899 15
900 *
901 +
902 RCL 12
903 RCL 03
904 -
905 51
906 -
907 COS
908 ST+ Y
909 RCL 05
910 RCL 12
911 +
912 STO 08
913 11
914 -
915 COS
916 +
917 RCL 05
918 62.2
919 -
920 SIN
921 22
922 *
923 +
924 39
925 ST+ 04
926 *
927 +
928 RCL 03
929 73
930 +
931 SIN
932 RCL 04
933 SIN
934 +
935 35
936 *
937 -
938 RCL 11          
939 19
940 +
941 STO 11
942 SIN
943 209
944 *
945 -
946 STO 03
947 RCL 11
948 14
949 +
950 COS
951 RCL 08
952 SIN
953 -
954 6
955 *
956 STO 04
957 RCL 11
958 COS
959 34
960 *
961 STO 05
962 429864
963 STO 06
964 19192
965 STO 09
966 45955
967 CHS
968 STO 12
969 4630
970 STO 08
971 3
972 CHS
973 STO 13
974 773
975 RCL 00
976 +
977 STO 11
978 1486
979 STO 07
980 172993
981 STO 10
982 74005
983 521
984 XEQ 14
985 RTN
986 LBL 08
987 XEQ 12
988 22
989 *
990 R^
991 STO 04
992 50
993 -
994 COS
995 24
996 *
997 +
998 RCL 06
999 RCL 13
1000 ST+ 04
1001 RCL 05
1002 +
1003 ST- 07
1004 -
1005 COS
1006 9
1007 *
1008 +
1009 RCL 04
1010 77
1011 -
1012 COS
1013 RCL 07
1014 16
1015 +
1016 COS
1017 +
1018 5
1019 *
1020 +
1021 RCL 13
1022 7
1023 +
1024 SIN
1025 71
1026 *
1027 +
1028 RCL 05
1029 27.8
1030 +
1031 COS
1032 584
1033 *
1034 +
1035 STO 03
1036 RCL 04
1037 69
1038 -
1039 SIN
1040 RCL 13
1041 36
1042 -
1043 COS
1044 -
1045 9
1046 *
1047 STO 04
1048 RCL 13
1049 6
1050 +
1051 COS
1052 17
1053 *
1054 CHS
1055 STO 05
1056 219883
1057 STO 06
1058 30070
1059 STO 09
1060 899
1061 STO 08
1062 CLX
1063 STO 13
1064 1770
1065 RCL 00
1066 9
1067 *
1068 -
1069 STO 11
1070 1426
1071 STO 07
1072 55654
1073 CHS
1074 STO 12
1075 48124
1076 STO 10
1077 83660
1078 +
1079 1102
1080 XEQ 14
1081 RTN
1082 LBL 09
1083 39489
1084 STO 09
1085 24900
1086 STO 08
1087 CLX
1088 STO 13
1089 17140
1090 STO 11
1091 224075
1092 STO 10
1093 238927
1094 STO 12
1095 146567
1096 STO 06
1097 36260
1098 -
1099 1397
1100 STO 07
1101 LBL 13
1102 0
1103 STO 03
1104 STO 04
1105 STO 05
1106 RDN
1107 LBL 14
1108 RCL 00
1109 ST* 06
1110 ST* 07
1111 ST* 13
1112 *
1113 +
1114 X<> 13
1115 ST+ 08
1116 5
1117 10^X
1118 ST/ 08
1119 RCL 06
1120 RCL 12
1121 +
1122 RCL 07        
1123 RCL 10
1124 +
1125 9
1126 STO 06
1127 SQRT
1128 10^X
1129 ST/ 03
1130 ST/ 04
1131 ST/ 05
1132 ST/ 09
1133 ST/ 11
1134 ST/ 13
1135 ST/ Z
1136 /
1137 STO 10
1138 -
1139 ENTER
1140 ENTER
1141 LBL 11
1142 SIN
1143 RCL 08
1144 R-D
1145 *
1146 +
1147 DSE 06
1148 GTO 11
1149 STO 07
1150 2
1151 /
1152 TAN
1153 RCL 08
1154 1
1155 ST+ Y
1156 RCL 08
1157 -
1158 /
1159 SQRT
1160 *
1161 ATAN
1162 ST+ X
1163 1
1164 RCL 07
1165 COS
1166 RCL 08
1167 *
1168 -
1169 RCL 09
1170 *
1171 RCL 05
1172 +
1173 X<> 10
1174 +
1175 RCL 13
1176 ST+ 03
1177 -
1178 COS
1179 RCL 11
1180 LASTX
1181 SIN
1182 P-R
1183 X<>Y
1184 ASIN
1185 RCL 04
1186 +
1187 STO 09
1188 X<> Z
1189 R-P
1190 CLX
1191 RCL 03
1192 +
1193 STO 08
1194 RCL 10
1195 XEQ 10
1196 RCL 02
1197 ST+ Z
1198 CLX
1199 RCL 01
1200 +
1201 XEQ 14
1202 STO 05
1203 RDN
1204 STO 03
1205 X<>Y
1206 STO 04
1207 LBL 00
1208 RCL 00
1209 77
1210 /
1211 549.4
1212 SQRT
1213 -
1214 RCL 04
1215 RCL 03
1216 XEQ 13
1217 STO 06
1218 RCL 16
1219 RCL Z
1220 STO 07
1221 90
1222 STO 11
1223 ST- Z
1224 R^
1225 +
1226 RCL 14
1227 -
1228 XEQ 13
1229 ST- 11
1230 X<>Y
1231 ENTER
1232 COS
1233 RCL 05
1234 409
1235 *
1236 /
1237 -
1238 STO 12
1239 10.3
1240 RCL 12
1241 5.11
1242 +
1243 /
1244 +
1245 TAN
1246 1/X
1247 59
1248 /
1249 +
1250 X<>Y
1251 RCL 11
1252 RTN
1253 LBL 13
1254 1
1255 XEQ 10
1256 RDN
1257 R-P
1258 X<> Z
1259 ST- Y
1260 X<> Z
1261 P-R
1262 R^
1263 XEQ 14
1264 RDN
1265 RTN
1266 LBL 10
1267 X<>Y
1268 RDN
1269 P-R
1270 R^
1271 X<>Y
1272 P-R
1273 RTN
1274 LBL 14
1275 R-P
1276 X<>Y
1277 RDN
1278 R-P
1279 R^
1280 X<>Y
1281 RTN
1282 LBL 12
1283 102
1284 RCL 00
1285 2.06
1286 *
1287 STO 04
1288 *
1289 STO 13
1290 4
1291 LASTX
1292 *
1293 STO 05
1294 52
1295 *
1296 STO 12
1297 19
1298 RCL 04
1299 *
1300 STO 11
1301 593
1302 LASTX
1303 *
1304 STO 07
1305 LASTX
1306 880
1307 *
1308 STO 08
1309 +
1310 STO 06
1311 LASTX
1312 RCL 07
1313 -
1314 STO 09
1315 +
1316 STO 10
1317 RCL 05
1318 ST+ X
1319 RCL 04        
1320 31
1321 *
1322 STO 03
1323 RCL 13
1324 RCL 05
1325 -
1326 11
1327 -
1328 COS
1329 END
 
 
    ( 1958 bytes / SIZE 017 )
 

    STACK           INPUTS                       OUTPUTS
         Z      YYYY.MNDD             h0 = apparent altitude  ( deg )
         Y      HH.MNSS (UT)                h  = true altitude  ( deg )
         X    DT = TT - UT ( sec)                 Az = Azimuth  ( deg )
 

   YYYY.MNDD   ENTER^
       HH.MNSS     ENTER^
        TT - UT        XEQ "EPH"   returns the azimuthal coordinates of the Sun.

-Then  R/S  or  XEQ 03   gives the coordinates of the Moon
-Then  R/S  or  XEQ 01   ------------------------  Mercury
     ............................................................................
-Then  R/S  or  XEQ 09   ------------------------- Pluto.


Example:         Seattle    Longitude = -122°333   Latitude = 47°600

   -122.333  STO 15
       47.6      STO 16

   Date      =  2024/05/01
   Time      =  12h41m  UT
 TT - UT  = 69.2 seconds


THE SUN
 

   2024.0501  ENTER^
      12.41       ENTER^
       69.2        XEQ "EPH"


               Az =   244°036   = R11                                                      
  RDN      h  =    -2°387    = R12
  RDN     h0 =    -1°691    = R13      ( often meaningless when < 0 )

and   R01 = X = 0.75315 AU
         R02 = Y = 0.66949 AU

  R03 = 41°635 = Geocentric longitude                 R06 = 39.200 = right-ascension
  R04 = 0  = Geocentric latitude                            R07 = 15.322  = declination
  R05 = 1.00770 AU = distance Earth-Sun


-And  R14 =  local true sidereal time = 8927°8538 = 287°8538 = 19h1903 = 19h11m25s


THE MOON

   R/S  or  XEQ 03     >>>>   Az =  -27.342                                
                                  RDN      h  =  15.328
                                  RDN     h0 =  15.388
                  

  R03 = 312.291 = Geocentric longitude                R06  = -43.838 = right-ascension
  R04 =  -4.682   = Geocentric latitude                   R07 =  -21.599  = declination
  R05 = 0.002496 AU = distance Earth-Moon


MERCURY


   R/S  or  XEQ 01   >>>>    Az =   267.532                         
                                RDN      h  =    3.975
                                RDN     h0 =    4.164
                              

  R03 =  17.420 = Geocentric longitude                  R06 = 16.988 = right-ascension           R08 = -100.006 = heliocentric longitude
  R04 = -2.421  = Geocentric latitude                     R07 =  4.600  = declination                 R09 = -3.661      = heliocentric latitude
  R05 = 0.70519 AU = distance Earth-Mercury                                                                R10 = 0.46660  AU = radius vector


VENUS


   R/S  or  XEQ 02    >>>>     Az =   252.963                                        
                                 RDN      h  =   -0.245
                                 RDN     h0 =     0.274

  R03 = 32.537 = Geocentric longitude                 R06 =  30.767 = right-ascension          R08 = 19.828 = heliocentric longitude
  R04 = -1.214  = Geocentric latitude                   R07 = 11.211  = declination                 R09 =  -2.851  = heliocentric latitude
  R05 = 1.7019 AU = distance Earth-Venus                                                                     R10 =  0.7251 AU = radius vector


MARS


   R/S  or  XEQ 04   >>>>    Az =  -76.902            
                                RDN      h  =    10.516
                                RDN     h0 =    10.602

  R03 =   0.690   = Geocentric longitude                 R06 = 1.136 = right-ascension            R08 = -27.876 = heliocentric longitude
  R04 =  -1.264   = Geocentric latitude                   R07 = -0.886  = declination                R09 = -1.807   = heliocentric latitude
  R05 = 1.9745 AU = distance Earth-Mars                                                                        R10 = 1.3817 AU = radius vector


JUPITER


   R/S  or  XEQ 05     >>>>    Az = 232.440
                                  RDN      h  = -7.538
                                  RDN     h0 = -7.619   ( often meaningless when < 0 )

  R03 = 54.259 = Geocentric longitude                   R06 = 52.086  = right-ascension             R08 = 56.778 = heliocentric longitude
  R04 = -0.759   = Geocentric latitude                    R07 =  18.097  = declination                   R09 = -0.908  = heliocentric latitude
  R05 =  5.989 AU = distance Earth-Jupiter                                                                           R10 = 5.010 AU = radius vector


SATURN


   R/S  or  XEQ 06     >>>>   Az = -62.993                                       
                                  RDN      h  =  13.916
                                  RDN     h0 =  13.982
  R03 = -13.315 = Geocentric longitude                   R06 = -11.565 = right-ascension            R08 = -18.194 = heliocentric longitude
  R04 = -1.752  = Geocentric latitude                      R07 =  -6.870  = declination                   R09 = -1.851  = heliocentric latitude
  R05 = 10.248 AU = distance Earth-Saturn                                                                          R10 = 9.704 AU = radius vector


URANUS


   R/S  or  XEQ 07   >>>>   Az =   234.003
                                RDN      h  =  -6.440
                                RDN     h0 =  -6.507   ( often meaningless when < 0 )

  R03 = 52.423 = Geocentric longitude                    R06 =  50.086 = right-ascension             R08 =  52.975 = heliocentric longitude
  R04 = -0.264  = Geocentric latitude                       R07 =  18.119  = declination                   R09 = -0.277  = heliocentric latitude
  R05 = 20.568 AU = distance Earth-Uranus                                                                           R10 = 19.579 AU = radius vector


NEPTUNE


   R/S  or  XEQ 08   >>>>   Az =   -75.228                                
                                 RDN      h  =  11.069
                                 RDN     h0 =  11.151

  R03 = -1.044  = Geocentric longitude                    R06 =  -0.468 = right-ascension                R08 = -2.354 = heliocentric longitude
  R04 = -1.231  = Geocentric latitude                       R07 =  -1.544   = declination                    R09 = -1.261  = heliocentric latitude
  R05 = 30.624 AU = distance Earth-Neptune                                                                          R10 = 29.891 AU = radius vector


PLUTO


   R/S  or  XEQ 09     >>>>   Az = -16.720                
                                  RDN      h  =  18.069
                                  RDN     h0 =  18.119

  R03 = -57.904 = Geocentric longitude                    R06 =  -54.908 = right-ascension             R08 =  -59.533 = heliocentric longitude
  R04 = -3.027  = Geocentric latitude                       R07 =  -22.638  = declination                    R09 = -3.011  = heliocentric latitude
  R05 = 34.825 AU = distance Earth-Pluto                                                                                R10 = 35.006 AU = radius vector



Notes:

-The number of days since 2000/01/01 are computed with less bytes than the first "EPH"
-But it doesn't work if the date is before 1900/03/01 or after 2100/02/28

-However, if the date is for instance:  1900/01/06  put  1900/01/07  in Z-register
-And with  2100/03/03  put  2100/03/02

-HP41 in mode DEG
-Z output = h0 = apparent altitude but it is not in a register.

-With this program, we can add another program with 154 bytes on HP41CV without module.


      c)  Program3  [ 1900/03/01 , 2100/02/28 ] ( 1834 bytes )


-Nutation, aberration and time-light are not computed
-Refraction is computed with Bennet formula.

Registers
 

 R00 = Time from 2000/01/01 0h TT ( unit = 100 Julian years )            ( Registers R14 & R15 are to be initialized before executing "EPH" )

 R01 = X   Rectangular
 R02 = Y   coordinates of the Sun,

 R03 = Geocentric longitude ( deg )
 R04 = Geocentric latitude ( deg )
 R05 = Distance Earth-Planet ( AU )

 R06 = Right-ascension ( deg )
 R07 = Declination ( deg )

 R08 = Heliocentric longitude ( deg )
 R09 = Heliocentric latitude ( deg )
 R10 = Radius vector ( AU )

 R11 = Azimuth ( deg )  ( reckoned clockwise from South )
 R12 = height ( deg )

 R13 = local sidereal time ( deg )

•  R14 = Longitude of the observer ( deg )  positive East
•  R15 = Latitude of the observer ( deg)

Flags: /
Subroutines:  /



 01 LBL "EPH"
 02 19
 03 10^X
 04 SQRT
 05 /
 06 STO 00
 07 RDN
 08 HR
 09 STO 13
 10 24
 11 /
 12 X<>Y
 13 INT
 14 STO 01
 15 X<> L
 16 FRC
 17 2
 18 10^X
 19 STO 02
 20 STO 03
 21 *
 22 INT
 23 LASTX
 24 FRC
 25 RCL 02
 26 *
 27 X<>Y
 28 3
 29 +
 30 6
 31 X>Y?
 32 DSE 01
 33 INT
 34 X<=Y?
 35 CHS
 36 +
 37 4
 38 +
 39 30.6
 40 *
 41 INT
 42 +
 43 RCL 01
 44 365.25
 45 ST* 02
 46 *
 47 INT
 48 +
 49 730563
 50 -
 51 +
 52 RCL 02
 53 /
 54 ST+ 00
 55 36000770
 56 STO 06
 57 *
 58 3
 59 10^X
 60 STO 09
 61 /
 62 RCL 03
 63 +
 64 31
 65 1/X
 66 -
 67 RCL 13
 68 15
 69 *
 70 +
 71 RCL 14
 72 +
 73 STO 13
 74 1719
 75 STO 07
 76 48
 77 -
 78 RCL 00
 79 4
 80 *
 81 -
 82 STO 08
 83 80026
 84 CHS
 85 STO 12
 86 2963
 87 +
 88 STO 10
 89 CLST
 90 STO 01
 91 STO 02
 92 STO 11
 93 XEQ 13
 94 CLX
 95 RCL 03
 96 RCL 05
 97 P-R
 98 STO 01
 99 RDN
100 STO 02
101 CLX
102 RCL 11
103 RTN
104 LBL 03
105 RCL 00
106 3599905
107 %
108 2.96
109 -
110 STO 11
111 STO 12
112 RCL 00
113 44526711
114 %
115 68.25
116 -
117 STO 03
118 ST+ X
119 STO 05
120 STO 10
121 ST- 11
122 STO 06
123 RCL 00
124 47719887
125 %
126 128.43
127 +
128 ST- 05
129 STO 07
130 ST+ 10
131 ST+ X
132 STO 08
133 RCL 00
134 48320202
135 %
136 86.66
137 +
138 STO 09
139 RCL 07
140 +
141 SIN
142 281
143 *
144 RCL 07
145 RCL 09
146 -
147 SIN
148 278
149 *
150 +
151 RCL 06
152 RCL 09
153 -
154 STO 04
155 SIN
156 173
157 *
158 +
159 RCL 05
160 RCL 09
161 +
162 SIN
163 55
164 *
165 +
166 RCL 05          
167 RCL 09
168 -
169 SIN
170 46
171 *
172 +
173 RCL 06
174 RCL 09
175 +
176 SIN
177 33
178 *
179 +
180 RCL 08
181 RCL 09
182 +
183 SIN
184 17
185 *
186 +
187 RCL 10
188 RCL 09
189 -
190 SIN
191 RCL 08
192 RCL 09
193 -
194 SIN
195 +
196 9
197 *
198 +
199 RCL 09
200 RCL 11
201 +
202 SIN
203 RCL 09
204 SIN
205 641
206 *
207 -
208 8
209 *
210 -
211 RCL 04
212 RCL 08
213 -
214 SIN
215 RCL 09
216 RCL 10
217 +
218 SIN
219 +
220 4
221 *
222 +
223 STO 04
224 RCL 07
225 SIN
226 6289
227 *
228 RCL 06
229 SIN
230 658
231 *
232 +
233 RCL 08
234 SIN
235 214
236 *
237 +
238 RCL 09
239 ST+ X
240 STO 09
241 SIN
242 114
243 *
244 -
245 RCL 06
246 RCL 08
247 -
248 SIN
249 59
250 *
251 +
252 RCL 05
253 RCL 12
254 -
255 SIN
256 57
257 *
258 +
259 RCL 10
260 SIN
261 53
262 *
263 +
264 RCL 11
265 SIN
266 46
267 *
268 -
269 RCL 12
270 RCL 07
271 -
272 SIN
273 41
274 *
275 -
276 RCL 05
277 SIN
278 98
279 *
280 RCL 07
281 RCL 09
282 +
283 SIN
284 -
285 13
286 *
287 +
288 RCL 07
289 RCL 09
290 -
291 SIN
292 RCL 06
293 ST+ X
294 STO 11
295 RCL 07
296 -
297 SIN
298 +
299 11
300 *
301 +
302 RCL 12
303 SIN
304 37
305 *
306 RCL 03
307 SIN
308 7
309 *
310 +
311 RCL 07
312 RCL 12
313 +
314 SIN
315 6
316 *
317 +
318 RCL 06
319 RCL 09
320 -
321 SIN
322 3
323 *
324 -
325 RCL 07
326 RCL 08
327 +
328 SIN
329 ST+ X
330 -
331 RCL 03
332 RCL 07
333 -
334 SIN
335 +
336 RCL 12
337 RCL 03
338 +
339 SIN
340 -
341 5
342 *
343 -
344 RCL 05          
345 ST+ X
346 SIN
347 9
348 *
349 +
350 RCL 05
351 RCL 12
352 +
353 SIN
354 8
355 *
356 -
357 RCL 06
358 RCL 12
359 +
360 SIN
361 7
362 *
363 -
364 481267881
365 RCL 00
366 *
367 211728
368 +
369 +
370 3
371 10^X
372 ST/ 04
373 /
374 360
375 MOD
376 STO 03
377 2574
378 RCL 07
379 COS
380 28
381 *
382 RCL 05
383 COS
384 5
385 *
386 +
387 RCL 06
388 COS
389 4
390 *
391 +
392 5
393 *
394 -
395 6
396 10^X
397 /
398 STO 05
399 XEQ 00
400 RTN
401 LBL 01
402 387.1
403 STO 09
404 20563
405 STO 08
406 7005
407 RCL 00
408 ST+ X
409 ST+ 08
410 +
411 STO 11
412 1556
413 STO 07
414 77456
415 STO 10
416 250205
417 STO 12
418 149474072
419 STO 06
420 48331
421 1186
422 XEQ 13
423 RTN
424 LBL 02
425 723.33
426 STO 09
427 677
428 RCL 00
429 5
430 *
431 -
432 STO 08
433 3395
434 RCL 00
435 +
436 STO 11
437 1402
438 STO 07
439 131564
440 STO 10
441 181179
442 STO 12
443 58519213
444 STO 06
445 76680
446 901
447 XEQ 13
448 RTN
449 LBL 04
450 XEQ 12
451 RCL 04
452 10
453 *
454 21
455 -
456 COS
457 14
458 *
459 RCL 09
460 22
461 *
462 RCL 03
463 +
464 ST+ 06
465 12
466 -
467 COS
468 6
469 *
470 -
471 RCL 06
472 49
473 +
474 COS
475 7
476 *
477 +
478 STO 03
479 CLX
480 STO 04
481 STO 05
482 1523.7
483 STO 09
484 9340
485 RCL 00
486 9
487 STO 11
488 *
489 +
490 STO 08
491 1841
492 STO 07
493 ST+ 11
494 4829
495 CHS
496 STO 12
497 19141696
498 STO 06
499 23940
500 CHS
501 STO 10
502 49558
503 772
504 XEQ 14
505 RTN
506 LBL 05
507 XEQ 12
508 RCL 08
509 73
510 +
511 COS
512 RCL 10
513 RCL 07
514 -
515 41
516 +
517 SIN
518 -
519 509
520 SQRT
521 *
522 RCL 06          
523 RCL 11
524 +
525 12
526 +
527 SIN
528 16
529 *
530 -
531 RCL 09
532 75
533 -
534 SIN
535 37
536 *
537 +
538 RCL 06
539 RCL 11
540 -
541 STO 06
542 61
543 -
544 SIN
545 6
546 *
547 RCL 10
548 33
549 -
550 SIN
551 8
552 *
553 +
554 RCL 11
555 7.3
556 -
557 SIN
558 47
559 *
560 -
561 7
562 *
563 +
564 STO 03
565 RCL 06
566 63
567 -
568 SIN
569 5
570 *
571 STO 04
572 CLX
573 STO 05
574 5203
575 STO 09
576 1303
577 RCL 00
578 5
579 *
580 -
581 STO 11
582 3036303
583 STO 06
584 4849
585 RCL 00
586 16
587 *
588 +
589 STO 08
590 1613
591 STO 07
592 34310
593 STO 12
594 14331
595 STO 10
596 86133
597 +
598 1021
599 XEQ 14
600 RTN
601 LBL 06
602 XEQ 12
603 STO 04
604 CLX
605 21
606 -
607 COS
608 RCL 11
609 ST+ X
610 STO 12
611 75
612 +
613 COS
614 -
615 RCL 05
616 ENTER
617 ST+ X
618 STO 05
619 +
620 SIN
621 -
622 8
623 *
624 RCL 08
625 30
626 -
627 COS
628 RCL 10
629 54
630 +
631 COS
632 +
633 RCL 09
634 RCL 11
635 +
636 STO 10
637 RCL 06
638 -
639 25
640 -
641 COS
642 +
643 RCL 07
644 RCL 11
645 +
646 STO 06
647 50
648 -
649 SIN
650 5
651 *
652 +
653 RCL 11
654 7.3
655 -
656 SIN
657 90
658 *
659 +
660 9
661 *
662 +
663 RCL 07
664 RCL 03
665 -
666 COS
667 RCL 08
668 RCL 04
669 -
670 RCL 05
671 -
672 79
673 -
674 COS
675 +
676 RCL 09
677 ST+ X
678 ST+ 07
679 30
680 +
681 STO 12
682 COS
683 38
684 *
685 -
686 6
687 *
688 -
689 RCL 10
690 SIN
691 7
692 *
693 RCL 11
694 32
695 -
696 COS
697 5
698 *
699 RCL 00          
700 *
701 -
702 RCL 09
703 14
704 +
705 COS
706 59
707 *
708 +
709 ST+ X
710 +
711 STO 03
712 RCL 12
713 3
714 +
715 COS
716 20
717 ST+ 06
718 ST- 08
719 *
720 RCL 07
721 28
722 -
723 COS
724 RCL 06
725 COS
726 3
727 *
728 -
729 6
730 *
731 +
732 STO 04
733 RCL 12
734 SIN
735 19
736 *
737 RCL 08
738 COS
739 8
740 *
741 +
742 STO 05
743 1223511
744 STO 06
745 9543
746 STO 09
747 5553
748 RCL 00
749 35
750 *
751 -
752 STO 08
753 2489
754 RCL 00
755 4
756 *
757 -
758 STO 11
759 1964
760 STO 07
761 50061
762 STO 12
763 93057
764 STO 10
765 20609
766 +
767 877
768 XEQ 14
769 RTN
770 LBL 07
771 XEQ 12
772 STO 09
773 RCL 05
774 ST+ X
775 STO 04
776 RCL 05
777 +
778 STO 08
779 10
780 STO 10
781 +
782 SIN
783 7
784 *
785 RCL 12
786 ST- 07
787 RCL 09
788 ST+ X
789 STO 11
790 +
791 20
792 -
793 SIN
794 ST+ X
795 RCL 07
796 28
797 +
798 COS
799 +
800 6
801 *
802 +
803 RCL 09
804 RCL 10
805 +
806 SIN
807 RCL 10
808 ST- 06
809 *
810 -
811 RCL 06
812 RCL 12
813 -
814 COS
815 RCL 12
816 RCL 04
817 -
818 COS
819 -
820 15
821 *
822 +
823 RCL 12
824 RCL 03
825 -
826 51
827 -
828 COS
829 ST+ Y
830 RCL 05
831 RCL 12
832 +
833 STO 08
834 11
835 -
836 COS
837 +
838 RCL 05
839 62.2
840 -
841 SIN
842 22
843 *
844 +
845 39
846 ST+ 04
847 *
848 +
849 RCL 03
850 73
851 +
852 SIN
853 RCL 04
854 SIN
855 +
856 35
857 *
858 -
859 RCL 11
860 19
861 +
862 STO 11
863 SIN
864 209
865 *
866 -
867 STO 03
868 RCL 11
869 14
870 +
871 COS
872 RCL 08          
873 SIN
874 -
875 6
876 *
877 STO 04
878 RCL 11
879 COS
880 34
881 *
882 STO 05
883 429864
884 STO 06
885 19192
886 STO 09
887 45955
888 CHS
889 STO 12
890 4630
891 RCL 00
892 3
893 *
894 -
895 STO 08
896 773
897 RCL 00
898 +
899 STO 11
900 1486
901 STO 07
902 172993
903 STO 10
904 74005
905 521
906 XEQ 14
907 RTN
908 LBL 08
909 XEQ 12
910 STO 10
911 RCL 05
912 -
913 11
914 -
915 COS
916 22
917 *
918 RCL 05
919 ST+ X
920 STO 04
921 50
922 -
923 COS
924 24
925 *
926 +
927 RCL 06
928 RCL 10
929 ST+ 04
930 RCL 05
931 +
932 -
933 COS
934 9
935 *
936 +
937 RCL 10
938 7
939 +
940 SIN
941 71
942 *
943 +
944 RCL 05
945 27.8
946 +
947 COS
948 584
949 *
950 +
951 STO 03
952 RCL 04
953 69
954 -
955 SIN
956 RCL 10
957 36
958 -
959 COS
960 -
961 9
962 *
963 STO 04
964 RCL 10
965 6
966 +
967 COS
968 17
969 *
970 CHS
971 STO 05
972 219883
973 STO 06
974 30070
975 STO 09
976 899
977 STO 08
978 1770
979 RCL 00
980 9
981 *
982 -
983 STO 11
984 1426
985 STO 07
986 55654
987 CHS
988 STO 12
989 48124
990 STO 10
991 83660
992 +
993 1102
994 XEQ 14
995 RTN
996 LBL 09
997 39489
998 STO 09
999 24900
1000 STO 08
1001 17140
1002 STO 11
1003 224075
1004 STO 10
1005 238927
1006 STO 12
1007 146567
1008 STO 06
1009 36260
1010 -
1011 1397
1012 STO 07
1013 LBL 13
1014 0
1015 STO 03
1016 STO 04
1017 STO 05
1018 RDN
1019 LBL 14
1020 RCL 00
1021 ST* 06
1022 ST* 07
1023 *
1024 +
1025 X<> 12
1026 5
1027 10^X
1028 ST/ 08
1029 CLX
1030 RCL 06
1031 +
1032 RCL 07
1033 RCL 10
1034 +
1035 9
1036 STO 06
1037 SQRT
1038 10^X
1039 ST/ 03
1040 ST/ 04
1041 ST/ 05
1042 ST/ 09
1043 ST/ 11
1044 ST/ 12
1045 ST/ Z
1046 /
1047 STO 10
1048 -
1049 ENTER
1050 ENTER
1051 LBL 11
1052 SIN
1053 RCL 08        
1054 R-D
1055 *
1056 +
1057 DSE 06
1058 GTO 11
1059 STO 07
1060 2
1061 /
1062 TAN
1063 RCL 08
1064 1
1065 ST+ Y
1066 RCL 08
1067 -
1068 /
1069 SQRT
1070 *
1071 ATAN
1072 ST+ X
1073 1
1074 RCL 07
1075 COS
1076 RCL 08
1077 *
1078 -
1079 RCL 09
1080 *
1081 RCL 05
1082 +
1083 X<> 10
1084 +
1085 RCL 12
1086 ST+ 03
1087 -
1088 COS
1089 RCL 11
1090 LASTX
1091 SIN
1092 P-R
1093 X<>Y
1094 ASIN
1095 RCL 04
1096 +
1097 STO 09
1098 X<> Z
1099 R-P
1100 CLX
1101 RCL 03
1102 +
1103 STO 08
1104 RCL 10
1105 XEQ 10
1106 RCL 02
1107 ST+ Z
1108 CLX
1109 RCL 01
1110 +
1111 XEQ 14
1112 STO 05
1113 RDN
1114 STO 03
1115 X<>Y
1116 STO 04
1117 LBL 00
1118 RCL 00
1119 77
1120 /
1121 549.4
1122 SQRT
1123 -
1124 RCL 04
1125 RCL 03
1126 XEQ 13
1127 STO 06
1128 RCL 15
1129 RCL Z
1130 STO 07
1131 90
1132 STO 11
1133 ST- Z
1134 R^
1135 +
1136 RCL 13
1137 -
1138 XEQ 13
1139 ST- 11
1140 X<>Y
1141 ENTER
1142 COS
1143 RCL 05
1144 409
1145 *
1146 /
1147 -
1148 STO 12
1149 10.3
1150 RCL 12
1151 5.11
1152 +
1153 /
1154 +
1155 TAN
1156 1/X
1157 59
1158 /
1159 +
1160 X<>Y
1161 RCL 11
1162 RTN
1163 LBL 13
1164 1
1165 XEQ 10
1166 RDN
1167 R-P
1168 X<> Z
1169 ST- Y
1170 X<> Z
1171 P-R
1172 R^
1173 XEQ 14
1174 RDN
1175 RTN
1176 LBL 10
1177 X<>Y
1178 RDN
1179 P-R
1180 R^
1181 X<>Y
1182 P-R
1183 RTN
1184 LBL 14
1185 R-P
1186 X<>Y
1187 RDN
1188 R-P
1189 R^
1190 X<>Y
1191 RTN
1192 LBL 12
1193 RCL 00
1194 2.06
1195 *
1196 STO 04
1197 4
1198 *
1199 STO 05
1200 52
1201 *
1202 STO 12
1203 19
1204 RCL 04
1205 *
1206 STO 11
1207 593
1208 LASTX
1209 *
1210 STO 07
1211 LASTX
1212 880
1213 *
1214 STO 08
1215 +
1216 STO 06
1217 LASTX
1218 RCL 07
1219 -
1220 STO 09
1221 +
1222 STO 10
1223 RCL 04
1224 31
1225 *
1226 STO 03
1227 RCL 04        
1228 102
1229 *
1230 END
 
 
    ( 1834 bytes / SIZE 016 )
 

    STACK           INPUTS                       OUTPUTS
         Z      YYYY.MNDD             h0 = apparent altitude  ( deg )
         Y      HH.MNSS (UT)                h  = true altitude  ( deg )
         X    DT = TT - UT ( sec)                 Az = Azimuth  ( deg )
 

   YYYY.MNDD   ENTER^
       HH.MNSS     ENTER^
        TT - UT        XEQ "EPH"   returns the azimuthal coordinates of the Sun.

-Then  R/S  or  XEQ 03   gives the coordinates of the Moon
-Then  R/S  or  XEQ 01   ------------------------  Mercury
     ............................................................................
-Then  R/S  or  XEQ 09   ------------------------- Pluto.


Example:         Seattle    Longitude = -122°333   Latitude = 47°600

   -122.333  STO 14
       47.6      STO 15

   Date      =  2024/05/01
   Time      =  12h41m  UT
 TT - UT  = 69.2 seconds


THE SUN
 

   2024.0501  ENTER^
      12.41       ENTER^
       69.2        XEQ "EPH"


               Az =   244°032   = R11
  RDN      h  =    -2°388    = R12
  RDN     h0 =    -1°693    = R13      ( often meaningless when < 0 )

and   R01 = X = 0.75310 AU
         R02 = Y = 0.66955 AU

  R03 = 41°639 = Geocentric longitude                 R06 = 39.205 = right-ascension
  R04 = 0  = Geocentric latitude                            R07 = 15.323  = declination
  R05 = 1.00770 AU = distance Earth-Sun


-And  R13 =  local true sidereal time = 8927°8538 = 287°8538 = 19h1903 = 19h11m25s


THE MOON

   R/S  or  XEQ 03     >>>>   Az =  -27.338
                                  RDN      h  =  15.330
                                  RDN     h0 =  15.390
                  

  R03 = 312.286 = Geocentric longitude                R06  = -43.843 = right-ascension
  R04 =  -4.680   = Geocentric latitude                   R07 =  -21.599  = declination
  R05 = 0.002496 AU = distance Earth-Moon


MERCURY


   R/S  or  XEQ 01   >>>>    Az =   267.526
                                RDN      h  =    3.973
                                RDN     h0 =    4.163
                              

  R03 =  17.426 = Geocentric longitude                  R06 = 16.993 = right-ascension           R08 = -100.006 = heliocentric longitude
  R04 = -2.421  = Geocentric latitude                     R07 =  4.602  = declination                 R09 = -3.661      = heliocentric latitude
  R05 = 0.70516 AU = distance Earth-Mercury                                                                R10 = 0.466560 AU = radius vector


VENUS


   R/S  or  XEQ 02    >>>>     Az =   252.961
                                 RDN      h  =   -0.246
                                 RDN     h0 =     0.273

  R03 = 32.539 = Geocentric longitude                 R06 =  30.770 = right-ascension          R08 = 19.828 = heliocentric longitude
  R04 = -1.214  = Geocentric latitude                   R07 = 11.212  = declination                 R09 =  -2.851  = heliocentric latitude
  R05 = 1.7019 AU = distance Earth-Venus                                                                     R10 =  0.7251 AU = radius vector


MARS


   R/S  or  XEQ 04   >>>>    Az =  -76.903
                                RDN      h  =    10.515
                                RDN     h0 =    10.601

  R03 =   0.692   = Geocentric longitude                 R06 = 1.138 = right-ascension            R08 = -27.876 = heliocentric longitude
  R04 =  -1.264   = Geocentric latitude                   R07 = -0.885 = declination                 R09 = -1.807   = heliocentric latitude
  R05 = 1.9745 AU = distance Earth-Mars                                                                        R10 = 1.3817 AU = radius vector


JUPITER


   R/S  or  XEQ 05     >>>>    Az = 232.437
                                  RDN      h  = -7.539
                                  RDN     h0 = -7.620   ( often meaningless when < 0 )

  R03 = 54.261 = Geocentric longitude                   R06 = 52.089  = right-ascension             R08 = 56.781 = heliocentric longitude
  R04 = -0.759   = Geocentric latitude                    R07 =  18.098  = declination                   R09 = -0.908  = heliocentric latitude
  R05 =  5.989 AU = distance Earth-Jupiter                                                                           R10 = 5.010 AU = radius vector


SATURN


   R/S  or  XEQ 06     >>>>   Az = -62.993
                                  RDN      h  =  13.915
                                  RDN     h0 =  13.981

  R03 = -13.314 = Geocentric longitude                   R06 = -11.564 = right-ascension            R08 = -18.194 = heliocentric longitude
  R04 = -1.753  = Geocentric latitude                      R07 =  -6.870  = declination                   R09 = -1.851  = heliocentric latitude
  R05 = 10.248 AU = distance Earth-Saturn                                                                          R10 = 9.704 AU = radius vector


URANUS


   R/S  or  XEQ 07   >>>>   Az =   234.000
                                RDN      h  =  -6.441
                                RDN     h0 =  -6.508   ( often meaningless when < 0 )

  R03 = 52.427 = Geocentric longitude                    R06 =  50.090 = right-ascension              R08 =  52.979 = heliocentric longitude
  R04 = -0.264  = Geocentric latitude                       R07 =  18.120 = declination                    R09 = -0.277  = heliocentric latitude
  R05 = 20.568 AU = distance Earth-Uranus                                                                           R10 = 19.579 AU = radius vector


NEPTUNE


   R/S  or  XEQ 08   >>>>   Az =   -75.225
                                 RDN      h  =  11.070
                                 RDN     h0 =  11.152

  R03 = -1.048  = Geocentric longitude                    R06 =  -0.472   = right-ascension              R08 = -2.358 = heliocentric longitude
  R04 = -1.231  = Geocentric latitude                       R07 =  -1.546   = declination                    R09 = -1.261  = heliocentric latitude
  R05 = 30.623 AU = distance Earth-Neptune                                                                          R10 = 29.891 AU = radius vector


PLUTO


   R/S  or  XEQ 09     >>>>   Az = -16.720
                                  RDN      h  =  18.069
                                  RDN     h0 =  18.119

  R03 = -57.904 = Geocentric longitude                    R06 =  -54.908 = right-ascension             R08 =  -59.533 = heliocentric longitude
  R04 = -3.027  = Geocentric latitude                       R07 =  -22.638  = declination                    R09 = -3.011  = heliocentric latitude
  R05 = 34.825 AU = distance Earth-Pluto                                                                                R10 = 35.006 AU = radius vector


Notes:

-The number of days since 2000/01/01 are computed with less bytes than the first "EPH"
-But it doesn't work if the date is before 1900/03/01 or after 2100/02/28

-However, if the date is for instance:  1900/01/06  put  1900/01/07  in Z-register
-And with  2100/03/03  put  2100/03/02

-HP41 in mode DEG
-Z output = h0 = apparent altitude but it is not in a register.


      d)  Program4  [ 1900/03/01 , 2100/02/28 ] ( 1517 bytes )


-Nutation, aberration and time-light are not computed
-Refraction is computed with Bennet formula.

-There are less decimals:  2 decimals with the angles ( deg ), 4 decimals with e ( eccentricity )

-I've found these maximum errors with heliocentric longitudes:

  Mars:     0°0152
  Jupiter:   0°0184
  Saturn:   0°0147
  Uranus:  0°0106
Neptune:  0°0087


Registers
 

 R00 = Time from 2000/01/01 0h TT ( unit = 100 Julian years )            ( Registers R14 & R15 are to be initialized before executing "EPH" )

 R01 = X   Rectangular
 R02 = Y   coordinates of the Sun,

 R03 = Geocentric longitude ( deg )
 R04 = Geocentric latitude ( deg )
 R05 = Distance Earth-Planet ( AU )

 R06 = Right-ascension ( deg )
 R07 = Declination ( deg )

 R08 = Heliocentric longitude ( deg )
 R09 = Heliocentric latitude ( deg )
 R10 = Radius vector ( AU )

 R11 = Azimuth ( deg )  ( reckoned clockwise from South )
 R12 = height ( deg )

 R13 = local sidereal time ( deg )

•  R14 = Longitude of the observer ( deg )  positive East
•  R15 = Latitude of the observer ( deg)

Flags: /
Subroutines:  /



 01 LBL "EPH"
 02 19
 03 10^X
 04 SQRT
 05 /
 06 STO 00
 07 RDN
 08 HR
 09 STO 13
 10 24
 11 /
 12 X<>Y
 13 INT
 14 STO 01
 15 X<> L
 16 FRC
 17 2
 18 10^X
 19 STO 09
 20 *
 21 INT
 22 LASTX
 23 FRC
 24 RCL 09
 25 *
 26 X<>Y
 27 3
 28 +
 29 6
 30 X>Y?
 31 DSE 01
 32 INT
 33 X<=Y?
 34 CHS
 35 +
 36 4
 37 +
 38 30.6
 39 *
 40 INT
 41 +
 42 RCL 01
 43 36525
 44 STO 02
 45 *
 46 RCL 09
 47 /
 48 INT
 49 +
 50 730563
 51 -
 52 +
 53 RCL 02
 54 /
 55 ST+ 00
 56 3600077
 57 STO 06
 58 *
 59 RCL 09
 60 ST/ Y
 61 +
 62 31
 63 1/X
 64 -
 65 RCL 13
 66 15
 67 *
 68 +
 69 RCL 14
 70 +
 71 STO 13
 72 167
 73 STO 08
 74 5
 75 +
 76 STO 07
 77 8003
 78 CHS
 79 STO 12
 80 297
 81 +
 82 STO 10
 83 CLST
 84 STO 01
 85 STO 02
 86 STO 11
 87 XEQ 13
 88 CLX
 89 RCL 03
 90 RCL 05
 91 P-R
 92 STO 01
 93 RDN
 94 STO 02
 95 CLX
 96 RCL 11
 97 RTN
 98 LBL 03
 99 RCL 00
100 3599905
101 %
102 3
103 -
104 STO 11
105 STO 12
106 RCL 00
107 44526711
108 %
109 68.25
110 -
111 STO 03
112 ST+ X
113 STO 05
114 STO 06
115 STO 10
116 ST- 11
117 RCL 00
118 47719887
119 %
120 51.57
121 -
122 ST- 05
123 STO 07
124 ST+ 10
125 ST+ X
126 STO 08
127 RCL 00
128 48320202
129 %
130 86.66
131 +
132 STO 09
133 RCL 06
134 +
135 SIN
136 PI
137 *
138 RCL 09
139 COS
140 RCL 07          
141 SIN
142 *
143 56
144 *
145 -
146 RCL 06
147 RCL 09
148 -
149 STO 04
150 SIN
151 17
152 *
153 +
154 RCL 05
155 SIN
156 RCL 09
157 COS
158 *
159 10
160 *
161 -
162 RCL 08
163 RCL 09
164 +
165 SIN
166 3
167 SQRT
168 *
169 +
170 RCL 10
171 RCL 09
172 -
173 SIN
174 -
175 RCL 08
176 RCL 09
177 -
178 SIN
179 +
180 RCL 09
181 RCL 11
182 +
183 SIN
184 -
185 RCL 09
186 SIN
187 513
188 *
189 +
190 RCL 04
191 RCL 08
192 -
193 SIN
194 RCL 09
195 RCL 10
196 +
197 SIN
198 -
199 2
200 /
201 +
202 STO 04
203 RCL 08
204 SIN
205 21
206 *
207 RCL 06
208 SIN
209 6
210 *
211 RCL 09
212 ST+ X
213 STO 09
214 SIN
215 -
216 11
217 *
218 +
219 RCL 07
220 SIN
221 629
222 *
223 -
224 RCL 06
225 RCL 08
226 -
227 SIN
228 RCL 05
229 RCL 12
230 -
231 SIN
232 -
233 6
234 *
235 +
236 RCL 10
237 SIN
238 RCL 11
239 SIN
240 +
241 5
242 *
243 -
244 RCL 12
245 RCL 07
246 -
247 SIN
248 4
249 *
250 +
251 RCL 05
252 SIN
253 127
254 *
255 -
256 RCL 07
257 COS
258 RCL 09
259 SIN
260 *
261 ST+ X
262 +
263 RCL 06
264 ST+ X
265 RCL 07
266 -
267 SIN
268 -
269 RCL 07
270 RCL 08
271 +
272 SIN
273 -
274 RCL 05
275 ST+ X
276 SIN
277 +
278 RCL 05
279 RCL 12
280 +
281 SIN
282 +
283 RCL 03
284 RCL 07
285 -
286 SIN
287 RCL 06
288 RCL 12
289 +
290 SIN
291 -
292 RCL 03          
293 RCL 12
294 +
295 SIN
296 +
297 RCL 10
298 RCL 12
299 -
300 SIN
301 -
302 RCL 12
303 SIN
304 37
305 *
306 -
307 RCL 03
308 SIN
309 7
310 *
311 -
312 RCL 07
313 RCL 12
314 +
315 SIN
316 6
317 *
318 +
319 RCL 06
320 RCL 09
321 -
322 SIN
323 3
324 *
325 +
326 2
327 /
328 +
329 48126788
330 RCL 00
331 *
332 21173
333 +
334 +
335 2
336 10^X
337 ST/ 04
338 /
339 360
340 MOD
341 257
342 RCL 07
343 COS
344 14
345 *
346 +
347 5
348 10^X
349 /
350 STO 05
351 XEQ 00
352 RTN
353 LBL 01
354 38.71
355 STO 09
356 2056
357 STO 08
358 700
359 STO 11
360 156
361 STO 07
362 7746
363 STO 10
364 25020
365 STO 12
366 14947407
367 STO 06
368 4833
369 119
370 XEQ 13
371 RTN
372 LBL 02
373 72.333
374 STO 09
375 68
376 STO 08
377 339
378 STO 11
379 140
380 STO 07
381 13156
382 STO 10
383 18118
384 STO 12
385 5851921
386 STO 06
387 7668
388 90
389 XEQ 13
390 RTN
391 LBL 04
392 XEQ 12
393 RCL 09
394 22
395 *
396 RCL 03
397 +
398 9
399 -
400 COS
401 CHS
402 RCL 06
403 47
404 +
405 COS
406 +
407 STO 03
408 CLX
409 STO 04
410 STO 05
411 152.37
412 STO 09
413 934
414 RCL 00
415 +
416 STO 08
417 185
418 STO 11
419 DSE X
420 STO 07
421 481.6
422 CHS
423 STO 12
424 1914170
425 STO 06
426 2394
427 CHS
428 STO 10
429 4956
430 77
431 XEQ 14
432 RTN
433 LBL 05
434 XEQ 12
435 RCL 08
436 73
437 +
438 COS
439 RCL 10          
440 RCL 07
441 -
442 41
443 +
444 SIN
445 -
446 5
447 SQRT
448 *
449 RCL 06
450 RCL 11
451 ST- 06
452 +
453 4
454 +
455 SIN
456 3
457 SQRT
458 *
459 -
460 RCL 09
461 76
462 -
463 SIN
464 RCL 06
465 60
466 -
467 SIN
468 +
469 4
470 *
471 +
472 RCL 10
473 33
474 -
475 SIN
476 6
477 *
478 +
479 RCL 11
480 7
481 -
482 SIN
483 33
484 *
485 -
486 STO 03
487 CLX
488 STO 04
489 STO 05
490 520.3
491 STO 09
492 4
493 /
494 STO 11
495 303630
496 STO 06
497 485
498 RCL 00
499 ST+ X
500 +
501 STO 08
502 161
503 STO 07
504 3431
505 STO 12
506 1433
507 STO 10
508 8613
509 +
510 102
511 XEQ 14
512 RTN
513 LBL 06
514 XEQ 12
515 RCL 06
516 RCL 09
517 RCL 11
518 +
519 STO 12
520 -
521 32
522 +
523 COS
524 RCL 11
525 ST+ X
526 76
527 +
528 COS
529 -
530 RCL 05
531 3
532 *
533 SIN
534 -
535 RCL 08
536 29
537 -
538 COS
539 +
540 RCL 10
541 53
542 +
543 COS
544 +
545 RCL 12
546 3
547 -
548 SIN
549 2
550 SQRT
551 *
552 +
553 RCL 11
554 40
555 -
556 COS
557 RCL 00
558 *
559 -
560 RCL 07
561 RCL 11
562 +
563 STO 04
564 49
565 -
566 SIN
567 5
568 *
569 +
570 RCL 09
571 14
572 +
573 COS
574 12
575 *
576 +
577 RCL 09
578 ST+ X
579 ST+ 07
580 30
581 +
582 STO 12
583 COS
584 23
585 *
586 +
587 RCL 11          
588 7
589 -
590 SIN
591 81
592 *
593 +
594 STO 03
595 RCL 12
596 3
597 +
598 COS
599 RCL 04
600 20
601 ST- 08
602 +
603 COS
604 -
605 ST+ X
606 STO 04
607 RCL 12
608 SIN
609 ST+ X
610 RCL 08
611 COS
612 +
613 STO 05
614 122351
615 STO 06
616 954
617 STO 09
618 555
619 RCL 00
620 PI
621 *
622 -
623 STO 08
624 249
625 STO 11
626 196
627 STO 07
628 5006
629 STO 12
630 9306
631 STO 10
632 2061
633 +
634 88
635 XEQ 14
636 RTN
637 LBL 07
638 XEQ 12
639 STO 07
640 ST+ X
641 STO 11
642 RCL 12
643 16
644 ST+ 07
645 -
646 STO 08
647 +
648 SIN
649 RCL 07
650 SIN
651 -
652 RCL 06
653 RCL 12
654 -
655 11
656 ST- 03
657 -
658 COS
659 RCL 12
660 RCL 04
661 -
662 4
663 +
664 COS
665 -
666 2
667 SQRT
668 *
669 +
670 RCL 05
671 RCL 08
672 +
673 COS
674 RCL 12
675 RCL 03
676 -
677 61
678 -
679 COS
680 +
681 4
682 *
683 +
684 RCL 03
685 COS
686 RCL 04
687 34
688 +
689 SIN
690 +
691 PI
692 *
693 -
694 RCL 11
695 19
696 +
697 STO 11
698 SIN
699 21
700 *
701 -
702 RCL 05
703 62
704 -
705 SIN
706 86
707 *
708 +
709 STO 03
710 29
711 RCL 12
712 -
713 SIN
714 STO 04
715 RCL 11
716 COS
717 PI
718 *
719 STO 05
720 42986
721 STO 06
722 1919
723 STO 09
724 4596
725 CHS
726 STO 12
727 463
728 STO 08
729 77
730 STO 11
731 149
732 STO 07
733 17299
734 STO 10
735 7400
736 52
737 XEQ 14
738 RTN
739 LBL 08
740 XEQ 12
741 STO 07
742 RCL 05          
743 -
744 COS
745 RCL 04
746 51
747 -
748 COS
749 +
750 5
751 SQRT
752 *
753 RCL 06
754 RCL 07
755 ST+ 04
756 RCL 05
757 +
758 -
759 COS
760 +
761 RCL 07
762 7
763 +
764 SIN
765 8
766 *
767 +
768 RCL 05
769 28
770 +
771 COS
772 58
773 *
774 +
775 STO 03
776 RCL 04
777 69
778 -
779 SIN
780 RCL 07
781 36
782 -
783 COS
784 -
785 STO 04
786 RCL 07
787 6
788 +
789 COS
790 3
791 SQRT
792 *
793 CHS
794 STO 05
795 21988
796 STO 06
797 3007
798 STO 09
799 90
800 STO 08
801 177
802 RCL 00
803 -
804 STO 11
805 143
806 STO 07
807 5565
808 CHS
809 STO 12
810 4812
811 STO 10
812 8366
813 +
814 110
815 XEQ 14
816 RTN
817 LBL 09
818 3949
819 STO 09
820 2490
821 STO 08
822 1714
823 STO 11
824 22407
825 STO 10
826 23893
827 STO 12
828 14657
829 STO 06
830 3626
831 -
832 140
833 STO 07
834 LBL 13
835 0
836 STO 03
837 STO 04
838 STO 05
839 RDN
840 LBL 14
841 RCL 00
842 ST* 06
843 ST* 07
844 *
845 +
846 X<> 12
847 9
848 X<> 06
849 +
850 RCL 07
851 RCL 10
852 +
853 4
854 10^X
855 ST/ 08
856 SQRT
857 ST/ 03
858 ST/ 04
859 ST/ 05
860 ST/ 09
861 ST/ 11
862 ST/ 12
863 ST/ Z
864 /
865 STO 10
866 -
867 ENTER
868 ENTER
869 LBL 11
870 SIN
871 RCL 08
872 R-D
873 *
874 +
875 DSE 06
876 GTO 11
877 STO 07
878 2
879 /
880 TAN
881 RCL 08
882 1
883 ST+ Y
884 RCL 08
885 -
886 /
887 SQRT
888 *
889 ATAN
890 ST+ X
891 1
892 RCL 07          
893 COS
894 RCL 08
895 *
896 -
897 RCL 09
898 *
899 RCL 05
900 +
901 X<> 10
902 +
903 RCL 12
904 ST+ 03
905 -
906 COS
907 RCL 11
908 LASTX
909 SIN
910 P-R
911 X<>Y
912 ASIN
913 RCL 04
914 +
915 STO 09
916 X<> Z
917 R-P
918 CLX
919 RCL 03
920 +
921 STO 08
922 RCL 10
923 XEQ 10
924 RCL 02
925 ST+ Z
926 CLX
927 RCL 01
928 +
929 XEQ 14
930 STO 05
931 X<> Z
932 STO 04
933 LBL 00
934 RCL 00
935 77
936 /
937 549.4
938 SQRT
939 -
940 RCL 04
941 R^
942 STO 03
943 XEQ 13
944 STO 06
945 RCL 15
946 RCL Z
947 STO 07
948 90
949 STO 11
950 ST- Z
951 R^
952 +
953 RCL 13
954 -
955 XEQ 13
956 ST- 11
957 X<>Y
958 ENTER
959 COS
960 RCL 05
961 409
962 *
963 /
964 -
965 STO 12
966 10.3
967 RCL 12
968 5.11
969 +
970 /
971 +
972 TAN
973 1/X
974 59
975 /
976 +
977 X<>Y
978 RCL 11
979 RTN
980 LBL 13
981 1
982 XEQ 10
983 RDN
984 R-P
985 X<> Z
986 ST- Y
987 X<> Z
988 P-R
989 R^
990 XEQ 14
991 RDN
992 RTN
993 LBL 10
994 X<>Y
995 RDN
996 P-R
997 R^
998 X<>Y
999 P-R
1000 RTN
1001 LBL 14
1002 R-P
1003 X<>Y
1004 RDN
1005 R-P
1006 R^
1007 X<>Y
1008 RTN
1009 LBL 12
1010 4
1011 RCL 00
1012 2.06
1013 *
1014 STO 04
1015 *
1016 STO 05
1017 52
1018 *
1019 STO 12
1020 19
1021 RCL 04
1022 *
1023 STO 11
1024 593
1025 LASTX
1026 *
1027 STO 07
1028 LASTX
1029 880
1030 *
1031 STO 08
1032 +
1033 STO 06
1034 LASTX
1035 RCL 07
1036 -
1037 STO 09        
1038 +
1039 STO 10
1040 31
1041 RCL 04
1042 *
1043 STO 03
1044 RCL 05
1045 ST+ X
1046 STO 04
1047 102
1048 LASTX
1049 *
1050 END
 
 
    ( 1517 bytes / SIZE 016 )
 

    STACK           INPUTS                       OUTPUTS
         Z      YYYY.MNDD             h0 = apparent altitude  ( deg )
         Y      HH.MNSS (UT)                h  = true altitude  ( deg )
         X    DT = TT - UT ( sec)                 Az = Azimuth  ( deg )
 

   YYYY.MNDD   ENTER^
       HH.MNSS     ENTER^
        TT - UT        XEQ "EPH"   returns the azimuthal coordinates of the Sun.

-Then  R/S  or  XEQ 03   gives the coordinates of the Moon
-Then  R/S  or  XEQ 01   ------------------------  Mercury
     ............................................................................
-Then  R/S  or  XEQ 09   ------------------------- Pluto.


Example:         Seattle    Longitude = -122°333   Latitude = 47°600

   -122.333  STO 14
       47.6      STO 15

   Date      =  2024/05/01
   Time      =  12h41m  UT
 TT - UT  = 69.2 seconds


THE SUN
 

   2024.0501  ENTER^
      12.41       ENTER^
       69.2        XEQ "EPH"


               Az =   244°035   = R11
  RDN      h  =    -2°387    = R12
  RDN     h0 =    -1°691    = R13      ( often meaningless when < 0 )

and   R01 = X = 0.75314 AU
         R02 = Y = 0.66950 AU

  R03 = 41°635 = Geocentric longitude                 R06 = 39.201 = right-ascension
  R04 = 0  = Geocentric latitude                            R07 = 15.322  = declination
  R05 = 1.00769 AU = distance Earth-Sun


-And  R13 =  local true sidereal time = 8927°8538 = 287°8538 = 19h1903 = 19h11m25s


THE MOON

   R/S  or  XEQ 03     >>>>   Az =  -27.347
                                  RDN      h  =  15.313
                                  RDN     h0 =  15.372
                  

  R03 = 312.295 = Geocentric longitude                R06  = -43.834 = right-ascension
  R04 =  -4.682   = Geocentric latitude                   R07 =  -21.599  = declination
  R05 = 0.002457 AU = distance Earth-Moon


MERCURY


   R/S  or  XEQ 01   >>>>    Az =   267.529
                                RDN      h  =    3.976
                                RDN     h0 =    4.166
                              

  R03 =  17.422 = Geocentric longitude                  R06 = 16.989 = right-ascension           R08 = -100.007 = heliocentric longitude
  R04 = -2.419  = Geocentric latitude                     R07 =  4.603  = declination                 R09 = -3.658      = heliocentric latitude
  R05 = 0.70517 AU = distance Earth-Mercury                                                                R10 = 0.46659  AU = radius vector


VENUS


   R/S  or  XEQ 02    >>>>     Az =   252.964
                                 RDN      h  =   -0.242
                                 RDN     h0 =     0.276

  R03 = 32.535 = Geocentric longitude                 R06 =  30.765 = right-ascension          R08 = 19.823 = heliocentric longitude
  R04 = -1.212  = Geocentric latitude                   R07 = 11.212  = declination                 R09 =  -2.847  = heliocentric latitude
  R05 = 1.7019 AU = distance Earth-Venus                                                                     R10 =  0.7252 AU = radius vector


MARS


   R/S  or  XEQ 04   >>>>    Az =  -76.901
                                RDN      h  =    10.516
                                RDN     h0 =    10.602

  R03 =   0.690   = Geocentric longitude                 R06 = 1.136 = right-ascension            R08 = -27.877 = heliocentric longitude
  R04 =  -1.264   = Geocentric latitude                   R07 = -0.886 = declination                 R09 = -1.807   = heliocentric latitude
  R05 = 1.9745 AU = distance Earth-Mars                                                                        R10 = 1.3817 AU = radius vector


JUPITER


   R/S  or  XEQ 05     >>>>    Az = 232.430
                                  RDN      h  = -7.537
                                  RDN     h0 = -7.619   ( often meaningless when < 0 )

  R03 = 54.267 = Geocentric longitude                   R06 = 52.094  = right-ascension             R08 = 56.789 = heliocentric longitude
  R04 = -0.755   = Geocentric latitude                    R07 =  18.103  = declination                   R09 = -0.903  = heliocentric latitude
  R05 =  5.988 AU = distance Earth-Jupiter                                                                           R10 = 5.010 AU = radius vector


SATURN


   R/S  or  XEQ 06     >>>>   Az = -62.993
                                  RDN      h  =  13.918
                                  RDN     h0 =  13.984
  R03 = -13.315 = Geocentric longitude                   R06 = -11.566 = right-ascension            R08 = -18.196 = heliocentric longitude
  R04 = -1.750  = Geocentric latitude                      R07 =  -6.868  = declination                   R09 = -1.849  = heliocentric latitude
  R05 = 10.245 AU = distance Earth-Saturn                                                                          R10 = 9.702 AU = radius vector


URANUS


   R/S  or  XEQ 07   >>>>   Az =   233.994
                                RDN      h  =  -6.442
                                RDN     h0 =  -6.509   ( often meaningless when < 0 )

  R03 = 52.432 = Geocentric longitude                    R06 =  50.095 = right-ascension              R08 =  52.985 = heliocentric longitude
  R04 = -0.262  = Geocentric latitude                       R07 =  18.123 = declination                   R09 = -0.275  = heliocentric latitude
  R05 = 20.567 AU = distance Earth-Uranus                                                                           R10 = 19.578 AU = radius vector


NEPTUNE


   R/S  or  XEQ 08   >>>>   Az =   -75.226
                                 RDN      h  =  11.069
                                 RDN     h0 =  11.150

  R03 = -1.045  = Geocentric longitude                    R06 =  -0.469 = right-ascension                R08 = -2.355 = heliocentric longitude
  R04 = -1.232  = Geocentric latitude                       R07 =  -1.546   = declination                    R09 = -1.262  = heliocentric latitude
  R05 = 30.623 AU = distance Earth-Neptune                                                                          R10 = 29.890 AU = radius vector


PLUTO


   R/S  or  XEQ 09     >>>>   Az = -16.726              
                                  RDN      h  =  18.068
                                  RDN     h0 =  18.119

  R03 = -57.899 = Geocentric longitude                    R06 =  -54.902 = right-ascension             R08 =  -59.528 = heliocentric longitude
  R04 = -3.027  = Geocentric latitude                       R07 =  -22.638  = declination                    R09 = -3.011  = heliocentric latitude
  R05 = 34.827 AU = distance Earth-Pluto                                                                                R10 = 35.008 AU = radius vector


Notes:

-The number of days since 2000/01/01 are computed with less bytes than the first "EPH"
-But it doesn't work if the date is before 1900/03/01 or after 2100/02/28

-However, if the date is for instance:  1900/01/06  put  1900/01/07  in Z-register
-And with  2100/03/03  put  2100/03/02

-HP41 in mode DEG
-Z output = h0 = apparent altitude but it is not in a register.



2°)  HP41CV  With X-Functions Module - Program  ( 867 bytes + 1516 or 1512 bytes )
 


Registers:
 

 R00 = Time from 2000/01/01 0h TT ( unit = 100 Julian years )            ( Registers R15-R16 are to be initialized before executing "EPH" )

 R01 = X   Rectangular
 R02 = Y   coordinates of the Sun,

 R03 = Geocentric longitude ( deg )
 R04 = Geocentric latitude ( deg )
 R05 = Distance Earth-Planet ( AU )

 R06 = Right-ascension ( deg )
 R07 = Declination ( deg )

 R08 = Heliocentric longitude ( deg )
 R09 = Heliocentric latitude ( deg )
 R10 = Radius vector ( AU )

 R11 = Azimuth ( deg )  ( reckoned clockwise from South )
 R12 = height ( deg )

 R13: temp

 R14 = local true sidereal time ( deg )
                            

•  R15 = Longitude of the observer ( deg )  positive East
•  R16 = Latitude of the observer ( deg)


-The first program computes the Sun, the Moon, Mercury, Venus and Mars.
-With this program ( if you call it "S" )   ,S  SAVEP  and   CLP



 01 -
 02 STO 03
 03 SIN
 04 RCL 00
 05 20
 06 *
 07 18
 08 -
 09 COS
 10 -
 11 RCL 00
 12 32964
 13 *
 14 22
 15 -
 16 COS
 17 -
 18 ST+ X
 19 X<> 03
 20 COS
 21 PI
 22 %
 23 STO 05
 24 CLST
 25 STO 01
 26 STO 02
 27 STO 04
 28 STO 11
 29 XEQ "X"
 30 CLX
 31 RCL 13
 32 RCL 05
 33 P-R
 34 STO 01
 35 RDN
 36 STO 02
 37 CLX
 38 RCL 11
 39 RTN
 40 LBL 03
 41 RCL 00
 42 3599905
 43 %
 44 2.96
 45 -
 46 STO 11
 47 STO 12
 48 RCL 00
 49 44526711
 50 %
 51 68.25
 52 -
 53 STO 13
 54 ST+ X
 55 STO 05
 56 STO 10
 57 ST- 11
 58 STO 06
 59 ST+ X
 60 STO 03
 61 RCL 00
 62 47719887
 63 %
 64 128.43
 65 +
 66 ST- 05
 67 STO 07
 68 ST+ 10
 69 ST+ X
 70 STO 08
 71 RCL 00
 72 48320202
 73 %
 74 86.66
 75 +
 76 STO 09
 77 RCL 07
 78 +
 79 SIN
 80 281
 81 *
 82 RCL 07          
 83 RCL 09
 84 -
 85 SIN
 86 278
 87 *
 88 +
 89 RCL 06
 90 RCL 09
 91 -
 92 SIN
 93 173
 94 *
 95 +
 96 RCL 05
 97 RCL 09
 98 +
 99 SIN
100 55
101 *
102 +
103 RCL 05
104 RCL 09
105 -
106 SIN
107 46
108 *
109 +
110 RCL 06
111 RCL 09
112 +
113 SIN
114 33
115 *
116 +
117 RCL 08
118 RCL 09
119 +
120 SIN
121 17
122 *
123 +
124 RCL 10
125 RCL 09
126 -
127 SIN
128 RCL 08
129 RCL 09
130 -
131 SIN
132 +
133 9
134 *
135 +
136 RCL 09
137 RCL 11
138 +
139 SIN
140 RCL 09
141 SIN
142 641
143 *
144 -
145 8
146 *
147 -
148 RCL 06
149 RCL 09
150 -
151 STO 04
152 RCL 08
153 -
154 SIN
155 RCL 09
156 RCL 10
157 +
158 SIN
159 +
160 4
161 *
162 +
163 RCL 04
164 RCL 12
165 +
166 SIN
167 PI
168 *
169 -
170 RCL 05          
171 RCL 12
172 -
173 STO 04
174 RCL 09
175 ST- 04
176 +
177 SIN
178 RCL 06
179 RCL 09
180 +
181 RCL 12
182 -
183 SIN
184 +
185 RCL 04
186 SIN
187 +
188 ST+ X
189 +
190 STO 04
191 RCL 07
192 SIN
193 6289
194 *
195 RCL 06
196 SIN
197 658
198 *
199 +
200 RCL 08
201 SIN
202 214
203 *
204 +
205 RCL 09
206 ST+ X
207 STO 09
208 SIN
209 114
210 *
211 -
212 RCL 06
213 RCL 08
214 -
215 SIN
216 59
217 *
218 +
219 RCL 05
220 RCL 12
221 -
222 SIN
223 57
224 *
225 +
226 RCL 10
227 SIN
228 53
229 *
230 +
231 RCL 11
232 SIN
233 46
234 *
235 -
236 RCL 12
237 RCL 07
238 -
239 SIN
240 41
241 *
242 -
243 RCL 05
244 SIN
245 98
246 *
247 RCL 07
248 RCL 09
249 +
250 SIN
251 -
252 13
253 *
254 +
255 RCL 07          
256 RCL 09
257 -
258 SIN
259 RCL 03
260 RCL 07
261 -
262 SIN
263 +
264 11
265 *
266 +
267 RCL 12
268 SIN
269 37
270 *
271 RCL 13
272 SIN
273 7
274 *
275 +
276 RCL 07
277 RCL 12
278 +
279 SIN
280 6
281 *
282 +
283 RCL 06
284 RCL 09
285 -
286 SIN
287 3
288 *
289 -
290 RCL 07
291 RCL 08
292 +
293 SIN
294 ST+ X
295 -
296 RCL 13
297 RCL 07
298 -
299 SIN
300 +
301 RCL 12
302 ST- 10
303 RCL 13
304 +
305 SIN
306 -
307 5
308 *
309 -
310 RCL 05
311 ST+ X
312 SIN
313 9
314 *
315 +
316 RCL 05
317 RCL 12
318 +
319 SIN
320 8
321 *
322 -
323 RCL 06
324 RCL 12
325 +
326 SIN
327 7
328 *
329 -
330 RCL 10
331 SIN
332 RCL 06
333 RCL 08
334 +
335 SIN
336 +
337 RCL 03
338 SIN
339 +
340 RCL 05          
341 RCL 08
342 -
343 SIN
344 +
345 RCL 00
346 132
347 *
348 60
349 -
350 SIN
351 -
352 4
353 *
354 +
355 RCL 08
356 RCL 12
357 ST+ 08
358 -
359 SIN
360 RCL 05
361 RCL 09
362 +
363 SIN
364 -
365 3
366 *
367 +
368 RCL 06
369 RCL 08
370 -
371 SIN
372 RCL 07
373 RCL 13
374 +
375 SIN
376 -
377 RCL 06
378 RCL 12
379 ST+ X
380 STO 09
381 -
382 SIN
383 +
384 RCL 08
385 SIN
386 -
387 RCL 09
388 SIN
389 -
390 RCL 05
391 RCL 09
392 -
393 SIN
394 +
395 ST+ X
396 +
397 481267881
398 RCL 00
399 ST+ X
400 -
401 RCL 00
402 *
403 211728
404 +
405 +
406 3
407 10^X
408 ST/ 04
409 /
410 360
411 MOD
412 STO 03
413 2574
414 RCL 07
415 COS
416 28
417 *
418 RCL 05
419 COS
420 5
421 *
422 +
423 RCL 06          
424 COS
425 4
426 *
427 +
428 5
429 *
430 -
431 6
432 10^X
433 /
434 STO 05
435 XEQ "Y"
436 RTN
437 LBL 01
438 387.1
439 STO 09
440 20563
441 STO 08
442 7005
443 RCL 00
444 2
445 STO 13
446 *
447 +
448 STO 11
449 1556
450 STO 07
451 77456
452 STO 10
453 250205
454 STO 12
455 149474072
456 STO 06
457 CLX
458 STO 03
459 STO 04
460 STO 05
461 48331
462 1186
463 XEQ "X"
464 RTN
465 LBL 02
466 RCL 00
467 45038
468 *
469 17
470 -
471 SIN
472 PI
473 *
474 STO 03
475 CLX
476 STO 04
477 STO 05
478 723.33
479 STO 09
480 677
481 STO 08
482 5
483 CHS
484 STO 13
485 3395
486 RCL 00
487 +
488 STO 11
489 1402
490 STO 07
491 131564
492 STO 10
493 181179
494 STO 12
495 58519213
496 STO 06
497 76680
498 901
499 XEQ "X"
500 RTN
501 LBL 04
502 XEQ "Z"
503 RCL 04
504 10
505 *
506 STO 05
507 SIN
508 RCL 00          
509 *
510 PI
511 *
512 RCL 03
513 264
514 *
515 13
516 +
517 COS
518 RCL 04
519 7077
520 *
521 SIN
522 -
523 RCL 04
524 532
525 *
526 49
527 +
528 COS
529 +
530 RCL 05
531 21
532 -
533 COS
534 7
535 *
536 +
537 ST+ X
538 +
539 RCL 09
540 22
541 *
542 RCL 03
543 +
544 ST+ 06
545 12
546 -
547 COS
548 6
549 *
550 -
551 RCL 06
552 49
553 +
554 COS
555 7
556 *
557 +
558 RCL 08
559 RCL 11
560 +
561 RCL 12
562 +
563 22
564 +
565 COS
566 RCL 06
567 ST+ X
568 12
569 +
570 COS
571 -
572 4
573 *
574 +
575 STO 03
576 CLX
577 STO 04
578 STO 05
579 1523.7
580 STO 09
581 9340
582 STO 08
583 1850
584 STO 11
585 9
586 STO 13
587 -
588 STO 07
589 4829
590 CHS
591 STO 12
592 19141696
593 STO 06          
594 23940
595 CHS
596 STO 10
597 49558
598 772
599 GTO "X"
600 END
 
 
    ( 867 bytes / SIZE 017 )



 01 LBL "EPH"
 02 19
 03 10^X
 04 SQRT
 05 /
 06 STO 00
 07 RDN
 08 HR
 09 STO 14
 10 24
 11 /
 12 X<>Y
 13 INT
 14 ST- L
 15 LASTX
 16 2
 17 10^X
 18 STO 02
 19 *
 20 ENTER
 21 INT
 22 STO 01
 23 -
 24 RCL 02
 25 *
 26 X<> 01
 27 1
 28 E^X
 29 -
 30 12
 31 /
 32 +
 33 INT
 34 ST- 01
 35 LASTX
 36 367
 37 *
 38 INT
 39 ST+ 01
 40 SIGN
 41 %
 42 INT
 43 .75
 44 CHS
 45 ST* Z
 46 *
 47 X<>Y
 48 RCL 01
 49 +
 50 INT
 51 +
 52 INT
 53 730434
 54 -
 55 +
 56 36525
 57 /
 58 ST+ 00
 59 36000770
 60 STO 06
 61 *
 62 3
 63 10^X
 64 STO 09
 65 /
 66 RCL 02
 67 +
 68 31
 69 1/X
 70 -
 71 RCL 14
 72 15
 73 *
 74 +
 75 RCL 15
 76 +
 77 STO 14
 78 1671
 79 STO 08
 80 48
 81 +
 82 STO 07
 83 4
 84 CHS
 85 STO 13
 86 80026
 87 CHS
 88 STO 12
 89 2963
 90 +
 91 STO 10
 92 CLA
 93 190
 94 XTOA
 95 RCL 00
 96 445267
 97 *
 98 68
 99 ASTO b
100 LBL 05
101 XEQ "Z"
102 RCL 08
103 73
104 +
105 COS
106 RCL 10
107 RCL 07
108 -
109 41
110 +
111 SIN
112 -
113 509
114 SQRT
115 *
116 RCL 08
117 3
118 *
119 28
120 +
121 COS
122 RCL 09
123 ST+ X
124 STO 12
125 25
126 +
127 COS
128 +
129 5
130 *
131 -
132 RCL 06
133 RCL 12
134 +
135 69
136 -
137 COS
138 RCL 06
139 RCL 11
140 +
141 STO 04
142 12
143 +
144 SIN
145 4
146 *
147 -
148 RCL 07          
149 COS
150 +
151 RCL 06
152 ST+ X
153 STO 12
154 RCL 09
155 +
156 79
157 +
158 COS
159 -
160 RCL 11
161 32
162 -
163 COS
164 RCL 00
165 *
166 +
167 4
168 *
169 +
170 RCL 11
171 ST- 06
172 ST- 12
173 ST+ X
174 75
175 ST- 09
176 +
177 COS
178 RCL 12
179 60
180 +
181 COS
182 -
183 RCL 13
184 3
185 ST* 05
186 *
187 76
188 +
189 COS
190 +
191 RCL 05
192 SIN
193 +
194 3
195 *
196 +
197 RCL 09
198 SIN
199 37
200 *
201 +
202 RCL 06
203 61
204 -
205 SIN
206 6
207 *
208 RCL 10
209 33
210 -
211 SIN
212 8
213 *
214 +
215 RCL 11
216 7.3
217 -
218 SIN
219 47
220 *
221 -
222 7
223 *
224 +
225 STO 03
226 RCL 04
227 18
228 +
229 COS
230 4
231 *
232 RCL 06
233 63
234 -
235 SIN
236 5
237 *
238 +
239 STO 04
240 RCL 10
241 COS
242 3
243 *
244 CHS
245 STO 05
246 5203
247 STO 09
248 1303
249 RCL 00
250 5
251 *
252 -
253 STO 11
254 3036303
255 STO 06
256 4849
257 STO 08
258 16
259 STO 13
260 1613
261 STO 07
262 34310
263 STO 12
264 14331
265 STO 10
266 86133
267 +
268 1021
269 XEQ "X"
270 RTN
271 LBL 06
272 XEQ "Z"
273 CLX
274 21
275 -
276 COS
277 RCL 11
278 ST+ X
279 STO 12
280 75
281 +
282 COS
283 -
284 RCL 05
285 R^
286 STO 05
287 +
288 SIN
289 -
290 8
291 *
292 RCL 06
293 14
294 +
295 COS
296 RCL 04          
297 7
298 *
299 62
300 +
301 COS
302 -
303 RCL 11
304 RCL 07
305 +
306 STO 04
307 LASTX
308 +
309 COS
310 -
311 RCL 07
312 RCL 12
313 -
314 36
315 +
316 COS
317 +
318 3
319 *
320 +
321 RCL 08
322 30
323 -
324 COS
325 RCL 10
326 54
327 +
328 COS
329 +
330 RCL 09
331 RCL 11
332 +
333 STO 10
334 RCL 06
335 -
336 25
337 -
338 COS
339 +
340 RCL 04
341 50
342 -
343 SIN
344 5
345 *
346 +
347 RCL 11
348 7.3
349 -
350 SIN
351 90
352 *
353 +
354 9
355 *
356 +
357 RCL 07
358 RCL 03
359 -
360 COS
361 RCL 08
362 RCL 13
363 -
364 RCL 05
365 -
366 79
367 -
368 COS
369 +
370 RCL 09
371 ST+ X
372 ST+ 07
373 30
374 +
375 STO 12
376 COS
377 38
378 *
379 -
380 6
381 *
382 -
383 RCL 10
384 SIN
385 7
386 *
387 RCL 11
388 32
389 -
390 COS
391 5
392 *
393 RCL 00
394 *
395 -
396 RCL 09
397 14
398 +
399 COS
400 59
401 *
402 +
403 ST+ X
404 +
405 STO 03
406 RCL 12
407 3
408 +
409 COS
410 20
411 ST+ 04
412 ST- 08
413 *
414 RCL 07
415 28
416 -
417 COS
418 RCL 04
419 COS
420 3
421 *
422 -
423 6
424 *
425 +
426 RCL 11
427 SIN
428 4
429 *
430 +
431 RCL 08
432 28
433 -
434 COS
435 RCL 13
436 3
437 *
438 76
439 -
440 COS
441 +
442 3
443 *
444 +
445 STO 04
446 RCL 12
447 SIN
448 19
449 *
450 RCL 08          
451 COS
452 8
453 *
454 +
455 RCL 09
456 73
457 -
458 COS
459 6
460 *
461 +
462 STO 05
463 1223511
464 STO 06
465 9543
466 STO 09
467 5553
468 STO 08
469 35
470 CHS
471 STO 13
472 2489
473 RCL 00
474 4
475 *
476 -
477 STO 11
478 1964
479 STO 07
480 50061
481 STO 12
482 93057
483 STO 10
484 20609
485 +
486 877
487 XEQ "X"
488 RTN
489 LBL 07
490 XEQ "Z"
491 5
492 *
493 CHS
494 R^
495 STO 04
496 RCL 05
497 +
498 STO 08
499 10
500 STO 10
501 +
502 SIN
503 7
504 *
505 +
506 RCL 12
507 ST- 07
508 RCL 13
509 ST+ X
510 STO 11
511 +
512 20
513 -
514 SIN
515 3
516 *
517 RCL 11
518 RCL 13
519 +
520 60
521 -
522 COS
523 +
524 4
525 *
526 +
527 RCL 07
528 28
529 +
530 COS
531 6
532 *
533 +
534 RCL 13
535 RCL 10
536 +
537 SIN
538 RCL 10
539 ST- 06
540 *
541 -
542 RCL 12
543 ST- 06
544 RCL 08
545 -
546 27
547 +
548 COS
549 RCL 12
550 RCL 04
551 -
552 COS
553 RCL 06
554 COS
555 -
556 5
557 *
558 +
559 3
560 *
561 -
562 RCL 12
563 RCL 03
564 -
565 51
566 -
567 COS
568 ST+ Y
569 RCL 05
570 RCL 12
571 +
572 STO 08
573 11
574 -
575 COS
576 +
577 RCL 05
578 62.2
579 -
580 SIN
581 22
582 *
583 +
584 39
585 ST+ 04
586 *
587 +
588 RCL 03
589 73
590 +
591 SIN
592 RCL 04
593 SIN
594 +
595 35
596 *
597 -
598 RCL 11          
599 19
600 +
601 STO 11
602 SIN
603 209
604 *
605 -
606 STO 03
607 RCL 11
608 14
609 +
610 COS
611 RCL 08
612 SIN
613 -
614 6
615 *
616 STO 04
617 RCL 11
618 COS
619 34
620 *
621 RCL 08
622 79
623 +
624 COS
625 7
626 *
627 -
628 STO 05
629 429864
630 STO 06
631 19192
632 STO 09
633 45955
634 CHS
635 STO 12
636 4630
637 STO 08
638 3
639 CHS
640 STO 13
641 773
642 RCL 00
643 +
644 STO 11
645 1486
646 STO 07
647 172993
648 STO 10
649 74005
650 521
651 XEQ "X"
652 RTN
653 LBL 08
654 XEQ "Z"
655 22
656 *
657 R^
658 STO 04
659 50
660 -
661 COS
662 24
663 *
664 +
665 RCL 06
666 RCL 13
667 ST+ 04
668 RCL 05
669 +
670 ST- 07
671 -
672 COS
673 9
674 *
675 +
676 RCL 04
677 77
678 -
679 COS
680 RCL 07
681 16
682 +
683 COS
684 +
685 5
686 *
687 +
688 RCL 13
689 7
690 +
691 SIN
692 71
693 *
694 +
695 RCL 05
696 27.8
697 +
698 COS
699 584
700 *
701 +
702 STO 03
703 RCL 04
704 69
705 -
706 SIN
707 RCL 13
708 36
709 -
710 COS
711 -
712 9
713 *
714 STO 04
715 RCL 05
716 63
717 -
718 COS
719 8
720 *
721 RCL 13
722 6
723 +
724 COS
725 17
726 *
727 -
728 STO 05
729 219883
730 STO 06
731 30070
732 STO 09
733 899
734 STO 08
735 CLX
736 STO 13
737 1770
738 RCL 00
739 9
740 *
741 -
742 STO 11
743 1426
744 STO 07
745 55654
746 CHS
747 STO 12
748 48124
749 STO 10
750 83660
751 +
752 1102
753 XEQ "X"
754 RTN
755 LBL 09
756 RCL 00          
757 2890
758 *
759 65
760 +
761 COS
762 7
763 *
764 RCL 00
765 1077
766 *
767 76
768 +
769 COS
770 4
771 *
772 +
773 STO 03
774 RCL 00
775 STO 05
776 ST+ 05
777 9
778 *
779 4
780 -
781 STO 04
782 39489
783 STO 09
784 24900
785 STO 08
786 CLX
787 STO 13
788 17140
789 STO 11
790 224075
791 STO 10
792 238927
793 STO 12
794 146567
795 STO 06
796 36260
797 -
798 1397
799 STO 07
800 LBL "X"
801 RCL 00
802 ST* 06
803 ST* 07
804 ST* 13
805 *
806 +
807 X<> 13
808 ST+ 08
809 5
810 10^X
811 ST/ 08
812 RCL 06
813 RCL 12
814 +
815 RCL 07
816 RCL 10
817 +
818 9
819 STO 06
820 SQRT
821 10^X
822 ST/ 03
823 ST/ 04
824 ST/ 05
825 ST/ 09
826 ST/ 11
827 ST/ 13
828 ST/ Z
829 /
830 STO 10
831 -
832 ENTER
833 ENTER
834 LBL 11
835 SIN
836 RCL 08
837 R-D
838 *
839 +
840 DSE 06
841 GTO 11
842 STO 07
843 2
844 /
845 TAN
846 RCL 08
847 1
848 ST+ Y
849 RCL 08
850 -
851 /
852 SQRT
853 *
854 ATAN
855 ST+ X
856 1
857 RCL 07
858 COS
859 RCL 08
860 *
861 -
862 RCL 09
863 *
864 RCL 05
865 +
866 X<> 10
867 +
868 RCL 13
869 ST+ 03
870 -
871 COS
872 RCL 11
873 LASTX
874 SIN
875 P-R
876 X<>Y
877 ASIN
878 RCL 04
879 +
880 STO 09
881 X<> Z
882 R-P
883 CLX
884 RCL 03
885 +
886 STO 08
887 RCL 10
888 XEQ 10
889 RCL 02
890 ST+ Z
891 CLX
892 RCL 01
893 +
894 XEQ 12
895 STO 05
896 RDN
897 STO 03
898 STO 13
899 X<>Y
900 STO 04
901 RCL 10
902 SQRT
903 1/X
904 RCL 10
905 ST/ Y
906 X^2
907 1
908 -
909 RCL 05          
910 X^2
911 ST+ Y
912 ST+ X
913 /
914 ST* Y
915 -
916 RCL 05
917 176
918 /
919 ST* Y
920 +
921 ST- 03
922 LBL "Y"
923 RCL 00
924 1934
925 *
926 55
927 +
928 1
929 P-R
930 391
931 /
932 549.4
933 SQRT
934 -
935 RCL 00
936 77
937 /
938 +
939 STO 12
940 COS
941 X<>Y
942 209
943 /
944 *
945 X<> 12
946 RCL 04
947 RCL 03
948 LASTX
949 -
950 STO 03
951 XEQ 13
952 STO 06
953 RCL 16
954 RCL Z
955 STO 07
956 90
957 STO 11
958 ST- Z
959 R^
960 +
961 RCL 12
962 +
963 RCL 14
964 -
965 XEQ 13
966 ST- 11
967 X<>Y
968 ENTER
969 COS
970 RCL 05
971 409
972 *
973 /
974 -
975 STO 12
976 10.3
977 RCL 12
978 5.11
979 +
980 /
981 +
982 TAN
983 1/X
984 59
985 /
986 +
987 X<>Y
988 RCL 11
989 RTN
990 LBL 13
991 1
992 XEQ 10
993 RDN
994 R-P
995 X<> Z
996 ST- Y
997 X<> Z
998 P-R
999 R^
1000 XEQ 12
1001 RDN
1002 RTN
1003 LBL 10
1004 X<>Y
1005 RDN
1006 P-R
1007 R^
1008 X<>Y
1009 P-R
1010 RTN
1011 LBL 12
1012 R-P
1013 X<>Y
1014 RDN
1015 R-P
1016 R^
1017 X<>Y
1018 RTN
1019 LBL "Z"
1020 102
1021 RCL 00
1022 2.06
1023 *
1024 STO 04
1025 *
1026 STO 13
1027 4
1028 LASTX
1029 *
1030 STO 05
1031 52
1032 *
1033 STO 12
1034 19
1035 RCL 04
1036 *
1037 STO 11
1038 593
1039 LASTX
1040 *
1041 STO 07
1042 LASTX
1043 880
1044 *
1045 STO 08
1046 +
1047 STO 06
1048 LASTX
1049 RCL 07
1050 -
1051 STO 09
1052 +
1053 STO 10
1054 RCL 05
1055 ST+ X
1056 RCL 04
1057 31
1058 *
1059 STO 03
1060 RCL 13        
1061 RCL 05
1062 -
1063 11
1064 -
1065 COS
1066 END
 
 
    ( 1516 bytes / SIZE 017 )


    STACK           INPUTS                       OUTPUTS
         Z      YYYY.MNDD             h0 = apparent altitude  ( deg )
         Y      HH.MNSS (UT)                h  = true altitude  ( deg )
         X    DT = TT - UT ( sec)                 Az = Azimuth  ( deg )
 

   YYYY.MNDD   ENTER^
       HH.MNSS     ENTER^
        TT - UT        XEQ "EPH"   returns the azimuthal coordinates of the Sun.

-Then  R/S  or  XEQ 03   gives the coordinates of the Moon
-Then  R/S  or  XEQ 01   ------------------------  Mercury

     ...........        XEQ 05  for Jupiter

-Then  R/S  or  XEQ 09   ------------------------- Pluto.


Example:         Seattle    Longitude = -122°333   Latitude = 47°600

   -122.333  STO 15
       47.6      STO 16

   Date      =  2024/05/01
   Time      =  12h41m  UT
 TT - UT  = 69.2 seconds


THE SUN
 

   2024.0501  ENTER^
      12.41       ENTER^
       69.2        XEQ "EPH"


               Az =   244°041   = R11                                                      
  RDN      h  =    -2°383    = R12
  RDN     h0 =    -1°687                  ( often meaningless when < 0 )

and   R01 = X = 0.75315 AU
         R02 = Y = 0.66949 AU

  R03 = 41°628 = Geocentric longitude                 R06 = 39.193 = right-ascension
  R04 = 0  = Geocentric latitude                            R07 = 15.321  = declination
  R05 = 1.00770 AU = distance Earth-Sun


-And  R14 =  local true sidereal time = 8927°8538 = 287°8538 = 19h1903 = 19h11m25s


THE MOON

  R/S  or  XEQ 03      >>>>   Az =  -27.344
                                  RDN      h  =  15.330
                                  RDN     h0 =  15.389
                  

  R03 = 312.291 = Geocentric longitude                R06  = -43.839 = right-ascension
  R04 =  -4.678   = Geocentric latitude                   R07 =  -21.598  = declination
  R05 = 0.002496 AU = distance Earth-Moon


MERCURY


   R/S  or  XEQ 01   >>>>    Az =   267.533                         
                                RDN      h  =    3.976
                                RDN     h0 =    4.165
                              

  R03 =  17.417 = Geocentric longitude                  R06 = 16.986 = right-ascension           R08 = -100.006 = heliocentric longitude
  R04 = -2.421  = Geocentric latitude                     R07 =  4.600  = declination                 R09 = -3.661      = heliocentric latitude
  R05 = 0.70519 AU = distance Earth-Mercury                                                                R10 = 0.46660  AU = radius vector


VENUS


   R/S  or  XEQ 02    >>>>     Az =   252.974                                        
                                 RDN      h  =   -0.240
                                 RDN     h0 =     0.278

  R03 = 32.524 = Geocentric longitude                 R06 =  30.755 = right-ascension          R08 = 19.830 = heliocentric longitude
  R04 = -1.214  = Geocentric latitude                   R07 = 11.208  = declination                 R09 =  -2.851  = heliocentric latitude
  R05 = 1.7019 AU = distance Earth-Venus                                                                     R10 =  0.7251 AU = radius vector


MARS


   R/S  or  XEQ 04   >>>>    Az =  -76.891
                                RDN      h  =    10.519
                                RDN     h0 =    10.604

  R03 =   0.679 = Geocentric longitude                   R06 = 1.126 = right-ascension            R08 = -27.879 = heliocentric longitude
  R04 =  -1.264   = Geocentric latitude                   R07 = -0.890  = declination                R09 = -1.807   = heliocentric latitude
  R05 = 1.9745 AU = distance Earth-Mars                                                                        R10 = 1.3817 AU = radius vector


JUPITER


  GTO "EPH" ( or CAT 1 ... )  and  XEQ 05     >>>>    Az = 232.446
                                                                         RDN      h  = -7.531
                                                                         RDN     h0 = -7.613   ( often meaningless when < 0 )

  R03 = 54.249 = Geocentric longitude                   R06 = 52.075  = right-ascension             R08 = 56.776 = heliocentric longitude
  R04 = -0.757   = Geocentric latitude                    R07 =  18.099  = declination                   R09 = -0.905  = heliocentric latitude
  R05 =  5.991 AU = distance Earth-Jupiter                                                                           R10 = 5.013 AU = radius vector


SATURN


   R/S  or  XEQ 06     >>>>   Az = -62.992
                                  RDN      h  =  13.918
                                  RDN     h0 =  13.984
  R03 = -13.317 = Geocentric longitude                   R06 = -11.568 = right-ascension            R08 = -18.189 = heliocentric longitude
  R04 = -1.749  = Geocentric latitude                      R07 =  -6.869  = declination                   R09 = -1.847  = heliocentric latitude
  R05 = 10.250 AU = distance Earth-Saturn                                                                          R10 = 9.706 AU = radius vector


URANUS


   R/S  or  XEQ 07   >>>>   Az =   234.008
                                RDN      h  =  -6.436
                                RDN     h0 =  -6.503   ( often meaningless when < 0 )

  R03 = 52.417 = Geocentric longitude                    R06 =  50.079 = right-ascension              R08 =  52.977 = heliocentric longitude
  R04 = -0.264  = Geocentric latitude                       R07 =  18.119  = declination                   R09 = -0.277  = heliocentric latitude
  R05 = 20.575 AU = distance Earth-Uranus                                                                           R10 = 19.586 AU = radius vector


NEPTUNE


   R/S  or  XEQ 08   >>>>   Az =   -75.223
                                 RDN      h  =  11.070
                                 RDN     h0 =  11.152

  R03 = -1.051  = Geocentric longitude                    R06 =  -0.474 = right-ascension                R08 = -2.354 = heliocentric longitude
  R04 = -1.231  = Geocentric latitude                       R07 =  -1.547   = declination                    R09 = -1.261  = heliocentric latitude
  R05 = 30.627 AU = distance Earth-Neptune                                                                          R10 = 29.895 AU = radius vector


PLUTO


   R/S  or  XEQ 09     >>>>   Az = -16.729
                                  RDN      h  =  18.065
                                  RDN     h0 =  18.115

  R03 = -57.897 = Geocentric longitude                    R06 =  -54.900 = right-ascension              R08 =  -59.525 = heliocentric longitude
  R04 = -3.028  = Geocentric latitude                       R07 =  -22.640  = declination                    R09 = -3.013  = heliocentric latitude
  R05 = 34.825 AU = distance Earth-Pluto                                                                                R10 = 35.006 AU = radius vector


Notes:

-HP41 in mode DEG

-Perhaps a better precision if we replace line 467 ( 5553 ) with 5554  and line 636 ( 4630 )  with  4634

-With geocentric coordinates, nutation, aberration and time-light are taken into account ( approximately )
-With this program, the apparent altitude is in Z-output but not in another register.

-We can add ( with Saturn ) after line 291:

 292 RCL 04
 293 771
 294 *
 295 78
 296 -
 297 COS
 298 RCL 04
 299 385
 300 *
 301 COS
 302 -
 303 RCL 04          
 304 177
 305 *
 306 51
 307 -
 308 COS
 309 -
 310 ST+ X
 311 +


-We can save 4 bytes if LBL 03 is in LBL "EPH":    LBL "Y"  becomes  LBL 00  and  XEQ "Y"  becomes  XEQ 00:


Registers:
 

 R00 = Time from 2000/01/01 0h TT ( unit = 100 Julian years )            ( Registers R15-R16 are to be initialized before executing "EPH" )

 R01 = X   Rectangular
 R02 = Y   coordinates of the Sun,

 R03 = Geocentric longitude ( deg )
 R04 = Geocentric latitude ( deg )
 R05 = Distance Earth-Planet ( AU )

 R06 = Right-ascension ( deg )
 R07 = Declination ( deg )

 R08 = Heliocentric longitude ( deg )
 R09 = Heliocentric latitude ( deg )
 R10 = Radius vector ( AU )

 R11 = Azimuth ( deg )  ( reckoned clockwise from South )
 R12 = height ( deg )

 R13: temp

 R14 = local true sidereal time ( deg )
                            

•  R15 = Longitude of the observer ( deg )  positive East
•  R16 = Latitude of the observer ( deg)


-The first program computes Mars, Jupiter, Saturn and Uranus.
-With this program ( if you call it "S" )   ,S  SAVEP  and   CLP



 01 COS
 02 +
 03 21
 04 RCL 05
 05 -
 06 COS
 07 7
 08 *
 09 +
 10 ST+ X
 11 +
 12 RCL 09
 13 22
 14 *
 15 RCL 03
 16 +
 17 ST+ 06
 18 12
 19 -
 20 COS
 21 6
 22 *
 23 -
 24 RCL 06
 25 49
 26 +
 27 COS
 28 7
 29 *
 30 +
 31 RCL 08
 32 RCL 11
 33 +
 34 RCL 12
 35 +
 36 22
 37 +
 38 COS
 39 RCL 06
 40 ST+ X
 41 12
 42 +
 43 COS
 44 -
 45 4
 46 *
 47 +
 48 STO 03
 49 CLX
 50 STO 04
 51 STO 05
 52 1523.7
 53 STO 09
 54 9340
 55 STO 08
 56 1850
 57 STO 11
 58 9
 59 STO 13
 60 -
 61 STO 07
 62 4829
 63 CHS
 64 STO 12
 65 19141696
 66 STO 06
 67 23940
 68 CHS
 69 STO 10
 70 49558
 71 772
 72 XEQ "X"
 73 RTN
 74 LBL 05
 75 XEQ "Z"
 76 RCL 08
 77 73
 78 +
 79 COS
 80 RCL 10
 81 RCL 07          
 82 -
 83 41
 84 +
 85 SIN
 86 -
 87 509
 88 SQRT
 89 *
 90 RCL 08
 91 3
 92 *
 93 28
 94 +
 95 COS
 96 RCL 09
 97 ST+ X
 98 STO 12
 99 25
100 +
101 COS
102 +
103 5
104 *
105 -
106 RCL 06
107 RCL 12
108 +
109 69
110 -
111 COS
112 RCL 06
113 RCL 11
114 +
115 STO 04
116 12
117 +
118 SIN
119 4
120 *
121 -
122 RCL 07
123 COS
124 +
125 RCL 06
126 ST+ X
127 STO 12
128 RCL 09
129 +
130 79
131 +
132 COS
133 -
134 RCL 11
135 32
136 -
137 COS
138 RCL 00
139 *
140 +
141 4
142 *
143 +
144 RCL 11
145 ST- 06
146 ST- 12
147 ST+ X
148 75
149 ST- 09
150 +
151 COS
152 RCL 12
153 60
154 +
155 COS
156 -
157 RCL 13
158 3
159 ST* 05
160 *
161 76
162 +
163 COS
164 +
165 RCL 05
166 SIN
167 +
168 3
169 *
170 +
171 RCL 09          
172 SIN
173 37
174 *
175 +
176 RCL 06
177 61
178 -
179 SIN
180 6
181 *
182 RCL 10
183 33
184 -
185 SIN
186 8
187 *
188 +
189 RCL 11
190 7.3
191 -
192 SIN
193 47
194 *
195 -
196 7
197 *
198 +
199 STO 03
200 RCL 04
201 18
202 +
203 COS
204 4
205 *
206 RCL 06
207 63
208 -
209 SIN
210 5
211 *
212 +
213 STO 04
214 RCL 10
215 COS
216 3
217 *
218 CHS
219 STO 05
220 5203
221 STO 09
222 1303
223 RCL 00
224 5
225 *
226 -
227 STO 11
228 3036303
229 STO 06
230 4849
231 STO 08
232 16
233 STO 13
234 1613
235 STO 07
236 34310
237 STO 12
238 14331
239 STO 10
240 86133
241 +
242 1021
243 XEQ "X"
244 RTN
245 LBL 06
246 XEQ "Z"
247 CLX
248 21
249 -
250 COS
251 RCL 11
252 ST+ X
253 STO 12
254 75
255 +
256 COS
257 -
258 RCL 05          
259 R^
260 STO 05
261 +
262 SIN
263 -
264 8
265 *
266 RCL 06
267 14
268 +
269 COS
270 RCL 04
271 7
272 *
273 62
274 +
275 COS
276 -
277 RCL 11
278 RCL 07
279 +
280 STO 04
281 LASTX
282 +
283 COS
284 -
285 RCL 07
286 RCL 12
287 -
288 36
289 +
290 COS
291 +
292 3
293 *
294 +
295 RCL 08
296 30
297 -
298 COS
299 RCL 10
300 54
301 +
302 COS
303 +
304 RCL 09
305 RCL 11
306 +
307 STO 10
308 RCL 06
309 -
310 25
311 -
312 COS
313 +
314 RCL 04
315 50
316 -
317 SIN
318 5
319 *
320 +
321 RCL 11
322 7.3
323 -
324 SIN
325 90
326 *
327 +
328 9
329 *
330 +
331 RCL 07
332 RCL 03
333 -
334 COS
335 RCL 08
336 RCL 13
337 -
338 RCL 05
339 -
340 79
341 -
342 COS
343 +
344 RCL 09          
345 ST+ X
346 ST+ 07
347 30
348 +
349 STO 12
350 COS
351 38
352 *
353 -
354 6
355 *
356 -
357 RCL 10
358 SIN
359 7
360 *
361 RCL 11
362 32
363 -
364 COS
365 5
366 *
367 RCL 00
368 *
369 -
370 RCL 09
371 14
372 +
373 COS
374 59
375 *
376 +
377 ST+ X
378 +
379 STO 03
380 RCL 12
381 3
382 +
383 COS
384 20
385 ST+ 04
386 ST- 08
387 *
388 RCL 07
389 28
390 -
391 COS
392 RCL 04
393 COS
394 3
395 *
396 -
397 6
398 *
399 +
400 RCL 11
401 SIN
402 4
403 *
404 +
405 RCL 08
406 28
407 -
408 COS
409 RCL 13
410 3
411 *
412 76
413 -
414 COS
415 +
416 3
417 *
418 +
419 STO 04
420 RCL 12
421 SIN
422 19
423 *
424 RCL 08
425 COS
426 8
427 *
428 +
429 RCL 09
430 73
431 -
432 COS
433 6
434 *
435 +
436 STO 05          
437 1223511
438 STO 06
439 9543
440 STO 09
441 5553
442 STO 08
443 35
444 CHS
445 STO 13
446 2489
447 RCL 00
448 4
449 *
450 -
451 STO 11
452 1964
453 STO 07
454 50061
455 STO 12
456 93057
457 STO 10
458 20609
459 +
460 877
461 XEQ "X"
462 RTN
463 LBL 07
464 XEQ "Z"
465 5
466 *
467 CHS
468 R^
469 STO 04
470 RCL 05
471 +
472 STO 08
473 10
474 STO 10
475 +
476 SIN
477 7
478 *
479 +
480 RCL 12
481 ST- 07
482 RCL 13
483 ST+ X
484 STO 11
485 +
486 20
487 -
488 SIN
489 3
490 *
491 RCL 11
492 RCL 13
493 +
494 60
495 -
496 COS
497 +
498 4
499 *
500 +
501 RCL 07
502 28
503 +
504 COS
505 6
506 *
507 +
508 RCL 13
509 RCL 10
510 +
511 SIN
512 RCL 10
513 ST- 06
514 *
515 -
516 RCL 12
517 ST- 06
518 RCL 08
519 -
520 27
521 +
522 COS
523 RCL 12
524 RCL 04
525 -
526 COS
527 RCL 06          
528 COS
529 -
530 5
531 *
532 +
533 3
534 *
535 -
536 RCL 12
537 RCL 03
538 -
539 51
540 -
541 COS
542 ST+ Y
543 RCL 05
544 RCL 12
545 +
546 STO 08
547 11
548 -
549 COS
550 +
551 RCL 05
552 62.2
553 -
554 SIN
555 22
556 *
557 +
558 39
559 ST+ 04
560 *
561 +
562 RCL 03
563 73
564 +
565 SIN
566 RCL 04
567 SIN
568 +
569 35
570 *
571 -
572 RCL 11
573 19
574 +
575 STO 11
576 SIN
577 209
578 *
579 -
580 STO 03
581 RCL 11
582 14
583 +
584 COS
585 RCL 08
586 SIN
587 -
588 6
589 *
590 STO 04
591 RCL 11
592 COS
593 34
594 *
595 RCL 08
596 79
597 +
598 COS
599 7
600 *
601 -
602 STO 05
603 429864
604 STO 06
605 19192
606 STO 09
607 45955
608 CHS
609 STO 12
610 4630
611 STO 08
612 3
613 CHS
614 STO 13
615 773
616 RCL 00          
617 +
618 STO 11
619 1486
620 STO 07
621 172993
622 STO 10
623 74005
624 521
625 GTO "X"
626 END
 
 
    ( 867 bytes / SIZE 017 )



 01 LBL "EPH"
 02 19
 03 10^X
 04 SQRT
 05 /
 06 STO 00
 07 RDN
 08 HR
 09 STO 14
 10 24
 11 /
 12 X<>Y
 13 INT
 14 ST- L
 15 LASTX
 16 2
 17 10^X
 18 STO 02
 19 *
 20 ENTER
 21 INT
 22 STO 01
 23 -
 24 RCL 02
 25 *
 26 X<> 01
 27 1
 28 E^X
 29 -
 30 12
 31 /
 32 +
 33 INT
 34 ST- 01
 35 LASTX
 36 367
 37 *
 38 INT
 39 ST+ 01
 40 SIGN
 41 %
 42 INT
 43 .75
 44 CHS
 45 ST* Z
 46 *
 47 X<>Y
 48 RCL 01
 49 +
 50 INT
 51 +
 52 INT
 53 730434
 54 -
 55 +
 56 36525
 57 /
 58 ST+ 00
 59 36000770
 60 STO 06
 61 *
 62 3
 63 10^X
 64 STO 09
 65 /
 66 RCL 02
 67 +
 68 31
 69 1/X
 70 -
 71 RCL 14
 72 15
 73 *
 74 +
 75 RCL 15
 76 +
 77 STO 14
 78 1671
 79 STO 08
 80 48
 81 +
 82 STO 07
 83 4
 84 CHS
 85 STO 13
 86 80026
 87 CHS
 88 STO 12
 89 2963
 90 +
 91 STO 10
 92 RCL 00
 93 445267
 94 *
 95 68
 96 -
 97 STO 03
 98 SIN
 99 RCL 00
100 20
101 *
102 18
103 -
104 COS
105 -
106 RCL 00
107 32964
108 *
109 22
110 -
111 COS
112 -
113 ST+ X
114 X<> 03
115 COS
116 PI
117 %
118 STO 05
119 CLST
120 STO 01
121 STO 02
122 STO 04
123 STO 11
124 XEQ "X"
125 CLX
126 RCL 13
127 RCL 05
128 P-R
129 STO 01
130 RDN
131 STO 02
132 CLX
133 RCL 11
134 RTN
135 LBL 03
136 RCL 00
137 3599905
138 %
139 2.96
140 -
141 STO 11
142 STO 12
143 RCL 00          
144 44526711
145 %
146 68.25
147 -
148 STO 13
149 ST+ X
150 STO 05
151 STO 10
152 ST- 11
153 STO 06
154 ST+ X
155 STO 03
156 RCL 00
157 47719887
158 %
159 128.43
160 +
161 ST- 05
162 STO 07
163 ST+ 10
164 ST+ X
165 STO 08
166 RCL 00
167 48320202
168 %
169 86.66
170 +
171 STO 09
172 RCL 07
173 +
174 SIN
175 281
176 *
177 RCL 07
178 RCL 09
179 -
180 SIN
181 278
182 *
183 +
184 RCL 06
185 RCL 09
186 -
187 SIN
188 173
189 *
190 +
191 RCL 05
192 RCL 09
193 +
194 SIN
195 55
196 *
197 +
198 RCL 05
199 RCL 09
200 -
201 SIN
202 46
203 *
204 +
205 RCL 06
206 RCL 09
207 +
208 SIN
209 33
210 *
211 +
212 RCL 08
213 RCL 09
214 +
215 SIN
216 17
217 *
218 +
219 RCL 10
220 RCL 09
221 -
222 SIN
223 RCL 08
224 RCL 09
225 -
226 SIN
227 +
228 9
229 *
230 +
231 RCL 09
232 RCL 11
233 +
234 SIN
235 RCL 09
236 SIN
237 641
238 *
239 -
240 8
241 *
242 -
243 RCL 06
244 RCL 09
245 -
246 STO 04
247 RCL 08
248 -
249 SIN
250 RCL 09
251 RCL 10
252 +
253 SIN
254 +
255 4
256 *
257 +
258 RCL 04
259 RCL 12
260 +
261 SIN
262 PI
263 *
264 -
265 RCL 05
266 RCL 12
267 -
268 STO 04
269 RCL 09
270 ST- 04
271 +
272 SIN
273 RCL 06
274 RCL 09
275 +
276 RCL 12
277 -
278 SIN
279 +
280 RCL 04
281 SIN
282 +
283 ST+ X
284 +
285 STO 04
286 RCL 07
287 SIN
288 6289
289 *
290 RCL 06          
291 SIN
292 658
293 *
294 +
295 RCL 08
296 SIN
297 214
298 *
299 +
300 RCL 09
301 ST+ X
302 STO 09
303 SIN
304 114
305 *
306 -
307 RCL 06
308 RCL 08
309 -
310 SIN
311 59
312 *
313 +
314 RCL 05
315 RCL 12
316 -
317 SIN
318 57
319 *
320 +
321 RCL 10
322 SIN
323 53
324 *
325 +
326 RCL 11
327 SIN
328 46
329 *
330 -
331 RCL 12
332 RCL 07
333 -
334 SIN
335 41
336 *
337 -
338 RCL 05
339 SIN
340 98
341 *
342 RCL 07
343 RCL 09
344 +
345 SIN
346 -
347 13
348 *
349 +
350 RCL 07
351 RCL 09
352 -
353 SIN
354 RCL 03
355 RCL 07
356 -
357 SIN
358 +
359 11
360 *
361 +
362 RCL 12
363 SIN
364 37
365 *
366 RCL 13
367 SIN
368 7
369 *
370 +
371 RCL 07
372 RCL 12
373 +
374 SIN
375 6
376 *
377 +
378 RCL 06
379 RCL 09
380 -
381 SIN
382 3
383 *
384 -
385 RCL 07
386 RCL 08
387 +
388 SIN
389 ST+ X
390 -
391 RCL 13
392 RCL 07
393 -
394 SIN
395 +
396 RCL 12
397 ST- 10
398 RCL 13
399 +
400 SIN
401 -
402 5
403 *
404 -
405 RCL 05
406 ST+ X
407 SIN
408 9
409 *
410 +
411 RCL 05
412 RCL 12
413 +
414 SIN
415 8
416 *
417 -
418 RCL 06
419 RCL 12
420 +
421 SIN
422 7
423 *
424 -
425 RCL 10
426 SIN
427 RCL 06
428 RCL 08
429 +
430 SIN
431 +
432 RCL 03
433 SIN
434 +
435 RCL 05
436 RCL 08
437 -
438 SIN
439 +
440 RCL 00          
441 132
442 *
443 60
444 -
445 SIN
446 -
447 4
448 *
449 +
450 RCL 08
451 RCL 12
452 ST+ 08
453 -
454 SIN
455 RCL 05
456 RCL 09
457 +
458 SIN
459 -
460 3
461 *
462 +
463 RCL 06
464 RCL 08
465 -
466 SIN
467 RCL 07
468 RCL 13
469 +
470 SIN
471 -
472 RCL 06
473 RCL 12
474 ST+ X
475 STO 09
476 -
477 SIN
478 +
479 RCL 08
480 SIN
481 -
482 RCL 09
483 SIN
484 -
485 RCL 05
486 RCL 09
487 -
488 SIN
489 +
490 ST+ X
491 +
492 481267881
493 RCL 00
494 ST+ X
495 -
496 RCL 00
497 *
498 211728
499 +
500 +
501 3
502 10^X
503 ST/ 04
504 /
505 360
506 MOD
507 STO 03
508 2574
509 RCL 07
510 COS
511 28
512 *
513 RCL 05
514 COS
515 5
516 *
517 +
518 RCL 06
519 COS
520 4
521 *
522 +
523 5
524 *
525 -
526 6
527 10^X
528 /
529 STO 05
530 XEQ 00
531 RTN
532 LBL 01
533 387.1
534 STO 09
535 20563
536 STO 08
537 7005
538 RCL 00
539 2
540 STO 13
541 *
542 +
543 STO 11
544 1556
545 STO 07
546 77456
547 STO 10
548 250205
549 STO 12
550 149474072
551 STO 06
552 CLX
553 STO 03
554 STO 04
555 STO 05
556 48331
557 1186
558 XEQ "X"
559 RTN
560 LBL 02
561 RCL 00
562 45038
563 *
564 17
565 -
566 SIN
567 PI
568 *
569 STO 03
570 CLX
571 STO 04
572 STO 05
573 723.33
574 STO 09
575 677
576 STO 08
577 5
578 CHS
579 STO 13
580 3395
581 RCL 00          
582 +
583 STO 11
584 1402
585 STO 07
586 131564
587 STO 10
588 181179
589 STO 12
590 58519213
591 STO 06
592 76680
593 901
594 XEQ "X"
595 RTN
596 LBL 04
597 CLA
598 190
599 XTOA
600 XEQ "Z"
601 RCL 04
602 10
603 *
604 STO 05
605 SIN
606 RCL 00
607 *
608 PI
609 *
610 RCL 03
611 264
612 *
613 13
614 +
615 COS
616 RCL 04
617 7077
618 *
619 SIN
620 -
621 RCL 04
622 532
623 *
624 49
625 +
626 ASTO b
627 LBL 08
628 XEQ "Z"
629 22
630 *
631 R^
632 STO 04
633 50
634 -
635 COS
636 24
637 *
638 +
639 RCL 06
640 RCL 13
641 ST+ 04
642 RCL 05
643 +
644 ST- 07
645 -
646 COS
647 9
648 *
649 +
650 RCL 04
651 77
652 -
653 COS
654 RCL 07
655 16
656 +
657 COS
658 +
659 5
660 *
661 +
662 RCL 13
663 7
664 +
665 SIN
666 71
667 *
668 +
669 RCL 05
670 27.8
671 +
672 COS
673 584
674 *
675 +
676 STO 03
677 RCL 04
678 69
679 -
680 SIN
681 RCL 13
682 36
683 -
684 COS
685 -
686 9
687 *
688 STO 04
689 RCL 05
690 63
691 -
692 COS
693 8
694 *
695 RCL 13
696 6
697 +
698 COS
699 17
700 *
701 -
702 STO 05
703 219883
704 STO 06
705 30070
706 STO 09
707 899
708 STO 08
709 CLX
710 STO 13
711 1770
712 RCL 00
713 9
714 *
715 -
716 STO 11
717 1426
718 STO 07
719 55654
720 CHS
721 STO 12
722 48124
723 STO 10
724 83660
725 +
726 1102
727 XEQ "X"
728 RTN
729 LBL 09
730 RCL 00          
731 2890
732 *
733 65
734 +
735 COS
736 7
737 *
738 RCL 00
739 1077
740 *
741 76
742 +
743 COS
744 4
745 *
746 +
747 STO 03
748 RCL 00
749 STO 05
750 ST+ 05
751 9
752 *
753 4
754 -
755 STO 04
756 39489
757 STO 09
758 24900
759 STO 08
760 CLX
761 STO 13
762 17140
763 STO 11
764 224075
765 STO 10
766 238927
767 STO 12
768 146567
769 STO 06
770 36260
771 -
772 1397
773 STO 07
774 LBL "X"
775 RCL 00
776 ST* 06
777 ST* 07
778 ST* 13
779 *
780 +
781 X<> 13
782 ST+ 08
783 5
784 10^X
785 ST/ 08
786 RCL 06
787 RCL 12
788 +
789 RCL 07
790 RCL 10
791 +
792 9
793 STO 06
794 SQRT
795 10^X
796 ST/ 03
797 ST/ 04
798 ST/ 05
799 ST/ 09
800 ST/ 11
801 ST/ 13
802 ST/ Z
803 /
804 STO 10
805 -
806 ENTER
807 ENTER
808 LBL 11
809 SIN
810 RCL 08
811 R-D
812 *
813 +
814 DSE 06
815 GTO 11
816 STO 07
817 2
818 /
819 TAN
820 RCL 08
821 1
822 ST+ Y
823 RCL 08
824 -
825 /
826 SQRT
827 *
828 ATAN
829 ST+ X
830 1
831 RCL 07
832 COS
833 RCL 08
834 *
835 -
836 RCL 09
837 *
838 RCL 05
839 +
840 X<> 10
841 +
842 RCL 13
843 ST+ 03
844 -
845 COS
846 RCL 11
847 LASTX
848 SIN
849 P-R
850 X<>Y
851 ASIN
852 RCL 04
853 +
854 STO 09
855 X<> Z
856 R-P
857 CLX
858 RCL 03
859 +
860 STO 08
861 RCL 10
862 XEQ 10
863 RCL 02
864 ST+ Z
865 CLX
866 RCL 01
867 +
868 XEQ 12
869 STO 05
870 RDN
871 STO 03
872 STO 13
873 X<>Y
874 STO 04
875 RCL 10
876 SQRT
877 1/X
878 RCL 10          
879 ST/ Y
880 X^2
881 1
882 -
883 RCL 05
884 X^2
885 ST+ Y
886 ST+ X
887 /
888 ST* Y
889 -
890 RCL 05
891 176
892 /
893 ST* Y
894 +
895 ST- 03
896 LBL 00
897 RCL 00
898 1934
899 *
900 55
901 +
902 1
903 P-R
904 391
905 /
906 549.4
907 SQRT
908 -
909 RCL 00
910 77
911 /
912 +
913 STO 12
914 COS
915 X<>Y
916 209
917 /
918 *
919 X<> 12
920 RCL 04
921 RCL 03
922 LASTX
923 -
924 STO 03
925 XEQ 13
926 STO 06
927 RCL 16
928 RCL Z
929 STO 07
930 90
931 STO 11
932 ST- Z
933 R^
934 +
935 RCL 12
936 +
937 RCL 14
938 -
939 XEQ 13
940 ST- 11
941 X<>Y
942 ENTER
943 COS
944 RCL 05
945 409
946 *
947 /
948 -
949 STO 12
950 10.3
951 RCL 12
952 5.11
953 +
954 /
955 +
956 TAN
957 1/X
958 59
959 /
960 +
961 X<>Y
962 RCL 11
963 RTN
964 LBL 13
965 1
966 XEQ 10
967 RDN
968 R-P
969 X<> Z
970 ST- Y
971 X<> Z
972 P-R
973 R^
974 XEQ 12
975 RDN
976 RTN
977 LBL 10
978 X<>Y
979 RDN
980 P-R
981 R^
982 X<>Y
983 P-R
984 RTN
985 LBL 12
986 R-P
987 X<>Y
988 RDN
989 R-P
990 R^
991 X<>Y
992 RTN
993 LBL "Z"
994 102
995 RCL 00
996 2.06
997 *
998 STO 04
999 *
1000 STO 13
1001 4
1002 LASTX
1003 *
1004 STO 05
1005 52
1006 *
1007 STO 12
1008 19
1009 RCL 04
1010 *
1011 STO 11
1012 593
1013 LASTX
1014 *
1015 STO 07
1016 LASTX
1017 880
1018 *
1019 STO 08
1020 +
1021 STO 06
1022 LASTX
1023 RCL 07
1024 -
1025 STO 09
1026 +
1027 STO 10
1028 RCL 05
1029 ST+ X
1030 RCL 04
1031 31
1032 *
1033 STO 03
1034 RCL 13        
1035 RCL 05
1036 -
1037 11
1038 -
1039 COS
1040 END
 
 
    ( 1512 bytes / SIZE 017 )


    STACK           INPUTS                       OUTPUTS
         Z      YYYY.MNDD             h0 = apparent altitude  ( deg )
         Y      HH.MNSS (UT)                h  = true altitude  ( deg )
         X    DT = TT - UT ( sec)                 Az = Azimuth  ( deg )
 

   YYYY.MNDD   ENTER^
       HH.MNSS     ENTER^
        TT - UT        XEQ "EPH"   returns the azimuthal coordinates of the Sun.

-Then  R/S  or  XEQ 03   gives the coordinates of the Moon
-Then  R/S  or  XEQ 01   ------------------------  Mercury

     ...........        XEQ 08  for Neptune

-Then  R/S  or  XEQ 09   ------------------------- Pluto.

References:

[1]  Jean Meeus - "Astronomical Algorithms" - Willmann-Bell  -  ISBN 0-943396-61-1
[2]  VSOP87D Series  ftp://ftp.imcce.fr/pub/ephem/planets/vsop87/
[3]  VSOP09 & TOP2010 Series  http://www.imcce.fr/~jlsimon
[4]  A fantastic software by Aldo Vitagliano: "SOLEX" which may be downloaded freely from http://chemistry.unina.it/~alvitagl/solex/