A=2*rand(2,2)-1;
%A = A*A';
%A=[1 -5;5, 1];
x = [1,1]';
z=[]
for i = 1:100
    x = A*x;
    z = [z,x/norm(x)];
end
plot(z(1,:),z(2,:),'r*')
eig(A)