function y = computehenon(N, a, b, y0)

% changing the values of a and b chnages the plot
% this shows that henon is sensitive to initial
% conditions
%

y = y0;

for i=3:N
    y=[y,a-y(i-1)^2+b*y(i-2)];
end
    