Exact L2 series solution of the Dirac-Coulomb problem for all energies

We obtain exact solution of the Dirac equation with the Coulomb potential as an infinite series of square integrable functions. This solution is for all energies, the discrete as well as the continuous. The spinor basis elements are written in terms of the confluent hypergeometric functions and chosen such that the matrix representation of the Dirac-Coulomb operator is tridiagonal. The wave equation results in a three-term recursion relation for the expansion coefficients of the wavefunction which is solved in terms of the Meixner-Pollaczek polynomials.