% -----------------------
% programme d'iteration d'un carre
% par la fonction d'Henon
% -----------------------
 

clear;

input(' pour Henon : a=1.4, b=0.3, x0=y0=0 ');

% ------------------
% declaration des constantes
% ------------------
a =1.4 ;  b= 0.3;
n = input('nombre d''iterations = ');

x=zeros(1,n+1);
y=zeros(1,n+1);
 

%---------------------------
% boucle de remplissage des vecteurs x et y
% ---------------------------
for  x0= -0.28 : 0.01 : 0.28
      x(1)=x0;
   for y0= -0.28 : 0.01 : 0.28
           y(1)=y0;
                 for i=1: n
                     y(i+1)=b*x(i);
                     x(i+1)=1+y(i)-a*(x(i)^2);
                 end
           hold on;
           figure2=plot(x(n), y(n), '.');
           set(figure2,'markersize', 4);
     end
end
grid;
zoom;              % permet de "zoomer et dezoomer" sur une partie du graphe avec la souris.