Given a bivariate normal distribution $(X_1,X_2)$ with parameters $\mu_1,\mu_2,\sigma_1,\sigma_2,\rho$ where $\rho\neq 0$.
I want to find $E(X_1X_2)$. Now $E(X_1X_2):=\int_{-\infty}^\infty \!\int_{-\infty}^\infty \! x_1x_2f(x_1,x_2) \, \mathrm{d}x_2\mathrm{d}x_1$ where $f(x_1,x_2)$ is the bivariate normal density function. Is there a short way to evaluating this integral?