hp41programs

INPOPHP41

INPOP-HP41 Ephemerides


Overview
 

 1°)  Heliocentric Coordinates ( Geocentric for the Moon )

  a)  Program#1
  b)  Program#2
  c)  Program#3  ( Clenshaw Algorithm )

 2°)  Data Files

  a)  2012/09/14 12h TDB - 2012/10/16 12h TDB ( 271 numbers )
  b)  2012/10/16 12h TDB - 2012/11/17 12h TDB ( 280 numbers )
 

-The programs listed in "INPOP Ephemerides for the HP-41" use a DATA file of - at least - 236 numbers  but many of them must be updated every 4 days.
-One could save all the coefficients for a time-span of 32 days, but it would require about 600 registers.

-In order to reduce the number of coefficients to be stored - and thanks to my HP-48 -  I've calculated Chebyshev coefficients that are valid over time-spans
  of 32 days to compute the rectangular equatorial coordinates of the major planets & the Moon.

-So, you will "only" have to store less than 300 numbers every 32 days. ( the exact number of coefficients may vary from month to month... )
-Moreover, the coordinates are directly heliocentric for the planets ( of course still geocentric for the Moon )

-For example, from 2012/09/14 12h TDB to 2012/10/16 12h TDB, we have:
 

ITEM      POINTER   Nb of Coeff / coordinate

JDi                 0                          1                        Initial julian day number
Mer          1   to  42                   14
Ven         43  to  69                    9
EMB       70  to  93                    8                       We could store directly the coefficients for the Earth, but it would require about 48 coefficients instead of 24
Mar         94 to 114                    7
Jup         115 to 129                   5
Sat         130 to 144                   5
Ura        145 to 156                   4
Nep       157 to 168                   4
Plu         169 to 180                   4
Moon    181 to 270                   30

TOTAL  ..............................   271 coeff.

-But for the next time-span of 32 days, a data file of 280 numbers must be created !
 

1°)   Heliocentric Coordinates ( Geocentric for the Moon )
 

     a)  Program#1
 

-First create a DATA file named "DATA" ( or another name, cf line 09 below ) and save the 271 numbers given in paragraph 2°)a)
-Many decimals may be neglected, thus saving several keystrokes.
 

Data Registers:   R00 & R33 to R65: temp                      R35 = number of Chebyshev coefficients per coordinate

                               R01 = XMercury     R04 = XVenus        R07 = XEarth        R10 = XMars        R13 = XJupiter        R16 = XSaturn
                               R02 = YMercury     R05 = YVenus        R08 = YEarth        R11 = YMars        R14 = YJupiter        R17 = YSaturn
                               R03 = ZMercury      R06 = ZVenus        R09 = ZEarth         R12 = ZMars        R15 = ZJupiter        R18 = ZSaturn

                                R19 = XUranus      R22 = XNeptune      R25 = XPluto      R28 = XMoon
                                R20 = YUranus      R23 = YNeptune      R26 = YPluto      R29 = YMoon       R31 & R32 are unused
                                R21 = ZUranus      R24 = ZNeptune       R27 = ZPluto      R30 = ZMoon
 

Flags:  /
Subroutines:  "J0" or "J1" or J2"  ( cf "Julian & Gregorian Calendars for the HP-41" )
                         "CdT" M-Code routine or focal program ( cf "Ephemerides & Chebyshev Polynomials for the HP-41 )
 
 

 01  LBL "INPHP"
 02  HR
 03  24
 04  /
 05  X<>Y          
 06  XEQ "J2"
 07  2451544.5
 08  +
 09  "DATA"
 10  0
 11  SEEKPTA
 12  CLX
 13  GETX
 14  -
 15  +
 16  16
 17  ST- Y
 18  /
 19  STO 00
 20  ACOS
 21  SIGN
 22  STO 33
 23  36.035
 24  STO 34         
 25  14
 26  XEQ 01
 27  9
 28  XEQ 01
 29  8
 30  XEQ 01
 31  7
 32  XEQ 01
 33  5
 34  XEQ 01
 35  5
 36  XEQ 01
 37  4
 38  XEQ 01
 39  4
 40  XEQ 01
 41  4
 42  XEQ 01
 43  30
 44  XEQ 01         
 45  GTO 00
 46  LBL 01
 47  STO 35
 48  3 E-3
 49  ST+ 33
 50  LBL 02
 51  RCL 34
 52  RCL 35
 53   E3
 54  /
 55  +
 56  GETRX
 57  RCL 00         
 58  CdT
 59  STO IND 33
 60  ISG 33
 61  GTO 02
 62  RTN
 63  LBL 00
 64  RCL 28
 65  82.30055986
 66  /
 67  ST- 07
 68  RCL 29         
 69  LASTX
 70  /
 71  ST- 08
 72  RCL 30
 73  LASTX
 74  /
 75  ST- 09
 76  END

 
    ( 165 bytes / SIZE 066 )
 
 

      STACK        INPUTS      OUTPUTS
           Y  YYYY.MNDD             /
           X      HH.MNSS             /

   where the time is expressed in TDB                                           ---Execution time = 50s---

Example:     2012/09/28  16h41m37s   TDB

    2012.0829    ENTER^
        16.4137    XEQ "INPHP"   returns

  •  Mercury

      R01 = -0.286284468  AU
      R02 = -0.322647018  AU
      R03 = -0.142670773  AU

  •  Venus

      R04 =  0.150015436  AU
      R05 =  0.646073594  AU
      R06 =  0.281193802  AU

  •  Earth

      R07 =  0.996686440  AU
      R08 =  0.092563987  AU
      R09 =  0.040123397  AU

  •  Mars

     R10 = -0.141329564  AU
     R11 = -1.328790361  AU
     R12 = -0.605660773  AU

  •  Jupiter

    R13 =  2.100486858  AU
    R14 =  4.224182200  AU
    R15 =  1.759464866  AU

  •  Saturn

    R16 = -8.336965937  AU
    R17 = -4.820164526  AU
    R18 = -1.631899460  AU

  •  Uranus

    R19 =  19.93790927  AU
    R20 =  2.138309326  AU
    R21 =  0.654582854  AU

  •  Neptune

    R22 =  26.42461537  AU
    R23 = -12.88997194  AU
    R24 = -5.933673859  AU

  •  Pluto

    R25 =  4.811481388  AU
    R26 = -30.05517705  AU
    R27 = -10.82812720  AU

  •  Moon

    R28 =  0.002554665534  AU
    R29 = -0.000545255583  AU
    R30 = -0.000000874463  AU

Notes:

-The execution time depends on the M-Code version of CdT that you are using.
-The M-Code CdT may of course be replaced by a focal program "CdT".

-Lines 64 to 75 replace the coordinates of the Earth-Moon-Barycenter by the coordinates of the Earth in registers R07-R08-R09.
-R31 & R32 are unused in case you want to calculate the nutations...

-This program must be modified according to the number of coefficients / coordinate:
-For instance, with the data file of paragraph 2°)b), replace line 25 by 17 ( instead of 14 )
 

     b)  Program#2
 

-This variant uses the same data file of 271 numbers.
 

Data Registers:   R00 & R33-R34-R35: temp                  R35 = -2 + number of Chebyshev coefficients per coordinate

                               R01 = XMercury     R04 = XVenus        R07 = XEarth        R10 = XMars        R13 = XJupiter        R16 = XSaturn
                               R02 = YMercury     R05 = YVenus        R08 = YEarth        R11 = YMars        R14 = YJupiter        R17 = YSaturn
                               R03 = ZMercury      R06 = ZVenus        R09 = ZEarth         R12 = ZMars        R15 = ZJupiter        R18 = ZSaturn

                                R19 = XUranus      R22 = XNeptune      R25 = XPluto      R28 = XMoon
                                R20 = YUranus      R23 = YNeptune      R26 = YPluto      R29 = YMoon       R31 & R32 are unused
                                R21 = ZUranus      R24 = ZNeptune       R27 = ZPluto      R30 = ZMoon
 

Flags:  /
Subroutine:  "J0" or "J1" or J2"  ( cf "Julian & Gregorian Calendars for the HP-41" )
 
 

 01  LBL "INPHP"
 02  HR
 03  24
 04  /
 05  X<>Y
 06  XEQ "J2"
 07  2451544.5
 08  +
 09  "DATA"
 10  0
 11  SEEKPTA
 12  CLX
 13  GETX
 14  -
 15  +
 16  16
 17  ST- Y
 18  /
 19  STO 00
 20  ACOS
 21  SIGN
 22  STO 33
 23  12
 24  XEQ 01
 25  7
 26  XEQ 01
 27  6
 28  XEQ 01         
 29  5
 30  XEQ 01
 31  3
 32  XEQ 01
 33  3
 34  XEQ 01
 35  2
 36  XEQ 01
 37  2
 38  XEQ 01
 39  2
 40  XEQ 01
 41  28
 42  XEQ 01         
 43  GTO 00
 44  LBL 01
 45  STO 35
 46  3 E-3
 47  ST+ 33
 48  LBL 02
 49  GETX
 50  RCL 00
 51  GETX
 52  *
 53  +
 54  STO IND 33
 55  RCL 35
 56  STO 34
 57  SIGN
 58  RCL 00
 59  LBL 03
 60  RCL 00         
 61  ST+ X
 62  RCL Y
 63  *
 64  R^
 65  -
 66  GETX
 67  RCL Y
 68  *
 69  ST+ IND 33
 70  RDN
 71  DSE 34
 72  GTO 03
 73  ISG 33
 74  GTO 02
 75  RTN
 76  LBL 00
 77  RCL 28         
 78  82.30055986
 79  /
 80  ST- 07
 81  RCL 29
 82  LASTX
 83  /
 84  ST- 08
 85  RCL 30
 86  LASTX
 87  /
 88  ST- 09
 89  END

 
   ( 182 bytes / SIZE 036 )
 
 

      STACK        INPUTS      OUTPUTS
           Y  YYYY.MNDD             /
           X      HH.MNSS             /

   where the time is expressed in TDB

Example:    Same example ->  same results, with perhaps a small difference in the last decimal              ---Execution time = 2m34s---

Notes:

-This program must be modified according to the number of coefficients / coordinate:
-For instance, with the data file of paragraph 2°)b), replace line 23 by 15 ( instead of 12 )
 

     c)  Program#3  ( Clenshaw Algorithm )
 

-This variant uses a similar data file of 271 numbers:
-But instead of storing the coefficients in the order given in paragraph 2°)

  JDi a0 ........... an  b0 ........... bn c0 ........... cn , a'0 ................. a'n' b'0 ............... b'n'  c'0 ............... c'n' , ..............

  the order must be reversed for each coordinate, i-e:

  JDi an .............. a0  bn ............... b0 cn ............... c0 , a'n' ................. a'0 b'n' ............... b'0  c'n' ............... c'0 , ..............
        X-coordinate     Y-coordinate      Z-coordinate    ,    X-coordinate       Y-coordinate          Z-coordinate   , ...............
 

---------------------------------------------------DATA FILE prgm 1,2 <> DATAFILE prgm 3-------------------------------------------------------
 

-If you have already stored the 271 numbers in the 1st way, use the following routine to get the proper data file for the 3rd version of "INPHP"
 

Data Registers:     R00 = number of coefficients / coordinate  ,  R01 thru R04: temp
Flags: /
Subroutines: /
 
 

 01  LBL "DSWAP"
 02  1
 03  STO 03           
 04  14
 05  STO 00
 06  XEQ 00
 07  9
 08  STO 00
 09  XEQ 00
 10  8
 11  STO 00
 12  XEQ 00
 13  7
 14  STO 00
 15  XEQ 00
 16  5
 17  STO 00           
 18  XEQ 00
 19  XEQ 00
 20  4
 21  STO 00
 22  XEQ 00
 23  XEQ 00
 24  XEQ 00
 25  30
 26  STO 00
 27  LBL 00
 28  3
 29  STO 04           
 30  LBL 01
 31  RCL 03
 32  STO 01
 33  RCL 00
 34  +
 35  STO 02
 36  DSE 02
 37  LBL 02
 38  RCL 02
 39  SEEKPTA
 40  GETX
 41  RCL 01           
 42  SEEKPT
 43  GETX
 44  RDN
 45  SEEKPT
 46  RDN
 47  SAVEX
 48  RDN
 49  SEEKPT
 50  X<>Y
 51  SAVEX
 52  1
 53  ST+ 01
 54  ST- 02
 55  RCL 01
 56  RCL 02
 57  X>Y?
 58  GTO 02
 59  RCL 00           
 60  ST+ 03
 61  DSE 04
 62  GTO 01
 63  END

 
    ( 108 bytes / SIZE 005 )
 
 

      STACK        INPUTS      OUTPUTS
           X             /             /
 
         ---Execution time = 3m35s---

-Simply place the name of the data file in the alpha register and XEQ "DSWAP"
-If the routine is executed twice, the original order is restored.

-This routine must also be modified according to the number of coefficients / coordinate:
-With the data file of §2°)b) replace line 04 by 17 ( instead of 14 )
 

-------------------------------------------------------------------------------------------------------------------------------------------------
 

Data Registers:   R00 & R33-R34-R35: temp                  R35 = -1 + number of Chebyshev coefficients per coordinate

                               R01 = XMercury     R04 = XVenus        R07 = XEarth        R10 = XMars        R13 = XJupiter        R16 = XSaturn
                               R02 = YMercury     R05 = YVenus        R08 = YEarth        R11 = YMars        R14 = YJupiter        R17 = YSaturn
                               R03 = ZMercury      R06 = ZVenus        R09 = ZEarth         R12 = ZMars        R15 = ZJupiter        R18 = ZSaturn

                                R19 = XUranus      R22 = XNeptune      R25 = XPluto      R28 = XMoon
                                R20 = YUranus      R23 = YNeptune      R26 = YPluto      R29 = YMoon       R31 & R32 are unused
                                R21 = ZUranus      R24 = ZNeptune       R27 = ZPluto      R30 = ZMoon
 

Flags:  /
Subroutine:  "J0" or "J1" or J2"  ( cf "Julian & Gregorian Calendars for the HP-41" )
 
 

 01  LBL "INPHP"
 02  HR
 03  24
 04  /
 05  X<>Y
 06  XEQ "J2"
 07  2451544.5
 08  +
 09  "DATA"
 10  0
 11  SEEKPTA
 12  CLX
 13  GETX
 14  -
 15  +
 16  16
 17  ST- Y
 18  /
 19  STO 00
 20  ACOS
 21  SIGN
 22  STO 33         
 23  13
 24  XEQ 01
 25  8
 26  XEQ 01
 27  7
 28  XEQ 01
 29  6
 30  XEQ 01
 31  4
 32  XEQ 01
 33  4
 34  XEQ 01
 35  3
 36  XEQ 01
 37  3
 38  XEQ 01
 39  3
 40  XEQ 01
 41  29
 42  XEQ 01         
 43  GTO 00
 44  LBL 01
 45  STO 35
 46  3 E-3
 47  ST+ 33
 48  LBL 02
 49  RCL 35
 50  STO 34
 51  CLST
 52  LBL 03
 53  RCL X
 54  RCL 00
 55  ST+ X
 56  *
 57  R^
 58  -
 59  GETX
 60  +
 61  DSE 34
 62  GTO 03
 63  RCL 00         
 64  *
 65  X<>Y
 66  -
 67  GETX
 68  +
 69  STO IND 33
 70  ISG 33
 71  GTO 02
 72  RTN
 73  LBL 00
 74  RCL 28
 75  82.30055986
 76  /
 77  ST- 07
 78  RCL 29         
 89  LASTX
 80  /
 81  ST- 08
 82  RCL 30
 83  LASTX
 84  /
 85  ST- 09
 86  END

 
   ( 172 bytes / SIZE 036 )
 
 

      STACK        INPUTS      OUTPUTS
           Y  YYYY.MNDD             /
           X      HH.MNSS             /

   where the time is expressed in TDB

Example:    Same example ->  same results                   ---Execution time = 2m14s---

Notes:

-This program is faster than the 2nd version and perhaps slightly more accurate.

-Like the first 2 versions, this program must also be modified according to the number of coefficients / coordinate:
-With the data file of paragraph 2°)b), replace line 23 by 16 ( instead of 13 )
 

2°)  Data Files
 

     a)   2012/09/14 12h TDB - 2012/10/16 12h TDB  ( 271 numbers )
 

ITEM      POINTER   Nb of Coeff / coordinate

JDi                 0                          1                        Initial julian day number
Mer          1   to  42                   14
Ven         43  to  69                    9
EMB       70  to  93                    8
Mar         94 to 114                    7
Jup         115 to 129                   5
Sat         130 to 144                   5
Ura        145 to 156                   4
Nep       157 to 168                   4
Plu         169 to 180                   4
Moon    181 to 270                   30

TOTAL  ..............................   271 coeff.
 

{   2456185                              Initial Julian Day number

{ -.204953673379                    X-Mercury
.248911474092
5.09637897868E-2
-.011310300414
1.54204585111E-4
-3.37874767222E-5
-1.45954910556E-5
3.00270594444E-6
-6.17101166667E-7
.000000096607
-9.91744444444E-9
5.32277777778E-10
2.10944444444E-10
-7.23888888889E-11}

{ -.280551606731                   Y-Mercury
-.180444555271
6.47351168172E-2
4.39374019872E-3
-8.58037940772E-4
1.20323983383E-4
-2.17892513222E-5
1.44599713889E-6
-1.10631644444E-7
-4.23739E-8
1.17876277778E-8
-2.78907777778E-9
4.66183333333E-10
-6.03944444444E-11}

{ -.128616380242                    Z-Mercury
-.122196906678
2.92965100716E-2
3.51968875593E-3
-4.74336357978E-4
6.77779250944E-5
-1.01261939778E-5
4.6111065E-7
4.8826E-9
-3.26515277778E-8
7.32488888889E-9
-1.54494444444E-9
2.27355555556E-10
-2.485E-11 }

{ .108325226011                    X-Venus
-.312599628174
-5.66317315161E-3
2.66730011733E-3
2.90081926667E-5
-6.97035183333E-6
-8.51121111111E-8
9.56066666667E-9
2.18777777778E-10 }

{ .618590535897                    Y-Venus
3.68126690602E-2
-3.24527622914E-2
-3.52397028556E-4
1.391237145E-4
.0000013155
-2.51675444444E-7
-3.75738888889E-9
2.89611111111E-10 }

{ .271466549159                    Z-Venus
3.63443356599E-2
-1.42430501526E-2
-3.27339968278E-4
6.07600668889E-5
1.03296572222E-6
-1.07850611111E-7
-2.29572222222E-9
1.16277777778E-10}

{ .973971406228                    X-Earth-Moon-Barycenter
-4.04065745627E-2
-1.86075136886E-2
8.54362785556E-5
2.95316265556E-5
8.09505E-8
-1.83922222222E-8
-3.31944444444E-10 }

{ .118465123483                    Y-Earth-Moon-Barycenter
.247063585786
-2.27410652333E-3
-7.85713013667E-4
8.98857222222E-7
7.44057333333E-7
6.20005555556E-9
-2.28111111111E-10}

{ 5.13548869036E-2               Z-Earth-Moon-Barycenter
.107106402636
-9.85834193261E-4
-3.40620676378E-4
3.89588722222E-7
3.22575333333E-7
2.69700555556E-9
-1.04422222222E-10}

{ -.114535085777                   X-Mars
.230991319782
6.83180879278E-4
-2.27626956167E-4
-2.68045894444E-6
6.34121666667E-8
2.66233333333E-9}

{ -1.32058315343                   Y-Mars
3.86562876278E-3
7.98734225111E-3
4.96853783333E-5
-4.06667055556E-6
-.00000008465
2.79444444444E-10}

{ -.602619988339                   Z-Mars
-4.46528391494E-3
3.64512153467E-3
2.89366763333E-5
-1.7928965E-6
-4.05360555556E-8
5.55E-11}

{ 2.08761811727                    X-Jupiter
-.111420298391
-3.10005737778E-4
2.90335944444E-6
2.42333333333E-9 }

{ 4.22924902075                    Y-Jupiter
5.03519610172E-2
-6.28010999444E-4
-9.61771666667E-7
9.41055555556E-9}

{ 1.76194995501                    Z-Jupiter
2.42949726044E-2
-2.61638677778E-4
-4.82973888889E-7
3.96222222222E-9 }

{ -8.33211586272                    X-Saturn
4.15442985772E-2
1.69074342778E-4
-1.70292222222E-7
-1.77777777778E-11}

{-4.82795777714                     Y-Saturn
-6.99452797944E-2
9.75231061111E-5
2.12277777778E-7
-1.33333333333E-11 }

{ -1.63532739461                    Z-Saturn
-.03068041031
3.29813905556E-5
9.50472222222E-8
-1.22222222222E-11 }

{19.9370129671                      X-Uranus
-7.54025537778E-3
-4.71276055556E-5
9.57222222222E-9}

{ 2.1444609077                       Y-Uranus
5.45953274944E-2
-5.67050444444E-6
-2.48322222222E-8}

{ .657289852944                      Z-Uranus
2.40187223943E-2
-1.81140388889E-6
-1.11069444444E-8}

{ 26.4272275671                      X-Neptune
2.33256858333E-2
-1.88808111111E-5
4.69444444444E-9}

{ -12.8852847947                      Y-Neptune
4.14884448611E-2
8.41400555556E-6
-8.15555555556E-9}

{ -5.93182035547                      Z-Neptune
1.64014222356E-2
3.89795666667E-6
-3.58833333333E-9}

{ 4.81717593956                       X-Pluto
5.05199940378E-2
-3.01985277778E-6
3.19111111111E-9 }

{ -30.05477167                         Y-Pluto
3.45400768333E-3
1.62541111111E-5
-6.13888888889E-9}

{ -10.8297191038                     Z-Pluto
-1.41656977278E-2
5.81723333333E-6
-1.22222222222E-9}

{-8.93465539411E-4                 X-Moon
-4.50104139444E-5
-2.2033887139E-3
2.5456011806E-4
1.15996437877E-3
-8.90660930339E-5
-1.25225428817E-4
2.76446729256E-5
-3.08393448E-6
-5.487007755E-6
9.67418168333E-7
3.42772954444E-7
1.79237124444E-7
5.58171222222E-9
-4.14951766667E-8
2.21858033333E-8
-2.07796777778E-9
-7.50548888889E-9
1.65045888889E-9
5.57391111111E-10
-2.52761666667E-10
7.24516666667E-11
1.39388888889E-11
8.01222222222E-12
-1.31666666667E-12
-8.23833333333E-12
2.11288888889E-12
8.59222222222E-13
-8.22777777778E-13
1.74111111111E-13 }

{ 5.55593140675E-5                 Y-Moon
3.17069326223E-4
-1.9441260971E-4
-1.98364406594E-3
1.17112196777E-4
4.32993209901E-4
-4.36863172056E-5
-1.84480303536E-5
1.27760400369E-5
-2.78454386833E-6
-1.63355389444E-6
-1.64266044444E-8
7.68652333333E-9
1.05490719444E-7
-1.48158254444E-8
-7.83113561111E-9
1.45043627222E-8
-2.43744338889E-9
-2.51114505556E-9
6.59055833333E-10
-1.46697777778E-11
-6.80536111111E-11
2.03123333333E-11
1.44433333333E-12
1.21227777778E-11
-1.79411111111E-12
-3.2855E-12
1.37111111111E-12
-1.06111111111E-14
-3.40666666667E-13}

{-5.15808009181E-5                 Z-Moon
1.15197052742E-4
-2.48193030574E-4
-7.26636090733E-4
1.38955885352E-4
1.55876856022E-4
-2.69434068284E-5
-4.69827843311E-6
4.58707098261E-6
-1.50020071217E-6
-5.34885814667E-7
2.28763673333E-8
1.71612764444E-8
4.02080659444E-8
-9.00195394444E-9
-1.17266594444E-9
5.30698333333E-9
-1.53185388889E-9
-8.11871055556E-10
2.95426333333E-10
-2.63738333333E-11
-1.98141111111E-11
8.83577777778E-12
1.15227777778E-12
4.46244444444E-12
-1.34594444444E-12
-1.06766666667E-12
5.89888888889E-13
-7.08666666667E-14
-1.15044444444E-13 }
}
 

     b)   2012/10/16 12h TDB - 2012/11/17 12h TDB  ( 280 numbers )
 

ITEM      POINTER   Nb of Coeff / coordinate

JDi                 0                          1                        Initial julian day number
Mer          1   to  51                   17
Ven         52  to  78                    9
EMB       79  to  102                  8
Mar        103 to 123                   7
Jup         124 to 138                   5
Sat         139 to 153                   5
Ura        154 to 165                   4
Nep       166 to 177                   4
Plu         178 to 189                   4
Moon    190 to 279                   30

TOTAL  ..............................   280 coeff.
 

{   2456217                               Initial Julian Day

 { .238401539586                     MERCURY
.062684303577
-.10766271724
-1.50616428904E-2
1.2350570862E-3
8.0542906096E-4
1.94596028525E-4
1.7753117685E-5
-5.335295645E-6
-2.87852997E-6
-6.4692435E-7
-3.553609E-8
3.1323995E-8
1.364001E-8
2.69218E-9
-5.535E-12
-1.98395E-10 }
{ -.106387413321
.329481319329
2.93322543032E-2
-1.72088117324E-2
-3.66540634375E-3
-3.261774041E-4
7.45087822E-5
4.154065815E-5
9.2582413E-6
5.896582E-7
-3.8773605E-7
-1.738078E-7
-3.505405E-8
-4.421E-10
2.32195E-9
8.7405E-10
1.497E-10 }
{ -8.15472134125E-2
.169504006302
2.68310151811E-2
-7.6310795048E-3
-2.0860426415E-3
-2.5774330575E-4
1.962594955E-5
2.034966955E-5
5.49873915E-6
6.134249E-7
-1.4005015E-7
-8.916075E-8
-2.19728E-8
-1.65045E-9
9.612E-10
4.6765E-10
1.008E-10 }

{ -.462391447423                      VENUS
-.232808613313
2.44408919858E-2
2.0051791906E-3
-1.068301295E-4
-5.4393461E-6
2.0821295E-7
.000000008524
-3.0245E-10 }
{ .445833250325
-.202105694458
-2.35661797563E-2
1.7419255234E-3
1.028586645E-4
-4.7381082E-6
-1.9936955E-7
.000000007496
2.8585E-10 }
{ .229853059288
-7.62009589495E-2
-1.21497078071E-2
6.5685395895E-4
5.303917775E-5
-1.78761795E-6
-1.028773E-7
2.83375E-9
1.4785E-10 }

{ .75319024866                      EARTH-MOON-BARYCENTER
-.177178246958
-1.47782761228E-2
.000544452248
2.616517335E-5
-4.3949225E-7
-2.37331E-8
7.175E-11 }
{ .565483577985
.192568027423
-1.11083255258E-2
-6.4767312895E-4
1.640167725E-5
7.550776E-7
-6.67355E-9
-4.9985E-10 }
{ .245145461574
8.34818393905E-2
-4.8156319045E-3
-2.8077881615E-4
7.11033495E-6
3.273363E-7
-2.90365E-9
-2.2265E-10 }

{ .343826757188                           MARS
.224921313164
-2.25215546895E-3
-2.567577503E-4
-7.4540665E-7
1.299867E-7
2.62065E-9 }
{ -1.2479274198
.068919825375
.008140136183
-2.79472715E-5
-5.5437085E-6
-5.51995E-8
.000000002288 }
{ -.581673721415
.025537193752
3.79447595625E-3
-5.8847065E-6
-2.5226661E-6
-2.884325E-8
9.783E-10 }

{ 1.86240820869                            JUPITER
-.113760412962
-2.749709735E-4
2.9333835E-6
.000000001651 }
{ 4.32489411788
4.52842713985E-2
-6.386496175E-4
-8.126455E-7
8.8765E-9 }
{ 1.80842920338
.02217976467
-.000267052438
-4.199115E-7
.000000003701 }

{ -8.2476811734                          SATURN
.042888672178
1.670111665E-4
-1.744255E-7
-1.795E-10 }
{ -4.9670600847
-.069154905806
.00010006853
2.106955E-7
-5.405E-10 }
{ -1.69642074985
-.030411993865
3.41222225E-5
.000000094484
-2.765E-10}

{ 19.9215558102                           URANUS
-.00791680424
-.00004701164
.00000000887 }
{ 2.25360529954
5.45488334895E-2
-.000005948509
-.000000022757 }
{ .70531240308
.024003724022
-1.9363108E-6
-1.03199E-8 }

{ 26.4737280785                          NEPTUNE
.02317487305
-.00001882439
.00000000393  }
{ -12.8022408585
.041555427025
.00000833602
-.0000000061}
{ -5.89898644423
1.64324606425E-2
.00000386352
-2.7905E-9 }

{ 4.91819189843                          PLUTO
5.04959976095E-2
-.000002981282
2.3845E-9 }
{ -30.0477338065
.003583812775
.00001620138
-.000000004005 }
{ -10.8580039886
-.01411919183
.00000581112
-.000000000445 }

{ -2.31665937163E-4                   MOON
-2.99454320647E-4
-7.70915082265E-4
1.98552658243E-3
4.16813684462E-4
-4.11074828487E-4
-5.08374985605E-5
3.6401878555E-6
-2.283171433E-6
5.097711403E-6
2.03690002E-6
4.238860075E-7
-3.05551458E-7
-2.53251568E-7
-3.13805395E-8
-3.408995E-10
1.226585E-8
1.0443432E-8
1.8680585E-9
-7.95988E-10
-9.698705E-10
-3.813555E-10
-1.43545E-11
5.76465E-11
5.18915E-11
1.71385E-11
-2.684E-12
-4.91855E-12
-2.62895E-12
-4.7725E-13 }
{ -7.22479516378E-4
1.4519869895E-4
-1.96134725696E-3
-7.677268893E-4
1.01515345107E-3
1.74995106323E-4
-9.07220116905E-5
-8.317441278E-6
-9.092517747E-6
-3.155784846E-6
8.88559542E-7
8.719996205E-7
4.7461499E-7
-3.0439228E-8
-6.84425895E-8
-2.6903416E-8
-1.4523717E-8
1.3418475E-9
3.9006351E-9
1.77929535E-9
3.157524E-10
-2.5788865E-10
-1.882112E-10
-7.28238E-11
-5.6243E-12
1.854105E-11
1.12219E-11
2.70705E-12
-7.826E-13
-1.1879E-12 }
{ -2.891536089E-4
3.06243962057E-5
-7.9462233269E-4
-1.24832986426E-4
4.13596660337E-4
3.15357897567E-5
-3.8134884062E-5
-2.76037656295E-6
-3.57858441825E-6
-7.619957368E-7
4.9923675655E-7
3.597559531E-7
1.521372023E-7
-3.215192195E-8
-2.81222721E-8
-1.004612265E-8
-4.4177363E-9
1.356442E-9
1.6089143E-9
5.9829705E-10
3.836095E-11
-1.275091E-10
-7.13827E-11
-2.24338E-11
2.15005E-12
8.32315E-12
3.96605E-12
6.0665E-13
-5.0745E-13
-4.823E-13 }
}
 

References:

[1]  http://www.imcce.fr/inpop
[2]  Press , Vetterling , Teukolsky , Flannery - "Numerical Recipes in C++" - Cambridge University Press - ISBN  0-521-75033-4