The task of obtaining the wavefunctions is solvable in almost any Quantum Chemistry (e.g., HyperChem, Gaussian, etc.) or DFT (e.g., VASP, Abinit, etc.) software package (open-source alternatives should be able to handle this task just as well). In Quantum Chemistry packages it should be as straightforward as placing an atom, selecting the desired method, basis, and clicking 'compute'. As for DFT packages, it can get trickier, for example you may have to simulate the single atom computation by considering a crystal with a very large unit cell.
As for obtaining the radial part of the wavefunction - you may need to export the wavefunction and integrate it elsewhere.