slender_column.m

%matrix of stiffness values
k= [10,20,20,10]; 
%matrix of mass values 
m= [1;2;4;];

%vector representation of force functions 
F1=[-1*(k(1)+k(2)),k(2),0];
F2=[k(2),-1*(k(2)+k(3)),k(3)];
F3=[0,k(3),-1*(k(3)+k(4))];

K= [F1;F2;F3];
K= -1*K;

M= diag(m);

Ktilde= (M^(-1/2))*(K)*(M^(-1/2)); %function 

e=eig(Ktilde) %solving for eigenvalues (lambda)

W_n= sqrt(e) %solving for omega