Assume G=1 for N-body simulation.