function g = g2x(x,y,s2)

sdnorm = 1 ./ (2*pi*s2.^4);
g = sdnorm .* (((x/s2).^2) - 1) .* exp(-(x.^2 + y.^2)./(2*s2.^2)); 
