
a=1.4;
b=-.3;
N=1000000;
x=zeros(N,1);
y=zeros(N,1);
x(1)=0;
y(1)=0;
   for n=2:N
      [x(n),y(n)]=henon(x(n-1),y(n-1),a,b);
   end % for n
       figure; grid on; hold on;
      scatter(x,y,'.');


return


