from visual import * # initial positions in au. # on June 28.0, 1969 (Julian day 2440400.50) r = array( [ 4.504795855674601e-03, # Sun 7.732544746890742e-04, 2.685039985573271e-04, 3.617650084026384e-01, # Mercury -9.077580081092690e-02, -8.571250013597834e-02, 6.127542336221157e-01, # Venus -3.483591895730089e-01, -1.952758215816601e-01, 1.205197129002194e-01, # Earth+moon -9.258323035784091e-01, -4.015377611192672e-01, -1.101837696063662e-01, # Mars -1.327593278911103e+00, -6.058866834707231e-01, -5.379703845550703e+00, # Jupiter -8.304767428789135e-01, -2.248295252614748e-01, 7.894394222587950e+00, # Saturn 4.596483181197301e+00, 1.558697670343092e+00, -1.826538634212984e+01, # Neptune -1.161959795438843e+00, -2.501080003472951e-01, -1.605499953831724e+01, # Uranus -2.394216805154041e+01, -9.400159225576589e+00, -3.048349217668479e+01, # Pluto -8.724432817555443e-01, 8.911620591248875e+00 ]) # initial velocities in au/day v = array( [ -3.524457445683394e-07, # SUN 5.177637780672221e-06, 2.229113252400401e-06, 3.367492758813184e-03, # Mercury 2.489452047536652e-02, 1.294630071041128e-02, 1.095206748972985e-02, # Venus 1.561768454932838e-02, 6.331105550175031e-03, 1.681126759384403e-02, # Earth+moon 1.748308899617618e-03, 7.582041930175932e-04, 1.448165239083299e-02, # Mars 2.424628123494023e-04, -2.815196436092114e-04, 1.092015004926191e-03, # Jupiter -6.518116265389094e-03, -2.820782997470719e-03, -3.217557586964646e-03, # Saturn 4.335809727271375e-03, 1.928645932522133e-03, 2.211920155513113e-04, # Neptune -3.762477278855801e-03, -1.651014777644869e-03, 2.642773506889557e-03, # Uranus -1.498309226803948e-03, -6.790394433396178e-04, 3.221893230538327e-04, # Pluto -3.143582327761257e-03, -1.077956399281471e-03 ] ) # Gm for bodies. Gm = array([ 2.959122082855911e-004, # SUN 4.912547451450812e-011, # Mercury 7.243456209632766e-010, # Venus 8.997011658557308e-010, # Earth 9.549528942224058e-011, # Mars 2.825342103445926e-007, # Jupiter 8.459468504830660e-008, # Saturn 1.288816238138035e-008, # Uranus 1.532112481284276e-008, # Neptune 2.276247751863699e-012 ]) # Pluto