// matlabfunction [Ke, Me, TLG] = MITC4(Young, Poisson, thick, density, node, Vn) nNPE = 4; % # of nodes per element nDPN = 5; % # of DOFs per node nDIM = 3; % Problem dimension nGP = 2; % # of Gauss points % Vn = Make_Vn(node); [V1,V2] = Director(Vn(1,:)); % norn vec 구하기 TLG = [V1' V2' Vn(1,:)']; Ke = Shell_Stiffness(Young, Poisson, thick, node, Vn); ..