clear
clc

k =[ 300  -200     0
  -200   500  -300
     0  -300   400]
     
     m =[  10     0     0
     0    15     0
     0     0     5]
     
     
c =[    15    -5     0
    -5    11    -6
     0    -6    10]

 
 A=[c m
     m zeros(3,3)]
 
 B=[k zeros(3,3)
     zeros(3,3) -m]
 
 [V,D]=eig(-B,A)
 
 DD=diag(D)
 
 
 V(1:3,1)
 
 V(1:3,2)
 
 pause
 V(4:6,1)*V(1,1)/V(4,1)
 
  S1=V(1:3,1)
  
  kr1=S1'*k*S1
  mr1=S1'*m*S1
  cr1=S1'*c*S1
  
  V'*A*V
  V'*B*V