Lagrange stiffness flux interpolation

Added investigation of accuracy of stiffness matrix comparing a poly
interpolation of the flux to an exact solution.