numits=80000;
a(1)=0.0;b(1)=-.08;c(1)=.62;d(1)=0.0;e(1)=-7.2;f(1)=5.0;
a(2)=-.07;b(2)=.05;c(2)=.39;d(2)=.005;e(2)=-7.3;f(2)=5.0;
a(3)=.07;b(3)=-.05;c(3)=.39;d(3)=.005;e(3)=-5.3;f(3)=5.0;
a(4)=0.0;b(4)=-.08;c(4)=.62;d(4)=0.0;e(4)=-4.4;f(4)=5.0;
a(5)=.10;b(5)=-.08;c(5)=.59;d(5)=.01;e(5)=-2.1;f(5)=4.7;
a(6)=.2;b(6)=0.0;c(6)=0.0;d(6)=.08;e(6)=-1.5;f(6)=3.9;
a(7)=-.10;b(7)=.08;c(7)=.59;d(7)=.01;e(7)=-1.1;f(7)=4.7;
a(8)=0.0;b(8)=-.08;c(8)=.59;d(8)=0.0;e(8)=2.4;f(8)=4.5;
a(9)=.215;b(9)=0.0;c(9)=0.0;d(9)=.08 ;e(9)=2.0;f(9)=9.1;
a(10)=0.0;b(10)=-.08;c(10)=.62;d(10)=0.0;e(10)=5.3;f(10)=5.0;
a(11)=.215;b(11)=0.0;c(11)=0.0;d(11)=.08;e(11)=6.4;f(11)=4.5;
a(12)=0.0;b(12)=-.08;c(12)=.625;d(12)=0.0;e(12)=8.0;f(12)=5.0;

y=zeros(2,numits);
x=rand(2,1);
z=zeros(2,1);

%We'll be throwing away the first few points...
for i=1:100
   k=ceil(12*rand);
      z(1)=a(k)*x(1)+b(k)*x(2)+e(k);
      z(2)=c(k)*x(1)+d(k)*x(2)+f(k);
      x=z;   
end

%Keep these for plotting
for i=1:numits
   k=ceil(12*rand);
   z(1)=a(k)*x(1)+b(k)*x(2)+e(k);
   z(2)=c(k)*x(1)+d(k)*x(2)+f(k);
   y(:,i)=z;
   x=z;   
end