A parallelized and efficient C++ implementation of the Barnes-Hut algorithm for simulating an N-body system.