  A relativistic extension of our pseudo-shifted $\ell$-expansion technique is
presented to solve for the eigenvalues of Dirac and Klein-Gordon equations.
Once more we show the numerical usefulness of its results via comparison with
available numerical integration data.
