Update sigmaMatrix

added function which computes sigmaMatrix for each element. At the
moment (as before they are each equal - but will allow functionality of
rigid loading conditions to be added)