How can ı find mode shape from natural frequency

My codes are in below. But the results are not that ı want. What can I do?

A = [ 1 2 3 4 ;
4 5 6 7 ;
7 8 9 10 ;
10 11 12 13 ];

B =  [ 5 7 8 4 ;
4 5 6 7 ;
7 8 9 10 ;
10 11 12 13 ];

[v,d]=eig(B,A);
w=sqrt(d)
w1=w(1,1)
w2=w(2,2)
w3=w(3,3)
w4=w(4,4)
disp([w1;w2;w3;w4])

for i=1:4
    v(:,i)=v(:,i)/v(2,i)
end
disp(v)

Answer 1

[KSSV]

[v,d]=eig(B,A);
plot(v)