Just though to point this out since I myself was wondering why my simulation crashes with many balls at high framerates. It is crucial for this verlet approach NOT to skip any frames since each will be needed for the calculation (especially with the naming convetion of this video).
