A slightly tweaked version of John Chambers' MERCURY N-body integrator, plus several IDL routines that are used to read and analyze the output.