Numerical solution of the radial Dirac equation in pseudopotential construction