function iterate(x0,y0,N)
x(1)=x0;
y(1)=y0;
	for t=1:N,
		[x(t+1),y(t+1)]=henon(x(t),y(t));
	end
	plot(x,y,'.');