
theta = 0:0.01:2*pi;

x = cos(theta);
y = sin(theta);


for tmp = 1:5
    
            
            xa = [];
            ya = [];        

            A = randn(2,2);

            for i = 1:length(theta)


                z = A*[x(i) ; y(i) ];
                xa = [xa ; z(1)];
                ya = [ya ; z(2)];

            end

            eig(A)
            plot(xa,ya,'.');
            grid on
            axis equal
            pause
            
end


