
k=1; m=1;

r = [ 1 1]';

lam= r/norm(r);


kbar  = k*lam*lam';

K = [ kbar  -kbar;
     -kbar  kbar];
 
M = m*eye(4);

[V  D] = eig(K,M)