Implemented analytic Jacobian evaluation for Guccione active contraction

[MarioHabenbacher]

Updated analytic evaluation of the Jacobian for the Guccione active contraction to incorporate length dependence:

  lambda(1) = 1/(SQRT(2.0_DP*E(1)+1)) !deriv of lambda wrt E
  lambda(2) = 1/(SQRT(2.0_DP*E(2)+1)) !deriv of lambda wrt E
  lambda(3) = 1/(SQRT(2.0_DP*E(3)+1)) !deriv of lambda wrt E

...

      STRESS_TENSOR(component_idx)=STRESS_TENSOR(component_idx)+VALUE*(1.0_DP+1.45_DP*(lambda(component_idx)-1.0_DP))