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: /
 
 
 

 
          ( 1465 bytes / SIZE 009 )
 
 

                    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: /
 
 

 
      ( 1827 bytes / SIZE 009 )
 
 

                    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: /
 
 

 
       ( 893 bytes / SIZE 009 )
 
 

                    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:

[1]   http://www.imcce.fr