def relax(self, steps=<NUM_LIT:1000>):
for step in range(steps):<EOL><INDENT>self.forces = self.force_hertzian() + self.force_damp()<EOL>self.integrate(self.forces)<EOL>self.boundary_condition()<EOL><DEDENT>
Relax the current configuration using just pair wise forces (no noise)