gravitational collapse of spongebob

N-body simulation made with Python, parallelized with numba, and animated with matplotlib. N=. Computation time around 5h for steps. Around 1s of compute time per step. Collisions are handled with a softening length. Code will be on my github, eventually