clc clear m1=5; m2=5; m3=5; % kg k1=100; k2=200; k3=300; % N/m M=[m1 0 0; 0 m2 0; 0 0 m3]; K=[k1+k2 -k2 0; -k2 k2+k3 -k3; 0 -k3 k3]; A=inv(M)*K; [eigVector, eigValue]=eig(A); Omega=sqrt(eigValue) ModeShapes=eigVector