I once wrote a similar simulation in C: https://github.com/kgabis/gravitysim https://github.com/kgabis/gravitysim It uses Barnes-Hut algorithm, so it runs in

I once wrote a similar simulation in C: https://github.com/kgabis/gravitysim https://github.com/kgabis/gravitysim It uses Barnes-Hut algorithm, so it runs in nlogn, instead of n^2.