int(int((1/(2*pi*sigma^2)*exp(-(i^2 + j^2)/(2*sigma^2)))^2, i = -inf..inf), j=-inf..inf)