%%%%%%%%%%%%%% L9_esim1 ti 21.11.06 B-sali%%%%%%%%%%%%%%%
close all
A=[.8 0;0 .64]

hold on
for i=-3:3
    for j=[-1 1]
      x0=[i; j]
      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;
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')
