function z=gxy(x,y)
%usage:  z=gxy(x,y);

%z=sqrt(x.^2+y.^2);
z=x.^2+y.^2;