function Z = fxyexp (X, Y);

  Z = X .* Y .* exp (-X.^2-Y.^2);

