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