TransNeptunian Objects for the HP-41
Overview
1°) Sedna - Hauméa - Makémaké
- Quaoar - Orcus - Varuna
2°) Altjira - Chaos - Deucalion - Ixion
3°) Salacia - Sila Nunam - Varda
-The positions of Pluto and Eris may be computed approximately by the
routines listed in "Astronomical Ephemeris(V)"
-The program hereunder calculates approximate heliocentric coordinates
of 13 other transneptunian objects over the time-span [1900-2100]
-A polynomial is used for each coordinate and the precision is about
0°010 for the longitudes
-Thus, we don't have to solve the equation of Kepler.
-Do not use these routines before 1900 or after 2100, the results would
rapidly become meaningless !
1°) Sedna - Hauméa - Makémaké
- Quaoar - Orcus - Varuna
Data Registers: • R00 = T in millenia since 2000/01/01 0h TT ( Register R00 is to be initialized before executing these programs )
R06 = L = heliocentric longitude ( deg )
R07 = B = heliocentric latitude ( deg )
referred to the mean ecliptic and equinox of the date
R08 = R = radius vector ( AU )
Flags: /
Subroutines: /
| 01 LBL "SE"
02 RCL 00 03 10 04 * 05 6 06 RCL Y 07 200 08 * 09 + 10 * 11 CHS 12 898 13 + 14 * 15 1081 16 + 17 * 18 735 19 - 20 * 21 3826 22 - 23 * 24 4621 25 - 26 * 27 2013 28 + 29 * 30 20250 31 + 32 * 33 54247 34 + 35 * 36 45876 37 + 38 STO 06 39 CLX 40 87 41 * 42 CHS 43 187 44 + 45 * 46 558 47 + 48 * 49 315 50 - 51 * 52 1487 53 - 54 * 55 999 56 - 57 * 58 1203 59 + 60 * 61 3866 62 + 63 * 64 4285 65 + 66 * 67 1575 68 - 69 * 70 11800 71 - 72 STO 07 73 CLX 74 570 75 * 76 CHS 77 221 78 + 79 * 80 1543 81 + 82 * 83 568 84 - 85 * 86 2249 87 - 88 * 89 639 90 - 91 * 92 2396 93 + 94 * 95 7923 96 + 97 * 98 13130 99 + 100 * 101 34263 102 - 103 * 104 90899 105 XEQ 10 106 RTN 107 LBL "HA" 108 RCL 00 109 10 110 * 111 1468 112 CHS 113 RCL Y 114 318 115 * 116 + 117 * 118 1332 119 - |
120 *
121 6387 122 + 123 * 124 4028 125 + 126 * 127 2644 128 - 129 * 130 12332 131 - 132 * 133 8029 134 + 135 * 136 18884 137 + 138 * 139 96832 140 + 141 * 142 186356 143 + 144 STO 06 145 CLX 146 691 147 * 148 151 149 + 150 * 151 365 152 - 153 * 154 463 155 - 156 * 157 2352 158 + 159 * 160 1313 161 + 162 * 163 5954 164 - 165 * 166 10961 167 - 168 * 169 30359 170 - 171 * 172 17920 173 + 174 * 175 25813 176 + 177 STO 07 178 CLX 179 106 180 * 181 188 182 + 183 * 184 283 185 - 186 * 187 320 188 - 189 * 190 1155 191 + 192 * 193 604 194 + 195 * 196 1332 197 + 198 * 199 429 200 + 201 * 202 14442 203 - 204 * 205 2377 206 - 207 * 208 51409 209 XEQ 10 210 RTN 211 LBL "MK" 212 RCL 00 213 10 214 * 215 661 216 RCL Y 217 37 218 * 219 + 220 * 221 1466 222 + 223 * 224 3414 225 - 226 * 227 5514 228 - 229 * 230 13208 231 + 232 * 233 3940 234 + 235 * 236 10596 237 - 238 * |
239 11812
240 - 241 * 242 104648 243 + 244 * 245 165249 246 + 247 STO 06 248 CLX 249 458 250 * 251 CHS 252 286 253 - 254 * 255 2471 256 + 257 * 258 506 259 - 260 * 261 1316 262 - 263 * 264 3352 265 - 266 * 267 3515 268 + 269 * 270 9308 271 + 272 * 273 39904 274 - 275 * 276 3217 277 + 278 * 279 28938 280 + 281 STO 07 282 CLX 283 114 284 * 285 CHS 286 465 287 + 288 * 289 29 290 + 291 * 292 752 293 - 294 * 295 501 296 + 297 * 298 463 299 - 300 * 301 2007 302 + 303 * 304 2196 305 - 306 * 307 10372 308 - 309 * 310 7513 311 + 312 * 313 51498 314 XEQ 10 315 RTN 316 LBL "QU" 317 RCL 00 318 10 319 * 320 800 321 RCL Y 322 119 323 * 324 - 325 * 326 CHS 327 833 328 - 329 * 330 1939 331 + 332 * 333 3085 334 + 335 * 336 1533 337 - 338 * 339 7803 340 - 341 * 342 688 343 + 344 * 345 12432 346 + 347 * 348 127015 349 + 350 * 351 247166 352 + 353 STO 06 354 CLX 355 278 356 * 357 210 |
358 +
359 * 360 1177 361 - 362 * 363 1440 364 - 365 * 366 1577 367 + 368 * 369 4776 370 + 371 * 372 3728 373 + 374 * 375 9966 376 - 377 * 378 15483 379 - 380 * 381 9137 382 + 383 * 384 6807 385 + 386 STO 07 387 CLX 388 4 389 * 390 21 391 + 392 * 393 53 394 + 395 * 396 CHS 397 82 398 + 399 * 400 74 401 - 402 * 403 628 404 - 405 * 406 458 407 + 408 * 409 2779 410 + 411 * 412 596 413 - 414 * 415 3494 416 - 417 * 418 43530 419 XEQ 10 420 RTN 421 LBL "OR" 422 RCL 00 423 10 424 * 425 21824 426 RCL Y 427 29314 428 * 429 - 430 * 431 100571 432 + 433 * 434 65613 435 - 436 * 437 139320 438 - 439 * 440 72419 441 + 442 * 443 108726 444 + 445 * 446 46643 447 - 448 * 449 45516 450 - 451 * 452 21082 453 + 454 * 455 3427 456 + 457 * 458 6715 459 + 460 * 461 12794 462 - 463 * 464 23195 465 + 466 * 467 4927 468 - 469 * 470 101439 471 + 472 * 473 140767 474 + 475 STO 06 476 CLX |
477 7290
478 * 479 2408 480 + 481 * 482 33892 483 - 484 * 485 11452 486 - 487 * 488 70649 489 + 490 * 491 10655 492 + 493 * 494 72527 495 - 496 * 497 1744 498 + 499 * 500 31147 501 + 502 * 503 3784 504 + 505 * 506 12375 507 - 508 * 509 2490 510 + 511 * 512 6673 513 + 514 * 515 1201 516 - 517 * 518 27101 519 + 520 * 521 21157 522 - 523 * 524 16511 525 - 526 STO 07 527 CLX 528 912 529 * 530 CHS 531 1450 532 - 533 * 534 4263 535 + 536 * 537 3893 538 + 539 * 540 8217 541 - 542 * 543 1770 544 - 545 * 546 7332 547 + 548 * 549 2097 550 - 551 * 552 3085 553 - 554 * 555 1643 556 + 557 * 558 3134 559 + 560 * 561 2927 562 - 563 * 564 2588 565 + 566 * 567 1567 568 - 569 * 570 18657 571 - 572 * 573 7508 574 + 575 * 576 47324 577 XEQ 10 578 RTN 579 LBL "VA" 580 RCL 00 581 10 582 * 583 1547 584 CHS 585 RCL Y 586 2055 587 * 588 + 589 * 590 4205 591 - 592 * 593 5273 594 + 595 * |
596 1753
597 - 598 * 599 1383 600 - 601 * 602 9237 603 + 604 * 605 14807 606 - 607 * 608 1276 609 - 610 * 611 20200 612 + 613 * 614 14080 615 - 616 * 617 123258 618 + 619 * 620 98643 621 + 622 STO 06 623 CLX 624 924 625 * 626 2658 627 - 628 * 629 728 630 - 631 * 632 7910 633 + 634 * 635 5975 636 - 637 * 638 2200 639 - 640 * 641 12494 642 + 643 * 644 27502 645 - 646 * 647 5519 648 - 649 * 650 37676 651 + 652 * 653 416 654 + 655 STO 07 656 CLX 657 18434 658 * 659 7536 660 - 661 * 662 74208 663 - 664 * 665 26869 666 + 667 * 668 120019 669 + 670 * 671 37271 672 - 673 * 674 99388 675 - 676 * 677 25233 678 + 679 * 680 44620 681 + 682 * 683 8579 684 - 685 * 686 11017 687 - 688 * 689 2255 690 + 691 * 692 2053 693 + 694 * 695 4324 696 - 697 * 698 450 699 - 700 * 701 5175 702 + 703 * 704 42987 705 LBL 10 706 + 707 E3 708 ST/ 06 709 ST/ 07 710 / 711 STO 08 712 RCL 07 713 RCL 06 714 END |
( 1465
bytes / SIZE 009 )
| STACK | INPUTS | OUTPUTS |
| Z | / | R |
| Y | / | B |
| X | / | L |
L = heliocentric longitude ( deg )
where B = heliocentric
latitude ( deg )
referred to the mean ecliptic and equinox of the date
R = radius vector ( AU )
Example: T = 0.1 i-e 2100/01/01 0 h TT
>>> 0.1 STO 00
• SEDNA
XEQ "SE" >>>> L = 114°977
RDN B = -6°164
RDN R = 77.823 AU
• HAUMEA
XEQ "HA" ( or R/S ) >>>>
L = 303°058
RDN B = 0°138
RDN R = 37.801 AU
• MAKEMAKE
XEQ "MK" ( or R/S ) >>>> L = 257°873
RDN B = 1°627
RDN R = 48.116 AU
• QUAOAR
XEQ "QU" ( or R/S ) >>>>
L = 381°475 = 21°475 ( add 360 MOD
after line 419 if you prefer a longitude between 0° and 360° )
RDN B = -1°553
RDN R = 41.979 AU
• ORCUS
XEQ "OR" ( or R/S ) >>>> L = 256°038
RDN B = -5°174
RDN R = 37.003 AU
• VARUNA
XEQ "VA" ( or R/S ) >>>> L = 219°615
RDN B = 14°838
RDN R = 44.872
2°) Altjira - Chaos - Deucalion - Ixion
Data Registers: • R00 = T in millenia since 2000/01/01 0h TT ( Register R00 is to be initialized before executing these programs )
R06 = L = heliocentric longitude ( deg )
R07 = B = heliocentric latitude ( deg )
referred to the mean ecliptic and equinox of the date
R08 = R = radius vector ( AU )
Flags: /
Subroutines: /
| 01 LBL "ALT"
02 RCL 00 03 10 04 * 05 1433 06 RCL Y 07 468 08 * 09 - 10 * 11 1249 12 + 13 * 14 3248 15 - 16 * 17 1612 18 - 19 * 20 2760 21 + 22 * 23 1710 24 + 25 * 26 2790 27 - 28 * 29 2324 30 + 31 * 32 5097 33 + 34 * 35 13535 36 - 37 * 38 119530 39 + 40 * 41 50748 42 + 43 STO 06 44 RDN 45 ST+ X 46 33 47 - 48 * 49 203 50 + 51 * 52 107 53 - 54 * 55 617 56 - 57 * 58 1050 59 + 60 * 61 104 62 - 63 * 64 1812 65 - 66 * 67 4081 68 + 69 * 70 2675 71 - 72 * 73 9165 74 - 75 * 76 7045 77 + 78 * 79 3922 80 + 81 STO 07 82 CLX 83 1501 84 * 85 CHS 86 33 87 - 88 * 89 4552 90 + 91 * 92 11 93 - 94 * 95 4965 96 - 97 * 98 192 99 + 100 * 101 1999 102 + 103 * 104 364 105 + 106 * 107 875 108 + 109 * 110 3770 |
111 -
112 * 113 2004 114 - 115 * 116 5552 117 + 118 * 119 45 E3 120 XEQ 10 121 RTN 122 LBL "CHA" 123 RCL 00 124 10 125 * 126 428 127 RCL Y 128 3066 129 * 130 + 131 * 132 CHS 133 13020 134 + 135 * 136 4260 137 + 138 * 139 23333 140 - 141 * 142 12636 143 - 144 * 145 21289 146 + 147 * 148 20113 149 + 150 * 151 5001 152 - 153 * 154 18129 155 - 156 * 157 14955 158 - 159 * 160 370 161 - 162 * 163 24734 164 + 165 * 166 133475 167 + 168 * 169 58063 170 + 171 STO 06 172 CLX 173 1039 174 * 175 CHS 176 1546 177 + 178 * 179 5606 180 + 181 * 182 4316 183 - 184 * 185 13689 186 - 187 * 188 1353 189 + 190 * 191 19312 192 + 193 * 194 11682 195 + 196 * 197 14582 198 - 199 * 200 28876 201 - 202 * 203 263 204 + 205 * 206 27872 207 + 208 * 209 1718 210 + 211 STO 07 212 CLX 213 1278 214 * 215 CHS 216 573 217 + 218 * 219 3656 220 + |
221 *
222 1796 223 - 224 * 225 3518 226 - 227 * 228 2977 229 + 230 * 231 1434 232 + 233 * 234 4654 235 - 236 * 237 2546 238 - 239 * 240 7665 241 + 242 * 243 7485 244 + 245 * 246 7315 247 - 248 * 249 42374 250 XEQ 10 251 RTN 252 LBL "DEU" 253 RCL 00 254 10 255 * 256 477 257 RCL Y 258 472 259 * 260 + 261 * 262 656 263 - 264 * 265 2996 266 - 267 * 268 817 269 - 270 * 271 5160 272 + 273 * 274 2154 275 + 276 * 277 2667 278 - 279 * 280 703 281 - 282 * 283 1702 284 - 285 * 286 6137 287 - 288 * 289 1667 290 + 291 * 292 16664 293 + 294 * 295 125754 296 + 297 * 298 205070 299 + 300 STO 06 301 CLX 302 16 303 * 304 28 305 - 306 * 307 12 308 - 309 * 310 107 311 + 312 * 313 6 314 + 315 * 316 181 317 - 318 * 319 145 320 - 321 * 322 77 323 + 324 * 325 414 326 + 327 * 328 460 329 + 330 * |
331 504
332 - 333 * 334 724 335 - 336 * 337 163 338 + 339 STO 07 340 CLX 341 1323 342 * 343 216 344 + 345 * 346 4064 347 - 348 * 349 453 350 - 351 * 352 4691 353 + 354 * 355 392 356 + 357 * 358 2892 359 - 360 * 361 1172 362 - 363 * 364 1297 365 + 366 * 367 4673 368 + 369 * 370 183 371 + 372 * 373 5893 374 - 375 * 376 43811 377 XEQ 10 378 RTN 379 LBL "IXI" 380 RCL 00 381 10 382 * 383 3 384 1/X 385 X<Y? 386 GTO 02 387 CHS 388 X>Y? 389 GTO 01 390 CLX 391 306238635 392 RCL Y 393 478257091 394 * 395 + 396 * 397 CHS 398 195257313 399 + 400 * 401 96608996 402 + 403 * 404 28872448 405 - 406 * 407 10722273 408 - 409 * 410 1867391 411 + 412 * 413 508599 414 + 415 * 416 1561 417 - 418 * 419 55507 420 + 421 * 422 47102 423 + 424 * 425 110342 426 + 427 * 428 243934 429 + 430 STO 06 431 CLX 432 139053289 433 * 434 CHS 435 20697861 436 - 437 * 438 44035780 439 + 440 * |
441 4501248
442 + 443 * 444 4957851 445 - 446 * 447 191679 448 - 449 * 450 290920 451 + 452 * 453 34716 454 + 455 * 456 13263 457 + 458 * 459 2344 460 + 461 * 462 22201 463 - 464 * 465 38503 466 - 467 * 468 2524 469 + 470 STO 07 471 CLX 472 309064846 473 * 474 127118960 475 - 476 * 477 141772391 478 - 479 * 480 40745781 481 + 482 * 483 22535490 484 + 485 * 486 4603584 487 - 488 * 489 1490043 490 - 491 * 492 225660 493 + 494 * 495 47135 496 + 497 * 498 5286 499 + 500 * 501 15457 502 - 503 * 504 20001 505 - 506 * 507 43607 508 GTO 10 509 LBL 01 510 X<> L 511 * 512 1 513 + 514 2 515 / 516 873889 517 RCL Y 518 891032 519 * 520 + 521 * 522 CHS 523 10320490 524 + 525 * 526 35952667 527 + 528 * 529 55195373 530 + 531 * 532 48792377 533 + 534 * 535 26684481 536 + 537 * 538 9112980 539 + 540 * 541 1875306 542 + 543 * 544 216361 545 + 546 * 547 11782 548 + 549 * 550 62934 |
551 +
552 * 553 210510 554 + 555 STO 06 556 CLX 557 177006 558 * 559 1200296 560 + 561 * 562 3193151 563 + 564 * 565 4267221 566 + 567 * 568 2807965 569 + 570 * 571 354218 572 + 573 * 574 795371 575 - 576 * 577 616205 578 - 579 * 580 213732 581 - 582 * 583 38755 584 - 585 * 586 12471 587 - 588 * 589 16200 590 - 591 * 592 12929 593 + 594 STO 07 595 CLX 596 1177551 597 * 598 CHS 599 6909976 600 - 601 * 602 16707636 603 - 604 * 605 21063195 606 - 607 * 608 13740303 609 - 610 * 611 2577142 612 - 613 * 614 2676211 615 + 616 * 617 2261746 618 + 619 * 620 793199 621 + 622 * 623 144897 624 + 625 * 626 3989 627 + 628 * 629 4946 630 - 631 * 632 48337 633 GTO 10 634 LBL 02 635 X<> L 636 * 637 1 638 - 639 2 640 / 641 107686049 642 RCL Y 643 18010537 644 * 645 - 646 * 647 283779508 648 - 649 * 650 434288432 651 + 652 * 653 427925202 654 - 655 * 656 283981142 657 + 658 * 659 128469451 660 - |
661 *
662 39071577 663 + 664 * 665 7791740 666 - 667 * 668 977929 669 + 670 * 671 5238 672 + 673 * 674 115259 675 + 676 * 677 71080 678 - 679 STO 06 680 CLX 681 5202308 682 * 683 CHS 684 31301881 685 + 686 * 687 81975451 688 - 689 * 690 122502246 691 + 692 * 693 115286136 694 - 695 * 696 71374552 697 + 698 * 699 29521121 700 - 701 * 702 7972809 703 + 704 * 705 1286169 706 - 707 * 708 160106 709 + 710 * 711 4449 712 + 713 * 714 30593 715 - 716 * 717 12232 718 - 719 STO 07 720 CLX 721 1367220 722 * 723 CHS 724 9597021 725 + 726 * 727 29380560 728 - 729 * 730 51378306 731 + 732 * 733 56504383 734 - 735 * 736 40571961 737 + 738 * 739 19103307 740 - 741 * 742 5734975 743 + 744 * 745 1014096 746 - 747 * 748 100376 749 + 750 * 751 1793 752 + 753 * 754 16814 755 - 756 * 757 35747 758 LBL 10 759 + 760 E3 761 ST/ 06 762 ST/ 07 763 / 764 STO 08 765 RCL 07 766 RCL 06 767 END |
( 1827 bytes / SIZE 009 )
| STACK | INPUTS | OUTPUTS |
| Z | / | R |
| Y | / | B |
| X | / | L |
L = heliocentric longitude ( deg )
where B = heliocentric
latitude ( deg )
referred to the mean ecliptic and equinox of the date
R = radius vector ( AU )
Example: T = 0.1 i-e 2100/01/01 0 h TT
>>> 0.1 STO 00
• ALTJIRA
XEQ "ALT" >>>> L =
163°198
RDN B = 1°790
RDN R = 46.250 AU
• CHAOS
XEQ "CHA" or R/S >>>>
L = 197°036
RDN B = 6°850
RDN R = 45.057 AU
• DEUCALION
XEQ "DEU" or R/S >>>>
L = 341°740
RDN B = -0°351
RDN R = 42.112 AU
• IXION
XEQ "IXI" or R/S >>>>
L = 78°108
RDN B = 2°033
RDN R = 33.799 AU
3°) Salacia - Sila Nunam - Varda
Data Registers: • R00 = T in millenia since 2000/01/01 0h TT ( Register R00 is to be initialized before executing these programs )
R06 = L = heliocentric longitude ( deg )
R07 = B = heliocentric latitude ( deg )
referred to the mean ecliptic and equinox of the date
R08 = R = radius vector ( AU )
Flags: /
Subroutines: /
| 01 LBL "SAL"
02 RCL 00 03 10 04 * 05 X<0? 06 GTO 01 07 400416 08 CHS 09 RCL Y 10 162107 11 * 12 + 13 * 14 17609 15 - 16 * 17 792336 18 + 19 * 20 495009 21 - 22 * 23 555412 24 - 25 * 26 853082 27 + 28 * 29 404470 30 - 31 * 32 74919 33 + 34 * 35 18737 36 - 37 * 38 6060 39 - 40 * 41 131311 42 + 43 * 44 21837 45 - 46 GTO 02 47 LBL 01 48 599565 49 RCL Y 50 83924 51 * 52 + 53 * 54 CHS 55 1449867 56 - 57 * 58 1390144 59 - 60 * 61 142845 |
62 -
63 * 64 604143 65 + 66 * 67 265749 68 + 69 * 70 144024 71 - 72 * 73 182142 74 - 75 * 76 55519 77 - 78 * 79 11402 80 - 81 * 82 131350 83 + 84 * 85 21842 86 - 87 LBL 02 88 STO 06 89 CLX 90 2772 91 * 92 3034 93 - 94 * 95 7706 96 - 97 * 98 10891 99 + 100 * 101 4134 102 + 103 * 104 8709 105 - 106 * 107 13 108 - 109 * 110 8591 111 - 112 * 113 24507 114 + 115 * 116 3713 117 - 118 * 119 54674 120 - 121 * 122 26733 |
123 +
124 * 125 20636 126 + 127 STO 07 128 CLX 129 238 130 * 131 CHS 132 1047 133 - 134 * 135 974 136 + 137 * 138 2118 139 + 140 * 141 1152 142 - 143 * 144 713 145 - 146 * 147 895 148 - 149 * 150 437 151 + 152 * 153 3898 154 + 155 * 156 7217 157 - 158 * 159 3964 160 - 161 * 162 9748 163 + 164 * 165 43175 166 XEQ 10 167 RTN 168 LBL "SIL" 169 RCL 00 170 10 171 * 172 489 173 RCL Y 174 1666 175 * 176 + 177 * 178 5110 179 - 180 * 181 1208 182 - 183 * |
184 5677
185 + 186 * 187 1134 188 + 189 * 190 2656 191 - 192 * 193 452 194 - 195 * 196 366 197 - 198 * 199 733 200 - 201 * 202 2373 203 + 204 * 205 126824 206 + 207 * 208 117729 209 + 210 STO 06 211 CLX 212 75 213 * 214 33 215 - 216 * 217 172 218 - 219 * 220 78 221 + 222 * 223 180 224 + 225 * 226 51 227 + 228 * 229 262 230 - 231 * 232 960 233 - 234 * 235 533 236 + 237 * 238 3892 239 + 240 * 241 732 242 - 243 * 244 4875 |
245 -
246 * 247 259 248 + 249 STO 07 250 CLX 251 736 252 * 253 CHS 254 1142 255 + 256 * 257 2112 258 + 259 * 260 2777 261 - 262 * 263 2218 264 - 265 * 266 2398 267 + 268 * 269 1077 270 + 271 * 272 1029 273 - 274 * 275 453 276 - 277 * 278 806 279 + 280 * 281 630 282 + 283 * 284 857 285 - 286 * 287 43589 288 XEQ 10 289 RTN 290 LBL "VRD" 291 RCL 00 292 10 293 * 294 3257 295 RCL Y 296 619 297 * 298 + 299 * 300 127 301 - 302 * 303 11235 304 - 305 * |
306 6314
307 - 308 * 309 10929 310 + 311 * 312 13495 313 + 314 * 315 4688 316 + 317 * 318 7315 319 - 320 * 321 16265 322 - 323 * 324 12850 325 - 326 * 327 15470 328 + 329 * 330 35322 331 + 332 * 333 105085 334 + 335 * 336 239496 337 + 338 STO 06 339 CLX 340 719 341 * 342 1108 343 + 344 * 345 2324 346 - 347 * 348 5648 349 - 350 * 351 48 352 - 353 * 354 9636 355 + 356 * 357 11120 358 + 359 * 360 1021 361 + 362 * 363 12243 364 - 365 * 366 22867 |
367 -
368 * 369 22865 370 - 371 * 372 20987 373 + 374 * 375 17973 376 + 377 STO 07 378 CLX 379 1677 380 * 381 475 382 - 383 * 384 4711 385 - 386 * 387 1110 388 + 389 * 390 4508 391 + 392 * 393 1535 394 - 395 * 396 2159 397 - 398 * 399 1558 400 + 401 * 402 3610 403 + 404 * 405 4549 406 + 407 * 408 6944 409 - 410 * 411 10853 412 - 413 * 414 48908 415 LBL 10 416 + 417 E3 418 ST/ 06 419 ST/ 07 420 / 421 STO 08 422 RCL 07 423 RCL 06 424 END |
( 893 bytes / SIZE 009 )
| STACK | INPUTS | OUTPUTS |
| Z | / | R |
| Y | / | B |
| X | / | L |
L = heliocentric longitude ( deg )
where B = heliocentric
latitude ( deg )
referred to the mean ecliptic and equinox of the date
R = radius vector ( AU )
Example: T = 0.1 i-e 2100/01/01 0 h TT
>>> 0.1 STO 00
• SALACIA
XEQ "SAL" >>>> L =
94°205
RDN B = 3°233
RDN R = 45.124 AU
• SILA-NUNAM
XEQ "SIL" or R/S >>>>
L = 245°367
RDN B = -1°966
RDN R = 43.684 AU
• VARDA
XEQ "VRD" or R/S >>>>
L = 374°255 = 14°255
RDN B = -3°431
RDN R = 39.243 AU
Notes:
-The polynomials are fitted to the French ephemerides INPOP13c ( cf
reference [1] )
-Registers R01 & R02 are unused, in case you want to use them to
store the rectangular coordinates of the Sun.
-A subroutine may also be added after line 713 to calculate the geocentric
coordinates in registers R03-R04-R05
Reference: