So I've modeled the interactions between the sun and all the planets (and the interactions between the planets) using Verlet integration.
I've used data from Wikipedia for masses, distance from the sun etc.
I initialized the initial velocities of the planets via the critical velocity equation.
This produces nice stable velocities.
I'm unsure of how to calculate the initial velocity of the moon so that it stays in orbit around the earth.