  A new representation for a regular solution of the radial Dirac system of a
special form is obtained. The solution is represented as a Neumann series of
Bessel functions uniformly convergent with respect to the spectral parameter.
For the coefficients of the series convenient for numerical computation
recurrent integration formulas are given. Numerical examples are presented.
