let $f_1(x,y)$ and $f_2(x,y)$ be two 2-dimensional Gaussian functions with means $(\mu_{1_x},\mu_{1_y}$} & $(\mu_{2_x},\mu_{2_y}$} and variances $(\sigma_{1_x},\sigma_{1_y}$} & $(\sigma_{2_x},\sigma_{2_y}$} respectively. Assuming that the Gaussian functions are not infinite , for example
$$f(x,y)=\frac{1}{2\pi\sigma^2}\int_{-a}^{a}\int_{-b}^{b}\exp\big(-({(x-\mu_x)^2+(y-\mu_y)^2})/{2\sigma^2}\big) \,\mathrm{d}x\,\mathrm{d}y$$
I am interested in finding the probability that the function $f_1$ lies within some distance $d$ of $f_2$, i.e, Pr(dist($f_1$,$f_2$)$\leq d$). Is there any way to find a Gaussian difference distribution expression for some finite Gaussian functions?
http://mathworld.wolfram.com/NormalDifferenceDistribution.html gives the Gaussian difference distribution of infinite Gaussian, but I am interested in bounded Gaussian. Please help.