[en] A program is presented to calculate particle energies and wave function tables for particles with a spin of 1/2 in an arbitrary central field at given quantum numbers n, l, j (n = principal quantum number, 1 = orbital moment, j = full moment). The Dirac equation is used as a wave equation. A possibility for the variation of the radial wave function tabulation method has been incorporated. The count time shows but a slight dependence on the selection of the initial point for energy. The interrelation between program operating time and energy accuracy value is such, that with the rational wave function table form selection it is possible to achieve an accuracy of 10^-9 for a time of about 1 second. A possibility is borseen for predicting normalized wave functions at a known energy level