I'm working on a visual simulation of the solar system for presentation to a homeschool group. I'm using GlowScript (very cool - a lot like VPython). The theme of the presentation is "Why don't things line up?" and I want to highlight how axial rotations and orbital planes don't really line up in our solar system. I can see from nasa fact pages the angle of axial tilt for each planet, which I assume is the angle between the rotation axis and the orbital plane of the planet (as opposed to the ecliptic). (ex: http://nssdc.gsfc.nasa.gov/planetary/factsheet/venusfact.html) Here's my question: what is the vector that represents these axes in 3D? Right now all my axes point in the same (roughtly celestial north) direction, which is less interesting than real-life, I assume. I see a section titled "North Pole of Rotation" which gives an RA and Dec for each planet. Does that mean if I draw a straight line through the axis of each planet it points to that part of the sky? (I don't really believe that because things move, so it seems those coordinates would be specific to this moment in time.) I'm looking for suggestions if anyone knows a good way to define the axes of our solar system in 3D. (I'm using an XYZ coordinate system centered on the center of the sun.) One other question, for those who may be inclined - I'm trying to build my model so that the motions are represented. I'm not building the model to scale so its easier to see, which means I'm not using gravity to model the motions. Do you know a good way to model the motions of the planets (and moons in particular) using only velocity vectors (no acceleration vectors)? There is probably a better way to ask the question I'm asking. I just want to model things moving in a realistic way without having to introduce gravity. Thanks! Dion