Another two functions

1. Assemble the stiffness matrix.
2. Implement the Gram_Schmidt orthonormalization for the use of
computing the warping parameter Re.