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 = X_Mercury     R04 = X_Venus        R07 = X_Earth        R10 = X_Mars        R13 = X_Jupiter        R16 = X_Saturn
                               R02 = Y_Mercury     R05 = Y_Venus        R08 = Y_Earth        R11 = Y_Mars        R14 = Y_Jupiter        R17 = Y_Saturn
                               R03 = Z_Mercury      R06 = Z_Venus        R09 = Z_Earth         R12 = Z_Mars        R15 = Z_Jupiter        R18 = Z_Saturn

                                R19 = X_Uranus      R22 = X_Neptune      R25 = X_Pluto      R28 = X_Moon
                                R20 = Y_Uranus      R23 = Y_Neptune      R26 = Y_Pluto      R29 = Y_Moon       R31 & R32 are unused
                                R21 = Z_Uranus      R24 = Z_Neptune       R27 = Z_Pluto      R30 = Z_Moon
 

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 = X_Mercury     R04 = X_Venus        R07 = X_Earth        R10 = X_Mars        R13 = X_Jupiter        R16 = X_Saturn
                               R02 = Y_Mercury     R05 = Y_Venus        R08 = Y_Earth        R11 = Y_Mars        R14 = Y_Jupiter        R17 = Y_Saturn
                               R03 = Z_Mercury      R06 = Z_Venus        R09 = Z_Earth         R12 = Z_Mars        R15 = Z_Jupiter        R18 = Z_Saturn

                                R19 = X_Uranus      R22 = X_Neptune      R25 = X_Pluto      R28 = X_Moon
                                R20 = Y_Uranus      R23 = Y_Neptune      R26 = Y_Pluto      R29 = Y_Moon       R31 & R32 are unused
                                R21 = Z_Uranus      R24 = Z_Neptune       R27 = Z_Pluto      R30 = Z_Moon
 

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°)

  JD_i a₀ ........... a_n  b₀ ........... b_n c₀ ........... c_n , a'₀ ................. a'_n' b'₀ ............... b'_n'  c'₀ ............... c'_n' , ..............

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

  JD_i a_n .............. a₀  b_n ............... b₀ c_n ............... c₀ , a'_n' ................. a'₀ b'_n' ............... b'₀  c'_n' ............... c'₀ , ..............
        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 = X_Mercury     R04 = X_Venus        R07 = X_Earth        R10 = X_Mars        R13 = X_Jupiter        R16 = X_Saturn
                               R02 = Y_Mercury     R05 = Y_Venus        R08 = Y_Earth        R11 = Y_Mars        R14 = Y_Jupiter        R17 = Y_Saturn
                               R03 = Z_Mercury      R06 = Z_Venus        R09 = Z_Earth         R12 = Z_Mars        R15 = Z_Jupiter        R18 = Z_Saturn

                                R19 = X_Uranus      R22 = X_Neptune      R25 = X_Pluto      R28 = X_Moon
                                R20 = Y_Uranus      R23 = Y_Neptune      R26 = Y_Pluto      R29 = Y_Moon       R31 & R32 are unused
                                R21 = Z_Uranus      R24 = Z_Neptune       R27 = Z_Pluto      R30 = Z_Moon
 

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