stiffness_matrix function, basis value function