% ---------------------------------------
% correlation des echantillons en x de la fonction de Henon
% a=1.4, b=0.3, x0=y0=0
% ---------------------------------------

clear;

n = input('nombre d''iterations = ');

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

x(1)=0;  y(1)=0;    % initialisation des parametres.
a=1.4; b=0.3;

% --------------------
% remplissage des vecteurs x et y
% --------------------
for i=1:n
    y(i+1)=b*x(i);
    x(i+1)=1+y(i)-a*(x(i)^2);
end

%------------------
% affichage des echantillons
%------------------
figure1=plot(x(1:length(x)-1), x(2:length(x)), 'r.');                    % ici, affiche x(i+1) en fonction de x(i).
% figure1=plot(x(1:length(x)-2), x(3:length(x)), 'r.');          -->  pour afficher  x(i+2) en fonction de x(i).

grid;
title('Correlation des echantillons x(i+1) et x(i)');
xlabel('x(i)');
ylabel('x(i+1)');
zoom;                                                                                         % permet de "zoomer et dezoomer" sur une partie du graphe avec la souris.