% Data for outer solar system simulation
% Units: mass in solar mass, time in Earth day, length in astronomical unit
% Source: J.A.Barth. Ahnerts Kalender fur Sternfreunde 1994. Verlag

g=2.95912208286e-4;         % the gravitational constant

m=zeros(6,1);                 % masses
m(1)=1.00000597682;         % Sun + inner planets
m(2)=9.54786104043e-4;      % Jupiter
m(3)=2.85583733151e-4;      % Saturn
m(4)=4.37273164546e-5;      % Uranus
m(5)=5.17759138449e-5;      % Neptune
m(6)=7.692307692307693e-9;  % Pluto

x=zeros(36,1);               

% Initial positions at 1994-09-05, UTC0:00
x(4)=-3.5023653;          % Jupiter
x(5)=-3.8169847;
x(6)=-1.5507963;

x(7)= 9.0755314;          % Saturn
x(8)=-3.0458353;
x(9)=-1.6483708;

x(10)=  8.3101420;          % Uranus
x(11)=-16.2901086;
x(12)= -7.2521278;

x(13)= 11.4707666;          % Neptune
x(14)=-25.7294829;
x(15)=-10.8169456;

x(16)=-15.5387357;          % Pluto
x(17)=-25.2225594;
x(18)= -3.1902382;

% Initial velocities at 1994-09-05, UTC0:00
x(22)= 0.00565429;          % Jupiter
x(23)=-0.00412490;
x(24)=-0.00190589;

x(25)= 0.00168318;          % Saturn
x(26)= 0.00483525;
x(27)= 0.00192462;

x(28)= 0.00354178;          % Uranus
x(29)= 0.00137102;
x(30)= 0.00055029;

x(31)= 0.00288930;          % Neptune
x(32)= 0.00114527;
x(33)= 0.00039677;

x(34)= 0.00276725;          % Pluto
x(35)=-0.00170702;
x(36)=-0.00136504;