Given $x \sim \mathcal{N}(x|0;1)$ and $y \sim \mathcal{N}(y|1;1)$,$x,y$ are two independent variables. How to find the expectation $\int \int (1+|x-y| )e^{-|x-y|}\mathcal{N}(x|0;1)\mathcal{N}(y|1;1) dx dy$. In other words, how to find the expecation of $f=(1+|x-y| )e^{-|x-y|}$ when $x,y$ are Gaussian random variables.
The absolute value makes the problem somewhat difficult.
Thank you so much