N-body solver, written for HPCSE I (2013)