tested the non-optimized n-body solver. It is working correctly. See movie