I am coding a simultion of a globular cluster in python, starting with particles of the same mass and $0$ initial velocity in a random spherical distribution. I am trying to inv