Advanced structural systems - mechanical engineering

clc clear m1=1; m2=1; k1=100; k3=100; k2=400; M=[m1 0; 0 m2]; K=[ k1+k2 -k2 ; -k2 k2+k3 ]; A=inv(M)*K; [eigVector, eigValue]=eig(A); Omega=sqrt(eigValue) ModeShapes=eigVector