Finished analytic computation of stress matrix.

Inner product $(\varphi^i, \varphi^j)$ was computed numerically
but the integrals were innacurate and higher precision gave me
floating point rounding errors. Switched to analytic solution.