clc
clear all

m1 = 0.4;
m2 = 0.01;
m3 = 0.078;

c1 = 1.00;
c2 = 36.00;
c3 = 1.53;

k1 = 0.01;
k2 = 3667;
k3 = 1152;

M = diag([m1 m2 m3]);

C = [-(c1+c2)    c2         0  ;
      c2      -(c2+c3)     c3  ;
      0          c3       -c3  ];
  
K =[-(k1+k2)     k2         0 ;
    k2       -(k2+k3)      k3  ;
      0          k3       -k3 ];