stiffness matrix
no errors
moved all the matrix related functions out of fem.c into their own place
started a gnuplot writer