Dado que la función de onda para la configuración base $1s^2$ es un producto de $Ne^{-\alpha r_1}, Ne^{-\alpha r_2}$, me gustaría calcular la energía promedio $\langle e^2/4\pi\epsilon_0r_{12}\rangle$ Entonces:
$$\left\langle \frac{e^2}{4\pi\epsilon_0r_{12}} \right\rangle = \int|Ne^{-\alpha r_1}|^2\cdot|Ne^{-\alpha r_1}|^2\frac{1}{|\vec{r}_1-\vec{r}_2|}d^3r_1\,d^3r_2$$
¿Alguna buena sugerencia sobre cómo resolver esta integral?