The paper results are derived before this commit. For making element_stiff_matrix.py a little faster, in physical_crd_xi function a line is added to avoid
calling function several times. In element_stiff_matrix function the location of lag_xi2 and der_lag_dxi2 changed and moved to outer loop to mae it faster.