I am faced with the following problem:
$$Q(t) = \int_0^t \int_0^1 \int_0^1 \exp\left(-\frac{(x-x_s)^2+(y-y_s)^2}{\omega^2}\right) \,dx \, dy \, d\tau.$$
I know the first trivial step, since the integrand is not dependent of $\tau$, then we get
$$Q(t) = t \int_0^1 \int_0^1 \exp\left(-\frac{(x-x_s)^2+(y-y_s)^2}{\omega^2}\right) \,dx\,dy,$$
but that just kicks the can a little bit further down the road. It reminds of a Gaussian integral, but this contains two variables $x$ and $y$, not just one. I have been trying with a variable transformation by changing it to polar coordinates. Then I end up with some unknown boundary. What is the trick to this? Could anybody give me hint/solution?
Best regards//