function henonmap(iterations)

X=zeros(1,iterations);
Y=X; 
X(1)=0; 
Y(1)=0;
 
for k=2:iterations    
        X(k)= Y(k-1)+1-1.4*(X(k-1))^2;
        Y(k) = 0.3*X(k-1);
        
        if iterations>100
            plot(X,Y,'b.','MarkerSize',5);
            xlabel('X')
            ylabel('Y')
            title('Henon Map')
            % grid;
            end               
end

% "Complex and Chaotic Nonlinear Dynamics. 
% Advances in Economics and Finance, 
% Mathematics and Statistics" 
% T.Vialar, Springer 2009
% Copyright(c).