We can now shear the material using Lees-Edwards boundaries. The program now shears gradually over an angle phi+dphi and outputs the stretching energy as a function of each step. It also outputs the new coordinates.