hp41programs

Astro5 Astronomical Ephemeris (VIII) for the HP41CV
 

Overview
 

1°)  HP41CV  without Module

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

 b)  Program2 [ 1900/03/01 , 2100/02/28 ]  ( 1945 bytes )
 c)
 Program3 [ 1900/03/01 , 2100/02/28 ]  ( 1821 bytes )
 d)  Program4 [ 1900/03/01 , 2100/02/28 ]  ( 1496 bytes )
 e)  
Program5 (?) [ 1900/03/01 , 2100/02/28 ]  ( 1083 bytes )

2°)  HP41CV  with X-Functions Module
- Program ( 867 bytes + 1503 or 1499 bytes )



1°)  HP41CV  Without Module
 

      a)  Program1  [ 1900 , 2100 ]  ( 2117 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 ST+ 12
529 *
530 RCL 03
531 +
532 ST+ 06
533 12
534 -
535 COS
536 6
537 *
538 -
539 RCL 06
540 49
541 +
542 COS
543 7
544 *
545 +
546 RCL 08
547 RCL 11
548 +
549 RCL 12
550 +
551 STO 10
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 ST+ 10
573 SIN
574 RCL 10
575 6
576 -
577 COS
578 -
579 ST+ X
580 -
581 STO 03
582 CLX
583 STO 04
584 STO 05
585 1523.7
586 STO 09
587 9340
588 STO 08
589 1850
590 STO 11
591 9
592 STO 13
593 -
594 STO 07
595 4829
596 CHS
597 STO 12
598 19141696
599 STO 06
600 23940
601 CHS
602 STO 10
603 49558
604 772
605 XEQ 14
606 RTN
607 LBL 05
608 XEQ 12
609 RCL 08
610 73
611 +
612 COS
613 RCL 10
614 RCL 07          
615 -
616 41
617 +
618 SIN
619 -
620 509
621 SQRT
622 *
623 RCL 08
624 3
625 *
626 28
627 +
628 COS
629 RCL 09
630 ST+ X
631 STO 12
632 25
633 +
634 COS
635 +
636 5
637 *
638 -
639 RCL 06
640 RCL 12
641 +
642 69
643 -
644 COS
645 RCL 06
646 RCL 11
647 +
648 STO 04
649 12
650 +
651 SIN
652 4
653 *
654 -
655 RCL 07
656 COS
657 +
658 RCL 06
659 ST+ X
660 RCL 09
661 +
662 79
663 +
664 COS
665 -
666 RCL 11
667 ST- 06
668 32
669 -
670 COS
671 RCL 00
672 *
673 +
674 4
675 *
676 +
677 RCL 09
678 75
679 -
680 SIN
681 37
682 *
683 +
684 RCL 06
685 61
686 -
687 SIN
688 6
689 *
690 RCL 10
691 33
692 -
693 SIN
694 8
695 *
696 +
697 RCL 11
698 7.3
699 -
700 SIN
701 47
702 *
703 -
704 7
705 *
706 +
707 STO 03
708 RCL 06
709 63
710 -
711 SIN
712 5
713 *
714 STO 04
715 CLX
716 STO 05
717 5203
718 STO 09
719 1303
720 RCL 00
721 5
722 *
723 -
724 STO 11
725 3036303
726 STO 06
727 4849
728 STO 08
729 16
730 STO 13
731 1613
732 STO 07
733 34310
734 STO 12
735 14331
736 STO 10
737 86133
738 +
739 1021
740 XEQ 14
741 RTN
742 LBL 06
743 XEQ 12
744 CLX
745 21
746 -
747 COS
748 RCL 11
749 ST+ X
750 75
751 +
752 COS
753 -
754 RCL 05
755 R^
756 STO 05
757 +
758 SIN
759 -
760 8
761 *
762 RCL 06
763 14
764 +
765 COS
766 PI
767 *
768 +
769 RCL 08
770 30
771 -
772 COS
773 RCL 10
774 54
775 +
776 COS
777 +
778 RCL 09
779 RCL 11
780 +
781 STO 10
782 RCL 06
783 -
784 25
785 -
786 COS
787 +
788 RCL 07
789 RCL 11
790 +
791 STO 04
792 50
793 -
794 SIN
795 5
796 *
797 +
798 RCL 11
799 7.3
800 -
801 SIN
802 90
803 *
804 +
805 9
806 *
807 +
808 RCL 07
809 RCL 03
810 -
811 COS
812 RCL 08
813 RCL 13
814 -
815 RCL 05
816 -
817 79
818 -
819 COS
820 +
821 RCL 09          
822 ST+ X
823 ST+ 07
824 30
825 +
826 STO 12
827 COS
828 38
829 *
830 -
831 6
832 *
833 -
834 RCL 10
835 SIN
836 7
837 *
838 RCL 11
839 32
840 -
841 COS
842 5
843 *
844 RCL 00
845 *
846 -
847 RCL 09
848 14
849 +
850 COS
851 59
852 *
853 +
854 ST+ X
855 +
856 STO 03
857 RCL 12
858 3
859 +
860 COS
861 20
862 ST+ 04
863 ST- 08
864 *
865 RCL 07
866 28
867 -
868 COS
869 RCL 04
870 COS
871 3
872 *
873 -
874 6
875 *
876 +
877 RCL 11
878 SIN
879 4
880 *
881 +
882 STO 04
883 RCL 12
884 SIN
885 19
886 *
887 RCL 08
888 COS
889 8
890 *
891 +
892 STO 05
893 1223511
894 STO 06
895 9543
896 STO 09
897 5553
898 STO 08
899 35
900 CHS
901 STO 13
902 2489
903 RCL 00
904 4
905 *
906 -
907 STO 11
908 1964
909 STO 07
910 50061
911 STO 12
912 93057
913 STO 10
914 20609
915 +
916 877
917 XEQ 14
918 RTN
919 LBL 07
920 XEQ 12
921 5
922 *
923 CHS
924 R^
925 STO 04
926 RCL 05
927 +
928 STO 08
929 10
930 STO 10
931 +
932 SIN
933 7
934 *
935 +
936 RCL 12
937 ST- 07
938 RCL 13
939 ST+ X
940 STO 11
941 +
942 20
943 -
944 SIN
945 ST+ X
946 RCL 07
947 28
948 +
949 COS
950 +
951 6
952 *
953 +
954 RCL 13
955 RCL 10
956 +
957 SIN
958 RCL 10
959 ST- 06
960 *
961 -
962 RCL 06
963 RCL 12
964 -
965 COS
966 RCL 12
967 RCL 04
968 -
969 COS
970 -
971 15
972 *
973 +
974 RCL 12
975 RCL 03
976 -
977 51
978 -
979 COS
980 ST+ Y
981 RCL 05
982 RCL 12
983 +
984 STO 08
985 11
986 -
987 COS
988 +
989 RCL 05
990 62.2
991 -
992 SIN
993 22
994 *
995 +
996 39
997 ST+ 04
998 *
999 +
1000 RCL 03
1001 73
1002 +
1003 SIN
1004 RCL 04
1005 SIN
1006 +
1007 35
1008 *
1009 -
1010 RCL 11
1011 19
1012 +
1013 STO 11
1014 SIN
1015 209
1016 *
1017 -
1018 STO 03
1019 RCL 11
1020 14
1021 +
1022 COS
1023 RCL 08
1024 SIN
1025 -
1026 6
1027 *
1028 STO 04
1029 RCL 11        
1030 COS
1031 34
1032 *
1033 STO 05
1034 429864
1035 STO 06
1036 19192
1037 STO 09
1038 45955
1039 CHS
1040 STO 12
1041 4630
1042 STO 08
1043 3
1044 CHS
1045 STO 13
1046 773
1047 RCL 00
1048 +
1049 STO 11
1050 1486
1051 STO 07
1052 172993
1053 STO 10
1054 74005
1055 521
1056 XEQ 14
1057 RTN
1058 LBL 08
1059 XEQ 12
1060 22
1061 *
1062 R^
1063 STO 04
1064 50
1065 -
1066 COS
1067 24
1068 *
1069 +
1070 RCL 06
1071 RCL 13
1072 ST+ 04
1073 RCL 05
1074 +
1075 ST- 07
1076 -
1077 COS
1078 9
1079 *
1080 +
1081 RCL 04
1082 77
1083 -
1084 COS
1085 RCL 07
1086 16
1087 +
1088 COS
1089 +
1090 5
1091 *
1092 +
1093 RCL 13
1094 7
1095 +
1096 SIN
1097 71
1098 *
1099 +
1100 RCL 05
1101 27.8
1102 +
1103 COS
1104 584
1105 *
1106 +
1107 STO 03
1108 RCL 04
1109 69
1110 -
1111 SIN
1112 RCL 13
1113 36
1114 -
1115 COS
1116 -
1117 9
1118 *
1119 STO 04
1120 RCL 13
1121 6
1122 +
1123 COS
1124 17
1125 *
1126 CHS
1127 STO 05
1128 219883
1129 STO 06
1130 30070
1131 STO 09
1132 899
1133 STO 08
1134 CLX
1135 STO 13
1136 1770
1137 RCL 00
1138 9
1139 *
1140 -
1141 STO 11
1142 1426
1143 STO 07
1144 55654
1145 CHS
1146 STO 12
1147 48124
1148 STO 10
1149 83660
1150 +
1151 1102
1152 XEQ 14
1153 RTN
1154 LBL 09
1155 RCL 00
1156 2890
1157 *
1158 65
1159 +
1160 COS
1161 7
1162 *
1163 STO 03
1164 RCL 00
1165 STO 05
1166 ST+ 05
1167 9
1168 *
1169 4
1170 -
1171 STO 04
1172 39489
1173 STO 09
1174 24900
1175 STO 08
1176 CLX
1177 STO 13
1178 17140
1179 STO 11
1180 224075
1181 STO 10
1182 238927
1183 STO 12
1184 146567
1185 STO 06
1186 36260
1187 -
1188 1397
1189 STO 07
1190 GTO 14
1191 LBL 12
1192 102
1193 RCL 00
1194 2.06
1195 *
1196 STO 04
1197 *
1198 STO 13
1199 4
1200 LASTX
1201 *
1202 STO 05
1203 52
1204 *
1205 STO 12
1206 19
1207 RCL 04
1208 *
1209 STO 11
1210 593
1211 LASTX
1212 *
1213 STO 07
1214 LASTX
1215 880
1216 *
1217 STO 08
1218 +
1219 STO 06
1220 LASTX
1221 RCL 07
1222 -
1223 STO 09
1224 +
1225 STO 10
1226 RCL 05
1227 ST+ X
1228 RCL 04
1229 31
1230 *
1231 STO 03
1232 RCL 13
1233 RCL 05
1234 -
1235 11
1236 -
1237 COS
1238 RTN
1239 LBL 13
1240 0
1241 STO 03
1242 STO 04
1243 STO 05
1244 RDN
1245 LBL 14
1246 RCL 00        
1247 ST* 06
1248 ST* 07
1249 ST* 13
1250 *
1251 +
1252 X<> 13
1253 ST+ 08
1254 5
1255 10^X
1256 ST/ 08
1257 RCL 06
1258 RCL 12
1259 +
1260 RCL 07
1261 RCL 10
1262 +
1263 9
1264 STO 06
1265 SQRT
1266 10^X
1267 ST/ 03
1268 ST/ 04
1269 ST/ 05
1270 ST/ 09
1271 ST/ 11
1272 ST/ 13
1273 ST/ Z
1274 /
1275 STO 10
1276 -
1277 ENTER
1278 ENTER
1279 LBL 11
1280 SIN
1281 RCL 08
1282 R-D
1283 *
1284 +
1285 DSE 06
1286 GTO 11
1287 STO 07
1288 2
1289 /
1290 TAN
1291 RCL 08
1292 1
1293 ST+ Y
1294 RCL 08
1295 -
1296 /
1297 SQRT
1298 *
1299 ATAN
1300 ST+ X
1301 1
1302 RCL 07
1303 COS
1304 RCL 08
1305 *
1306 -
1307 RCL 09
1308 *
1309 RCL 05
1310 +
1311 X<> 10
1312 +
1313 RCL 13
1314 ST+ 03
1315 -
1316 COS
1317 RCL 11
1318 LASTX
1319 SIN
1320 P-R
1321 X<>Y
1322 ASIN
1323 RCL 04
1324 +
1325 STO 09
1326 X<> Z
1327 R-P
1328 CLX
1329 RCL 03
1330 +
1331 STO 08
1332 X<>Y
1333 RCL 10
1334 P-R
1335 RCL 08
1336 X<>Y
1337 P-R
1338 RCL 02
1339 ST+ Z
1340 CLX
1341 RCL 01
1342 +
1343 R-P
1344 X<>Y
1345 STO 03
1346 STO 13
1347 RDN
1348 R-P
1349 STO 05
1350 X<>Y
1351 STO 04
1352 RCL 10
1353 SQRT
1354 1/X
1355 RCL 10
1356 ST/ Y
1357 X^2
1358 1
1359 -
1360 RCL 05
1361 X^2
1362 ST+ Y
1363 ST+ X
1364 /
1365 ST* Y
1366 -
1367 RCL 05
1368 176
1369 /
1370 ST* Y
1371 +
1372 ST- 03
1373 LBL 00
1374 RCL 00
1375 1934
1376 *
1377 55
1378 +
1379 1
1380 P-R
1381 391
1382 /
1383 549.4
1384 SQRT
1385 -
1386 RCL 00
1387 77
1388 /
1389 +
1390 STO 12
1391 COS
1392 X<>Y
1393 209
1394 /
1395 *
1396 X<> 12
1397 RCL 03
1398 LASTX
1399 -
1400 STO 03
1401 RCL 04
1402 XEQ 13
1403 STO 06
1404 RCL 16
1405 R^
1406 STO 07
1407 90
1408 STO 11
1409 ST- Z
1410 R^
1411 +
1412 RCL 12
1413 +
1414 RCL 14
1415 -
1416 X<>Y
1417 XEQ 13
1418 ST- 11
1419 RCL Z
1420 ENTER
1421 COS
1422 RCL 05
1423 409
1424 *
1425 /
1426 -
1427 STO 12
1428 10.3
1429 RCL 12
1430 5.11
1431 +
1432 /
1433 +
1434 TAN
1435 1/X
1436 59
1437 /
1438 +
1439 X<>Y
1440 RCL 11        
1441 RTN
1442 LBL 13
1443 SIN
1444 X<>Y
1445 LASTX
1446 COS
1447 P-R
1448 RDN
1449 R-P
1450 X<> Z
1451 -
1452 X<>Y
1453 P-R
1454 R^
1455 R-P
1456 X<>Y
1457 RDN
1458 R-P
1459 R^
1460 END
 
 
    ( 2117 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.890
                                RDN      h  =    10.519
                                RDN     h0 =    10.605

  R03 =   0.678 = Geocentric longitude                   R06 = 1.125 = right-ascension            R08 = -27.881 = heliocentric longitude
  R04 =  -1.264   = Geocentric latitude                   R07 = -0.891  = 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.990
                                  RDN      h  =  13.916
                                  RDN     h0 =  13.981
  R03 = -13.318 = Geocentric longitude                   R06 = -11.568 = right-ascension            R08 = -18.191 = heliocentric longitude
  R04 = -1.752  = Geocentric latitude                      R07 =  -6.872  = 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.

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

  Mars:     0°0114
  Jupiter:   0°0128
  Saturn:   0°0115
  Uranus:  0°0107
Neptune:  0°0066
  Pluto:     0°0110


      b)  Program2  [ 1900/03/01 , 2100/02/28 ] ( 1945 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 X<>Y
1195 RCL 10
1196 P-R
1197 RCL 08
1198 X<>Y
1199 P-R
1200 RCL 02
1201 ST+ Z
1202 CLX
1203 RCL 01
1204 +
1205 R-P
1206 X<>Y
1207 STO 03
1208 RDN
1209 R-P
1210 STO 05
1211 X<>Y
1212 STO 04
1213 LBL 00
1214 RCL 00
1215 77
1216 /
1217 549.4
1218 SQRT
1219 -
1220 RCL 03
1221 RCL 04
1222 XEQ 13
1223 STO 06
1224 RCL 16
1225 R^
1226 STO 07
1227 90
1228 STO 11
1229 ST- Z
1230 R^
1231 +
1232 RCL 14
1233 -
1234 X<>Y
1235 XEQ 13
1236 ST- 11
1237 RCL Z
1238 ENTER
1239 COS
1240 RCL 05
1241 409
1242 *
1243 /
1244 -
1245 STO 12
1246 10.3
1247 RCL 12
1248 5.11
1249 +
1250 /
1251 +
1252 TAN
1253 1/X
1254 59
1255 /
1256 +
1257 X<>Y
1258 RCL 11
1259 RTN
1260 LBL 13
1261 SIN
1262 X<>Y
1263 LASTX
1264 COS
1265 P-R
1266 RDN
1267 R-P
1268 X<> Z
1269 -
1270 X<>Y
1271 P-R
1272 R^
1273 R-P
1274 X<>Y
1275 RDN
1276 R-P
1277 R^
1278 RTN
1279 LBL 12
1280 102
1281 RCL 00
1282 2.06
1283 *
1284 STO 04
1285 *
1286 STO 13
1287 4
1288 LASTX
1289 *
1290 STO 05
1291 52
1292 *
1293 STO 12
1294 19
1295 RCL 04
1296 *
1297 STO 11
1298 593
1299 LASTX
1300 *
1301 STO 07
1302 LASTX
1303 880
1304 *
1305 STO 08
1306 +
1307 STO 06
1308 LASTX
1309 RCL 07        
1310 -
1311 STO 09
1312 +
1313 STO 10
1314 RCL 05
1315 ST+ X
1316 RCL 04
1317 31
1318 *
1319 STO 03
1320 RCL 13
1321 RCL 05
1322 -
1323 11
1324 -
1325 COS
1326 END
 
 
    ( 1945 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.


      c)  Program3  [ 1900/03/01 , 2100/02/28 ] ( 1821 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 X<>Y
1105 RCL 10
1106 P-R
1107 RCL 08
1108 X<>Y
1109 P-R
1110 RCL 02
1111 ST+ Z
1112 CLX
1113 RCL 01
1114 +
1115 R-P
1116 X<>Y
1117 STO 03
1118 RDN
1119 R-P
1120 STO 05
1121 X<>Y
1122 STO 04
1123 LBL 00
1124 RCL 00
1125 77
1126 /
1127 549.4
1128 SQRT
1129 -
1130 RCL 03
1131 RCL 04
1132 XEQ 13
1133 STO 06
1134 RCL 15
1135 R^
1136 STO 07
1137 90
1138 STO 11
1139 ST- Z
1140 R^
1141 +
1142 RCL 13
1143 -
1144 X<>Y
1145 XEQ 13
1146 ST- 11
1147 RCL Z
1148 ENTER
1149 COS
1150 RCL 05
1151 409
1152 *
1153 /
1154 -
1155 STO 12
1156 10.3
1157 RCL 12
1158 5.11
1159 +
1160 /
1161 +
1162 TAN
1163 1/X
1164 59
1165 /
1166 +
1167 X<>Y
1168 RCL 11
1169 RTN
1170 LBL 13
1171 SIN
1172 X<>Y
1173 LASTX
1174 COS
1175 P-R
1176 RDN
1177 R-P
1178 X<> Z
1179 -
1180 X<>Y
1181 P-R
1182 R^
1183 R-P
1184 X<>Y
1185 RDN
1186 R-P
1187 R^
1188 RTN
1189 LBL 12
1190 RCL 00
1191 2.06
1192 *
1193 STO 04
1194 4
1195 *
1196 STO 05
1197 52
1198 *
1199 STO 12
1200 19
1201 RCL 04
1202 *
1203 STO 11
1204 593
1205 LASTX
1206 *
1207 STO 07
1208 LASTX
1209 880
1210 *
1211 STO 08
1212 +
1213 STO 06
1214 LASTX
1215 RCL 07        
1216 -
1217 STO 09
1218 +
1219 STO 10
1220 RCL 04
1221 31
1222 *
1223 STO 03
1224 RCL 04
1225 102
1226 *
1227 END
 
 
    ( 1821 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 ] ( 1496 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
  Pluto:     0°0159


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

  R03 =   0.694   = Geocentric longitude                 R06 = 1.139 = right-ascension            R08 = -27.871 = heliocentric longitude
  R04 =  -1.264   = Geocentric latitude                   R07 = -0.884 = declination                 R09 = -1.807   = heliocentric latitude
  R05 = 1.9746 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.891 AU = radius vector


PLUTO


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

  R03 = -57.901 = Geocentric longitude                    R06 =  -54.905 = right-ascension             R08 =  -59.530 = heliocentric longitude
  R04 = -3.026  = Geocentric latitude                       R07 =  -22.637  = 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.


      e)  Program5 (?)  [ 1900/03/01 , 2100/02/28 ] ( 1083 bytes )


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

-Even less decimals... and even less precision !

-Maximum errors with heliocentric longitudes:

  Mars:     0°035
  Jupiter:   0°084
  Saturn:   0°070
  Uranus:  0°072
Neptune:  0°053
  Pluto:     0°046


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 360007.7
 57 STO 06
 58 ST+ 06
 59 *
 60 10
 61 /
 62 RCL 09
 63 +
 64 31
 65 1/X
 66 -
 67 RCL 13
 68 15
 69 *
 70 +
 71 RCL 14
 72 +
 73 STO 13
 74 167
 75 STO 08
 76 34.4
 77 STO 07
 78 1600.6
 79 CHS
 80 STO 12
 81 59.4
 82 +
 83 STO 10
 84 10
 85 ST/ 09
 86 CLST
 87 STO 01
 88 STO 02
 89 STO 11
 90 XEQ 13
 91 CLX
 92 RCL 03          
 93 RCL 05
 94 P-R
 95 STO 01
 96 RDN
 97 STO 02
 98 CLX
 99 RCL 11
100 RTN
101 LBL 03
102 RCL 00
103 3599905
104 %
105 3
106 -
107 STO 11
108 STO 12
109 RCL 00
110 44526711
111 %
112 68.25
113 -
114 STO 03
115 ST+ X
116 STO 05
117 STO 06
118 STO 10
119 ST- 11
120 RCL 00
121 47719887
122 %
123 51.57
124 -
125 ST- 05
126 STO 07
127 ST+ 10
128 RCL 00
129 48320202
130 %
131 86.66
132 +
133 STO 09
134 COS
135 RCL 07
136 SIN
137 *
138 5.6
139 *
140 CHS
141 RCL 06
142 RCL 09
143 -
144 SIN
145 3
146 SQRT
147 *
148 +
149 RCL 05
150 SIN
151 RCL 09
152 COS
153 *
154 -
155 RCL 09
156 SIN
157 51.3
158 *
159 +
160 STO 04
161 RCL 07
162 ST+ X
163 STO 08
164 SIN
165 ST+ X
166 RCL 06
167 SIN
168 6.6
169 *
170 +
171 RCL 09
172 ST+ X
173 SIN
174 -
175 RCL 07
176 SIN
177 63
178 *
179 -
180 RCL 06
181 RCL 08
182 -
183 SIN
184 RCL 05
185 RCL 12
186 -
187 SIN
188 -
189 RCL 10          
190 SIN
191 -
192 RCL 11
193 SIN
194 -
195 RCL 12
196 RCL 07
197 -
198 SIN
199 +
200 2
201 /
202 +
203 RCL 05
204 SIN
205 12.7
206 *
207 -
208 RCL 12
209 SIN
210 ST+ X
211 -
212 4812679
213 RCL 00
214 *
215 2117.3
216 +
217 +
218 10
219 ST/ 04
220 /
221 360
222 MOD
223 STO 03
224 26
225 RCL 07
226 COS
227 2
228 SQRT
229 *
230 -
231 4
232 10^X
233 /
234 STO 05
235 XEQ 00
236 RTN
237 LBL 01
238 3.87
239 STO 09
240 2056
241 STO 08
242 140
243 STO 11
244 31
245 STO 07
246 1549
247 STO 10
248 5004
249 STO 12
250 2989481
251 STO 06
252 967
253 24
254 XEQ 13
255 RTN
256 LBL 02
257 7.233
258 STO 09
259 68
260 STO 08
261 68
262 STO 11
263 28
264 STO 07
265 2631
266 STO 10
267 3624
268 STO 12
269 1170384
270 STO 06
271 1534
272 18
273 XEQ 13
274 RTN
275 LBL 04
276 15.24
277 STO 09
278 934
279 STO 08
280 37
281 STO 11
282 STO 07
283 96.4
284 CHS
285 STO 12
286 382834
287 STO 06
288 479
289 CHS
290 STO 10          
291 991
292 15
293 XEQ 13
294 RTN
295 LBL 05
296 XEQ 12
297 RCL 11
298 7
299 -
300 SIN
301 7
302 *
303 CHS
304 RCL 10
305 33
306 -
307 SIN
308 RCL 06
309 RCL 11
310 -
311 60
312 -
313 SIN
314 +
315 RCL 09
316 76
317 -
318 SIN
319 +
320 +
321 STO 03
322 52.03
323 STO 09
324 26
325 STO 11
326 60726
327 STO 06
328 485
329 STO 08
330 32
331 STO 07
332 686
333 STO 12
334 287
335 STO 10
336 2009
337 20
338 XEQ 14
339 RTN
340 LBL 06
341 XEQ 12
342 RCL 11
343 7
344 -
345 SIN
346 16
347 *
348 RCL 09
349 ST+ X
350 30
351 +
352 COS
353 5
354 *
355 +
356 RCL 09
357 14
358 +
359 COS
360 ST+ X
361 +
362 RCL 07
363 RCL 11
364 +
365 49
366 -
367 SIN
368 +
369 STO 03
370 24470
371 STO 06
372 95
373 STO 09
374 555
375 STO 08
376 50
377 STO 11
378 39
379 STO 07
380 1001
381 STO 12
382 1861
383 STO 10          
384 2273
385 18
386 XEQ 14
387 RTN
388 LBL 07
389 XEQ 12
390 RCL 05
391 62
392 -
393 SIN
394 18
395 *
396 RCL 03
397 ST+ X
398 19
399 +
400 SIN
401 4
402 *
403 -
404 RCL 04
405 177
406 *
407 51
408 -
409 COS
410 +
411 RCL 03
412 RCL 05
413 +
414 ST+ X
415 11
416 -
417 COS
418 +
419 STO 03
420 8597
421 STO 06
422 192
423 STO 09
424 919
425 CHS
426 STO 12
427 463
428 STO 08
429 15
430 STO 11
431 30
432 STO 07
433 3460
434 STO 10
435 1480
436 10
437 XEQ 14
438 RTN
439 LBL 08
440 XEQ 12
441 RCL 05
442 28
443 +
444 COS
445 12
446 *
447 RCL 03
448 7
449 +
450 SIN
451 2
452 SQRT
453 *
454 +
455 STO 03
456 4398
457 STO 06
458 301
459 STO 09
460 90
461 STO 08
462 35
463 STO 11
464 29
465 STO 07
466 1113
467 CHS
468 STO 12
469 962
470 STO 10
471 2636
472 22
473 XEQ 14
474 RTN
475 LBL 09
476 395
477 STO 09
478 2490
479 STO 08
480 4481
481 STO 10          
482 343
483 STO 11
484 4778
485 STO 12
486 2931
487 STO 06
488 2206
489 28
490 STO 07
491 LBL 13
492 0
493 STO 03
494 RDN
495 LBL 14
496 RCL 00
497 ST* 06
498 ST* 07
499 *
500 +
501 X<> 12
502 9
503 X<> 06
504 +
505 RCL 07
506 RCL 10
507 +
508 4
509 10^X
510 ST/ 08
511 SQRT
512 SQRT
513 ST/ 09
514 ST+ X
515 ST/ 03
516 ST/ 11
517 ST/ 12
518 ST/ Z
519 /
520 STO 10
521 -
522 ENTER
523 ENTER
524 LBL 11
525 SIN
526 RCL 08
527 R-D
528 *
529 +
530 DSE 06
531 GTO 11
532 STO 07
533 2
534 /
535 TAN
536 RCL 08
537 1
538 ST+ Y
539 RCL 08
540 -
541 /
542 SQRT
543 *
544 ATAN
545 ST+ X
546 1
547 RCL 07
548 COS
549 RCL 08
550 *
551 -
552 RCL 09
553 *
554 X<> 10
555 +
556 RCL 12
557 ST+ 03
558 -
559 COS
560 RCL 11
561 LASTX
562 SIN
563 P-R
564 X<>Y
565 ASIN
566 STO 09
567 X<> Z
568 R-P
569 CLX
570 RCL 03
571 +
572 STO 08
573 X<>Y
574 RCL 10
575 P-R
576 RCL Z
577 X<>Y
578 P-R
579 RCL 02          
580 ST+ Z
581 CLX
582 RCL 01
583 +
584 R-P
585 X<>Y
586 STO 03
587 RDN
588 R-P
589 STO 05
590 X<>Y
591 STO 04
592 LBL 00
593 RCL 00
594 77
595 /
596 549.4
597 SQRT
598 -
599 RCL 03
600 RCL 04
601 XEQ 13
602 STO 06
603 RCL 15
604 R^
605 STO 07
606 90
607 STO 11
608 ST- Z
609 R^
610 +
611 RCL 13
612 -
613 X<>Y
614 XEQ 13
615 ST- 11
616 RCL Z
617 ENTER
618 COS
619 RCL 05
620 409
621 *
622 /
623 -
624 STO 12
625 10.3
626 RCL 12
627 5.11
628 +
629 /
630 +
631 TAN
632 1/X
633 59
634 /
635 +
636 X<>Y
637 RCL 11
638 RTN
639 LBL 13
640 SIN
641 X<>Y
642 LASTX
643 COS
644 P-R
645 RDN
646 R-P
647 X<> Z
648 -
649 X<>Y
650 P-R
651 R^
652 R-P
653 X<>Y
654 RDN
655 R-P
656 R^
657 RTN
658 LBL 12
659 4
660 RCL 00
661 2.06
662 *
663 STO 04
664 *
665 STO 05
666 19
667 RCL 04
668 *
669 STO 11
670 593
671 LASTX
672 *
673 STO 07
674 LASTX
675 880
676 *
677 STO 08
678 +
679 STO 06
680 LASTX
681 RCL 07          
682 -
683 STO 09
684 +
685 STO 10
686 102
687 RCL 04
688 *
689 STO 03
690 END
 
 
    ( 1083 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.29
                                  RDN      h  =  15.39
                                  RDN     h0 =  15.45
                  

  R03 = 312.24 = Geocentric longitude                R06  = -43.88 = right-ascension
  R04 =  -4.69   = Geocentric latitude                   R07 =  -21.62  = declination
  R05 = 0.00271 AU = distance Earth-Moon


MERCURY


   R/S  or  XEQ 01   >>>>    Az =   267.52
                                RDN      h  =    3.97
                                RDN     h0 =    4.16
                              

  R03 =  17.43 = Geocentric longitude                  R06 = 17.00 = right-ascension           R08 = -100.01 = heliocentric longitude
  R04 = -2.42  = Geocentric latitude                     R07 =  4.61  = declination                  R09 = -3.65      = heliocentric latitude
  R05 = 0.7052 AU = distance Earth-Mercury                                                                R10 = 0.4665  AU = radius vector


VENUS


   R/S  or  XEQ 02    >>>>     Az =   252.96
                                 RDN      h  =   -0.25
                                 RDN     h0 =     0.27

  R03 = 32.54 = Geocentric longitude                 R06 =  30.77 = right-ascension          R08 = 19.84 = heliocentric longitude
  R04 = -1.22  = Geocentric latitude                   R07 = 11.21  = declination                 R09 =  -2.85  = heliocentric latitude
  R05 = 1.702 AU = distance Earth-Venus                                                                   R10 =  0.725 AU = radius vector


MARS


   R/S  or  XEQ 04   >>>>    Az =  -76.90
                                RDN      h  =    10.52
                                RDN     h0 =    10.60

  R03 =   0.68   = Geocentric longitude                 R06 = 1.13 = right-ascension            R08 = -27.88 = heliocentric longitude
  R04 =  -1.26   = Geocentric latitude                   R07 = -0.89 = declination                 R09 = -1.81   = heliocentric latitude
  R05 = 1.975 AU = distance Earth-Mars                                                                      R10 = 1.382 AU = radius vector


JUPITER


   R/S  or  XEQ 05     >>>>    Az = 232.44
                                  RDN      h  = -7.53
                                  RDN     h0 = -7.61   ( often meaningless when < 0 )

  R03 = 54.25 = Geocentric longitude                   R06 = 52.08  = right-ascension             R08 = 56.77 = heliocentric longitude
  R04 = -0.75   = Geocentric latitude                    R07 =  18.10  = declination                   R09 = -0.90  = heliocentric latitude
  R05 =  5.989 AU = distance Earth-Jupiter                                                                       R10 = 5.010 AU = radius vector


SATURN


   R/S  or  XEQ 06     >>>>   Az = -63.02
                                  RDN      h  =  13.91
                                  RDN     h0 =  13.97

  R03 = -13.29 = Geocentric longitude                   R06 = -11.54 = right-ascension            R08 = -18.18 = heliocentric longitude
  R04 = -1.76  = Geocentric latitude                      R07 =  -6.86  = declination                   R09 = -1.85  = heliocentric latitude
  R05 = 10.213 AU = distance Earth-Saturn                                                                      R10 = 9.670 AU = radius vector


URANUS


   R/S  or  XEQ 07   >>>>   Az =   233.95
                                RDN      h  =  -6.44
                                RDN     h0 =  -6.51   ( often meaningless when < 0 )

  R03 = 52.48 = Geocentric longitude                    R06 =  50.14 = right-ascension              R08 =  53.03 = heliocentric longitude
  R04 = -0.25  = Geocentric latitude                       R07 =  18.15 = declination                   R09 = -0.26  = heliocentric latitude
  R05 = 20.594 AU = distance Earth-Uranus                                                                      R10 = 19.605 AU = radius vector


NEPTUNE


   R/S  or  XEQ 08   >>>>   Az =   -75.21
                                 RDN      h  =  11.10
                                 RDN     h0 =  11.18

  R03 = -1.07  = Geocentric longitude                    R06 =  -0.50 = right-ascension                R08 = -2.38 = heliocentric longitude
  R04 = -1.21  = Geocentric latitude                       R07 =  -1.53   = declination                    R09 = -1.24 = heliocentric latitude
  R05 = 30.662 AU = distance Earth-Neptune                                                                      R10 = 29.930 AU = radius vector


PLUTO


   R/S  or  XEQ 09     >>>>   Az = -16.69
                                  RDN      h  =  18.07
                                  RDN     h0 =  18.12

  R03 = -57.94 = Geocentric longitude                    R06 =  -54.94 = right-ascension             R08 =  -59.57 = heliocentric longitude
  R04 = -3.02  = Geocentric latitude                       R07 =  -22.64  = declination                    R09 = -3.00  = heliocentric latitude
  R05 = 34.833 AU = distance Earth-Pluto                                                                            R10 = 35.014 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 + 1503 or 1499 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 X<>Y
888 RCL 10
889 P-R
890 RCL 08
891 X<>Y
892 P-R
893 RCL 02
894 ST+ Z
895 CLX
896 RCL 01
897 +
898 R-P
899 X<>Y
900 STO 03
901 STO 13
902 RDN
903 R-P
904 STO 05
905 X<>Y
906 STO 04
907 RCL 10
908 SQRT
909 1/X
910 RCL 10          
911 ST/ Y
912 X^2
913 1
914 -
915 RCL 05
916 X^2
917 ST+ Y
918 ST+ X
919 /
920 ST* Y
921 -
922 RCL 05
923 176
924 /
925 ST* Y
926 +
927 ST- 03
928 LBL "Y"
929 RCL 00
930 1934
931 *
932 55
933 +
934 1
935 P-R
936 391
937 /
938 549.4
939 SQRT
940 -
941 RCL 00
942 77
943 /
944 +
945 STO 12
946 COS
947 X<>Y
948 209
949 /
950 *
951 X<> 12
952 RCL 03
953 LASTX
954 -
955 STO 03
956 RCL 04
957 XEQ 13
958 STO 06
959 RCL 16
960 R^
961 STO 07
962 90
963 STO 11
964 ST- Z
965 R^
966 +
967 RCL 12
968 +
969 RCL 14
970 -
971 X<>Y
972 XEQ 13
973 ST- 11
974 RCL Z
975 ENTER
976 COS
977 RCL 05
978 409
979 *
980 /
981 -
982 STO 12
983 10.3
984 RCL 12
985 5.11
986 +
987 /
988 +
989 TAN
990 1/X
991 59
992 /
993 +
994 X<>Y
995 RCL 11
996 RTN
997 LBL 13
998 SIN
999 X<>Y
1000 LASTX
1001 COS
1002 P-R
1003 RDN
1004 R-P
1005 X<> Z
1006 -
1007 X<>Y
1008 P-R
1009 R^
1010 R-P
1011 X<>Y
1012 RDN
1013 R-P
1014 R^
1015 RTN
1016 LBL "Z"
1017 102
1018 RCL 00
1019 2.06
1020 *
1021 STO 04
1022 *
1023 STO 13
1024 4
1025 LASTX
1026 *
1027 STO 05
1028 52
1029 *
1030 STO 12
1031 19
1032 RCL 04
1033 *
1034 STO 11
1035 593
1036 LASTX
1037 *
1038 STO 07
1039 LASTX
1040 880
1041 *
1042 STO 08
1043 +
1044 STO 06
1045 LASTX
1046 RCL 07
1047 -
1048 STO 09
1049 +
1050 STO 10
1051 RCL 05
1052 ST+ X
1053 RCL 04
1054 31
1055 *
1056 STO 03
1057 RCL 13
1058 RCL 05        
1059 -
1060 11
1061 -
1062 COS
1063 END
 
 
    ( 1503 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 X<>Y
862 RCL 10
863 P-R
864 RCL 08
865 X<>Y
866 P-R
867 RCL 02
868 ST+ Z
869 CLX
870 RCL 01
871 +
872 R-P
873 X<>Y
874 STO 03
875 STO 13
876 RDN
877 R-P
878 STO 05
879 X<>Y
880 STO 04
881 RCL 10          
882 SQRT
883 1/X
884 RCL 10
885 ST/ Y
886 X^2
887 1
888 -
889 RCL 05
890 X^2
891 ST+ Y
892 ST+ X
893 /
894 ST* Y
895 -
896 RCL 05
897 176
898 /
899 ST* Y
900 +
901 ST- 03
902 LBL 00
903 RCL 00
904 1934
905 *
906 55
907 +
908 1
909 P-R
910 391
911 /
912 549.4
913 SQRT
914 -
915 RCL 00
916 77
917 /
918 +
919 STO 12
920 COS
921 X<>Y
922 209
923 /
924 *
925 X<> 12
926 RCL 03
927 LASTX
928 -
929 STO 03
930 RCL 04
931 XEQ 13
932 STO 06
933 RCL 16
934 R^
935 STO 07
936 90
937 STO 11
938 ST- Z
939 R^
940 +
941 RCL 12
942 +
943 RCL 14
944 -
945 X<>Y
946 XEQ 13
947 ST- 11
948 RCL Z
949 ENTER
950 COS
951 RCL 05
952 409
953 *
954 /
955 -
956 STO 12
957 10.3
958 RCL 12
959 5.11
960 +
961 /
962 +
963 TAN
964 1/X
965 59
966 /
967 +
968 X<>Y
969 RCL 11
970 RTN
971 LBL 13
972 SIN
973 X<>Y
974 LASTX
975 COS
976 P-R
977 RDN
978 R-P
979 X<> Z
980 -
981 X<>Y
982 P-R
983 R^
984 R-P
985 X<>Y
986 RDN
987 R-P
988 R^
989 RTN
990 LBL "Z"
991 102
992 RCL 00
993 2.06
994 *
995 STO 04
996 *
997 STO 13
998 4
999 LASTX
1000 *
1001 STO 05
1002 52
1003 *
1004 STO 12
1005 19
1006 RCL 04
1007 *
1008 STO 11
1009 593
1010 LASTX
1011 *
1012 STO 07
1013 LASTX
1014 880
1015 *
1016 STO 08
1017 +
1018 STO 06
1019 LASTX
1020 RCL 07
1021 -
1022 STO 09
1023 +
1024 STO 10
1025 RCL 05
1026 ST+ X
1027 RCL 04        
1028 31
1029 *
1030 STO 03
1031 RCL 13
1032 RCL 05
1033 -
1034 11
1035 -
1036 COS
1037 END
 
 
    ( 1499 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/