%%%%%%%%%%%%%% L9_esim12  ti 21.11.06 B-sali%%%%%%%%%%%%%%%
close all
A=[1.25 -0.75;-0.75 1.25];
[V,D]=eig(A)

hold on
      
      x0=[5; 5.1]
      x=x0;
      n=6;
       for k=1:n
         y=A*x(:,end);
         x=[x y];
       end;   
       plot(x(1,:),x(2,:),'o')
       plot(x(1,:),x(2,:))
       plot(x(1,1),x(2,1),'rx')
     end;
end;
for i=[-3 3]
      x0=[i; 0]
      x=x0;
      n=20;
       for k=1:n
         y=A*x(:,end);
         x=[x y];
       end;   
       plot(x(1,:),x(2,:),'o')
       plot(x(1,:),x(2,:))
       plot(x(1,1),x(2,1),'rx')
end;
title('Attraktori')
