 Although they differ in how they compute or apply the constraints themselves, the constraints are still modelled using Lagrange multipliers which are computed using the Gauss–Seidel method.The original SHAKE algorithm is capable of constraining both rigid and flexible molecules (eg.water, benzene and biphenyl) and introduces negligible error or energy drift into a molecular dynamics simulation.