Changed everything to numerical simulation

Now, the strain at each time step and the comparison to the nonaffine
case is calculated numerically. This should fix everything!