REBOUND from the ground up: the outer planets for a million years

Afterwards you can build and run an N-body system in REBOUND with WHFast or IAS15, record its orbital elements, and pick the time step by the energy error.