simple N-body solver integrated into Python