Lot of modifications to improve stability.

Notable changes are:
  - friction follows Coulomb model
  - friction and restitution coefficients are combined by a multiplication (instead of a mean)
  - restitution is enabled
  - there are now two integration passes: one before constraints resolution and one after