% Positions U, velocities V and masses M of the nine planets,
% the Sun and the Moon in the following order:
%
%	Mercury
%	Venus
%	Earth
%	Mars
%	Jupiter
%	Saturn
%	Uranus
%	Neptune
%	Pluto
%	Sun
%	Moon
%
% at 2000-01-01 00.00 Universal Time (UT1, approx GMT)
%
% Units: Length: AU, 1 Astronomical Unit
%        Time:   a,  1 year approximately
%        Mass:   M,  1 solar mass
%
% Comment: a would be the length of one year if earths orbit were
%          a circle...
%
% Positions obtained from U.S. Naval Observatory, www.usno.navy.mil.
%
% Velocities obtained by fitting a 10:th degree polynomial to the 
% positions of December 1999, or more exactly the 32 positions
% 1999-12-01, 00:00 - 1999-12-02, 00:00 - ... - 2000-01-01, 00:00.
%
% Masses obtained from "Fundamental Astronomy", 3:rd edition, by
% Karttunen, Kröger, Oja, Poutanen, Donner.
%
% Anders Logg, logg@math.chalmers.se
% 2000-11-14

% Positions, X Y Z, barycentric coordinates
U = [
  -0.147853935  -0.400627944  -0.198916163
  -0.725771746  -0.039677000   0.027897127
  -0.175679599   0.886201933   0.384435698
   1.383219717  -0.008134314  -0.041033184
   3.996313003   2.731004338   1.073280866
   6.401404019   6.170259699   2.273032684
  14.423408013 -12.510136707  -5.683124574
  16.803677095 -22.983473914  -9.825609566
  -9.884656563 -27.981265594  -5.753969974
  -0.007141917  -0.002638933  -0.000919462
  -0.177802714   0.884620944   0.384016593 ];

% Velocities, X Y Z, barycentric coordinates
V = [
   7.733816715  -2.014137426  -1.877564183
   0.189682646  -6.762413869  -3.054194695
  -6.292645274  -1.010423954  -0.438086386
   0.275092348   5.042903370   2.305658434
  -1.664796930   2.146870503   0.960782651
  -1.565320566   1.286649577   0.598747577
   0.980209400   0.896663122   0.378850106
   0.944045755   0.606863295   0.224889959
   1.108139341  -0.414389073  -0.463196118
   0.001962209  -0.002469700  -0.001108260
  -6.164023246  -1.164502534  -0.506131880 ];

% Comment: Earths velocity is 2*pi... :)

% Masses
M = [
   1.0 / 6023600
   1.0 / 408523.5
   1.0 / 328900.5
   1.0 / 3098710
   1.0 / 1047.355
   1.0 / 3498.5
   1.0 / 22869
   1.0 / 19314
   1.0 / 150000000
   1.0 / 1.0
   1.0 / (328900.5*81.30) ];

% Some constants
Ms = 1.989e30;      % kg, mass of Sun
AU = 1.49597870e11; % m, Astronomical Unit
a  = 365.24;        % days
G  = 4*pi^2;        % AU^3 / (M a^2), gravitational constant

% written by Anders Logg, November 2000. 

