function z=gxy(x,y)
%usage:  z=gxy(x,y);
z=ones(size(x));
%z=1+x*0+y*0;