Trying to solve $\frac{1}{\sigma\sqrt{2\pi}}\int_r^\infty xe^{-\frac{1}{2}(\frac{x-\mu}{\sigma})^2}dx$ I'm looking to simplify/solve $$\frac{1}{\sigma\sqrt{2\pi}}\int_r^\infty xe^{-\frac{1}{2}(\frac{x-\mu}{\sigma})^2}dx$$
, where $\frac{1}{\sigma\sqrt{2\pi}}e^{-\frac{1}{2}(\frac{x-\mu}{\sigma})^2}$ is the pdf of the normal distribution.
I tried to look at solving it with substitution, but I struggle at a certain point.
So if I define $u=-\frac{1}{2}(\frac{x-\mu}{\sigma})^2$, then $du = -\frac{x-\mu}{\sigma^2}dx$
Now, I want to come to:
$$\frac{1}{\sigma\sqrt{2\pi}}\int_r^\infty e^udu$$ but is the following correct?
$$\mu(1-F(r))- \frac{1}{\sigma^3 \sqrt{2\pi}}\int_r^\infty e^udu$$, where $F(.)$ is the CDF of the normal distribution.
I believe I'm making a mistake here, but just don't see what. Any help or hints are appreciated.
EDIT
I’m not looking for a closed form solution but a solution in terms of the pdf and cdf of the normal distribution.

Let $\phi(x)=\frac{1}{\sqrt{2\pi}}e^{-x^2/2}$ and $\Phi(x)=\int_{-\infty}^x\phi(t)\,\mathrm dt$ denote the standard normal density and distribution, respectively. Now consider the slightly more general problem
$$
I_b=\int_r^bx\frac{1}{\sigma}\phi\left(\frac{x-\mu}{\sigma}\right)\,\mathrm dx.
$$
If we define $\rho=(r-\mu)/\sigma$, $\beta=(b-\mu)/\sigma$, and $Z=\Phi(\beta)-\Phi(\rho)$ then the problem can be written as
$$
I_b=Z\int_r^bx\frac{\frac{1}{\sigma}\phi\left(\frac{x-\mu}{\sigma}\right)}{Z}\,\mathrm dx=Z\mathsf E(X|r\leq X\leq b),
$$
which is just the first moment of the truncated gaussian distribution scaled by $Z$.  Using the first moment of the trucnated gaussian we obtain
$$
I_b=\mu(\Phi(\beta)-\Phi(\rho))+\sigma(\phi(\rho)-\phi(\beta)).
$$
For the specific problem at hand, $b=\infty$ and upon taking taking the limit we obtain
$$
I_\infty=\mu(1-\Phi(\rho))+\sigma\phi(\rho).
$$
