This is an old revision of the document!
* Check numerically that the constructed HO wf's are orthonormal. I.e., evaluate $\int r^2dr R_{nl}(r)R_{n'l'}(r)$. You will want to use Gaussian quadrature to discretize the integrals. Gaussian quadrature is discussed some in . While Gauss-Laguerre quadrature is ideal for this problem, plain Gauss-Legendre quadrature, which is widely available in canned routines you can easily find via Google, is sufficient.
* Next, write a function/subroutine that calculates the matrix elements of the Coulomb potential, $\langle nl|V|n'l\rangle$. Use atomic units ($e=m_e=1$) where $V(r)=1/r$.
\begin{equation} \langle nl|V|n'l\rangle = \int_0^{\infty}r^2dr R_{nl}(r}\frac{1}{r}R_{n'l}(r) \end{equation}
Here are