function [x,y]=henon(x0,y0)
%Henon map
%a=1.4;
%b=0.3;
x=1.4-x0*x0+0.3*y0;
y=x0;