Add nbody simulations

How should I test this? When I run it (`python benchmark_nbody.py`), it gets stuck in an infinite loop because of the r**3 call. Otherwise, this looks like correct code.