I need to show that $\int \limits_{\mathbb{R}^2}|x+ty|e^{-(x^2+y^2)/{2}}dxdy=\sqrt{8\pi (t^2+1)}$, for any given $t \in \mathbb{R}$.
I've tried the substitution $\left\{\begin{matrix} u=x+ty & \\ v= x-ty & \end{matrix}\right.$ that didn't work.
Also, polar coordinates yield: $ I=\left (\int_{0}^{\infty}r^2e^{-r^2/2}dr \right ) \left (\int_{0}^{2\pi}|\cos\theta + t\sin \theta| d\theta \right ) $ which is already a better expression but I still don't think that it can lead to an answer.
Any other ideas?
The substitution $$ u = \frac{x+t y}{\sqrt{1+t^2}},\qquad v = \frac{t x-y}{\sqrt{1+t^2}} $$ works far better, leading to
$$ \sqrt{1+t^2}\int_{\mathbb{R}^2} |u| e^{-\frac{u^2+v^2}{2}}\,du\,dv = \sqrt{8\pi(1+t^2)}\int_{0}^{+\infty}u e^{-u^2/2}\,du = \color{red}{\sqrt{8\pi (1+t^2)}} $$ as wanted, by Fubini's theorem and parity.