% Section 1.7  Brin & Stuck
   close all;
   A=[2,1;1,1];
   N=1000;
   x=.1*rand(N,1)+.5;
   y=.1*rand(N,1)+.5;
   figure; hold on;
   plot([0;1;1;0;0],[0;0;1;1;0]);
   scatter(x,y,'.');
   for iter=1:10
      for k=1:N
         v=A*[x(k);y(k)];
	 x(k)=v(1)-floor(v(1));
	 y(k)=v(2)-floor(v(2));
      end % for k
      figure; hold on;
      plot([0;1;1;0;0],[0;0;1;1;0]);
      scatter(x,y,'.');
   end % for iter


