Program1_1
% -----------------------
% 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.