From his reply tweet:

N-body simulation made with Python, parallelized with numba, and animated with matplotlib.

N=100.000. Computation time around 5h for 2.000 steps. Around 1s of compute time per step. Collisions are handled with a softening length.