function Kloc=get_local_stiff(ien,xyz,etype)
% Kloc=get_local_stiff(id,ien,xyz)
% compute local stiffness matrices
[mnod,nel]=size(ien);
for i=1:nel