c1 = 5; c2 = 25e3; m1 = [1 0 0 ; 0 2 0 ; 0 0 1]; M = c1*m1; k1 = [4 -3 0 ; -3 4 -1 ; 0 -1 3]; K = c2*k1; [v,d] = eig(K,M); w1 = sqrt(d(1,1​))*sqrt(c2​/c1); w2 = sqrt(d(2,2))*sqr

c1 = 5; 
c2 = 25e3; 
m1 = [1 0 0 ; 0 2 0 ; 0 0 1]; 
M = c1*m1; 
k1 = [4 -3 0 ; -3 4 -1 ; 0 -1 3]; 
K = c2*k1; 
[v,d] = eig(K,M); 
w1 = sqrt(d(1,1))*sqrt(c2/c1); 
w2 = sqrt(d(2,2))*sqrt(c2/c1); 
w3 = sqrt(d(3,3))*sqrt(c2/c1); %to plot modes of vib.&frequencies: 
hold on
x=[0 1 2 3 4]; 
y1=[0 v(1,1) v(2,1) v(3,1) 0]; %1st mode of vib. 
plot(x,y1) 
y2=[0 v(1,2) v(2,2) v(3,2) 0]; %2nd mode of vib. 
plot(x,y2) 
y3=[0 v(1,3) v(2,3) v(3,3) 0]; %3rd mode of vib. 
plot(x,y3)
hold off

[James Tursa]

Please edit your post to make it readable, and ask a question. Also fix the title to something meaningful.