The radial part of the Schrödinger equation for a system of two particles, of
charges +e and -e, interacting through a Coulombic potential is:
$$-\frac{\hbar^2}{2\mu r^2} \frac{d}{dr} \left( r^2 \frac{dR(r)}{dr} \right) - \frac{e^2 R(r)}{4\pi \epsilon_0 r} + \frac{\hbar^2}{2\mu r^2} l(l+1)R(r) = ER(r).$$
(a) Show that the function:
$$R(r) = \frac{1}{4\sqrt{2\pi}} \left( \frac{1}{a_0} \right)^{3/2} \left( 2 - \frac{r}{a_0} \right) exp \left( -\frac{r}{2a_0} \right),$$
with l = 0, is a solution of the radial Schrödinger equation with the energy
eigenvalue $$E = -\frac{\mu e^4}{8h^2 (4\pi \epsilon_0)^2}.$$ Take $$a_0 = \frac{\hbar^2 (4\pi \epsilon_0)}{\mu e^2}.$$