Changed read_nbody to use astropy.io.fits. Minor docs changes.