I was interested in how to come up with constraints for Verlet integration. For example, one constraint to prevent collisions is to model infinitely stiff springs between particles and the penetration points of the surface. http://www.gotoandplay.it/_articles/2005/08/advCharPhysics.php