% define mass and stiffness matrices




m = 1;
k = 1e5;
M = diag([m;m;m]);
K = [2*k -k 0;-k 2*k -k;0 -k 2*k];
C = .00005*K;

%M = eye(3,3)*2;
%C = [3.8313 -1.1234 -.2440;--1.1234 3.5873 -1.3674;-.244 -1.3674 2.4639];
%K = [20000 -10000 0 ;-10000 20000 -10000;0 -10000 10000];