Mathematica Stack Exchange is a question and answer site for users of Mathematica. It's 100% free, no registration required.

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

This question already has an answer here:

I recently made an n-body simulation and thought it would be interesting to try and model it on the interior of our solar system (Sun to Mars), but I cannot find initial conditions for such an undertaking. My biggest problem seems to be that I'm not working with a spherical coordinate system and instead have opted to using the 3D Cartesian system. Does anyone know of any data that gives initial positions and velocities of the Solar System's planets and their moons in 3D Cartesian coordinates, so I can plug this into a system of differential equations as initial conditions?

share|improve this question

marked as duplicate by Sjoerd C. de Vries, whuber, m_goldberg, Michael E2, Artes May 21 '13 at 6:17

This question has been asked before and already has an answer. If those answers do not fully address your question, please ask a new question.

FWIW, AstronomicalData[] has a nice bounty of information, including orbit paths and eccentricities. See e.g. this question. On the other hand, there's a certain appeal in using numerical integration instead of merely taking coordinates from an ephemeris... – J. M. May 20 '13 at 19:25
Just in case you are interested in accuracy, it might be worth noting that the system of ODEs described in your previous question at… is not what is usually meant by an "n-body simulation," because it does not account for interactions among the bodies: it's just a collection of independent central field solutions. It might be a fine approximation for times close to the initial time but will be observably incorrect. – whuber May 20 '13 at 19:36
@J.M.: Yeah, there's something very satisfying about using a system of ODEs to calculate each planet's path, so I'd love to have initial conditions to work with and see how the system evolves :) – user7388 May 20 '13 at 20:08
You were asking for positions and velocities. These can be found using AstronomicalData as demonstrated in that question. – Sjoerd C. de Vries May 20 '13 at 21:44
Since your question got closed... if you're interested in seeing how JPL does it, you might want to see this. – J. M. May 21 '13 at 7:44