function [s,t]=henon(x,y,a,b)
   s=a-b*y-x^2;
   t=x;
return

