Program1




% --------------------------------------
% programme de resolution de la fonction d'Henon
% possibilite d'ajuster d'autres valeurs pour a, b, n , x0 et y0
% --------------------------------------
 

clear;

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

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

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

%---------------------------
% boucle de 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

figure1=plot(x, y, '.');

title('x0=0 , y0=0 , a=1.4, b=0.3, n=100000');     % ajuster a, b, n,  x(0) et y(0)  selon les valeurs choisies
xlabel('x');
ylabel('y');
grid;
zoom;              % permet de "zoomer et dezoomer" sur une partie du graphe avec la souris.


back to current page