function k = LinearSpringElementStiffness(k) % LinearSpringElementStiffness k = [k 0 -k 0; 0 0 0 0; -k 0 k 0; 0 0 0 0]; endfunction