Accurate numerical solution of the Schrödinger and Dirac wave equations for central fields. Nature of problem: This subroutine package provides numerical solutions of the Schrodinger and Dirac equations for central fields such that V(r)->0 when r-> infinity and rV(r) is finite for r=0. Radial wave functions, eigenvalues for bound states and phase shifts for free states are evaluated with a prescribed accuracy. Results are delivered to the main program after calling a single subroutine.