% Graficki prikaz prvih nekoliko iteracija Henonovog preslikavanja (grafik
% 1), kao i Henonov atraktor

% pocetni skup: skup tacaka (kruznica)
fi=0:0.1:2*pi;
r=0.3;
xx=r*cos(fi);
yy=r*sin(fi);

% pocetni skup: tacka
%xx=0;
%yy=0;

a=1.4;
b=0.3;

figure(1)
plot(xx,yy)
axis equal
hold on

%broj iteracija za crtanje
n=3;

xtekuce=xx;
ytekuce=yy;

% crtamo prvih n iteracija
for i=1:n
    xnovo=ytekuce+1-a*xtekuce.^2;
    ynovo=b*xtekuce;
    plot(xnovo,ynovo)
    xtekuce=xnovo;
    ytekuce=ynovo;
end

legend('(x0,y0)', '(x1,y1)', '(x2,y2)', '(x3,y3)')

% preskacemo narednih 200 iteracija
for i=1:200
    xnovo=ytekuce+1-a*xtekuce.^2;
    ynovo=b*xtekuce;
    
    xtekuce=xnovo;
    ytekuce=ynovo;
end

figure(2)
plot(xnovo,ynovo,'.r')
axis equal
hold on
N=5000; %max broj iteracija

%crtamo tacke atraktora (pocevsi od 200+n iteracije do N)
for i=1:N
    xnovo=ytekuce+1-a*xtekuce.^2;
    ynovo=b*xtekuce;
    plot(xnovo,ynovo,'.r')
    xtekuce=xnovo;
    ytekuce=ynovo;
end

