The Yukawa or Debye-Hueckel potential, by a suitable choice of scale factor (d) and the use of the corresponding energy unit (h/sup 2//2md/sup 2/), may be written in the form V(r:Z) = -(2Z/r)e/sup -r/. For Z much less than 1 or l much greater than 1 or E = k/sup 2/ sufficiently large the partial wave phase shifts are given by the Born approximation which mathematically vanishes at E = 0. On the other hand, Levinson's theorem indicates that a phase shift, in units of ..pi.., is equal to n(Z,l), the integral number of bound states sustainable by the potential. In this work an optical model code written by Relyea for molecular potentials is adapted to be suitable for singular attractive potentials such as the Yukawa potential. This modified code is used to obtain the numerical Schroedinger phase shifts for a broad range of Z and l values. By the addition of integers to assure continuity of the phase shifts and the satisfaction of Levinson's theorem the authors obtain an array of continuous functions which, in effect, interpolate between the Levinson limit and the Born limit. Next, they find an accurate analytic characterization of the number of bound states n/submore » a/ (Z,l) suggested by the Sommerfeld model. Finally, they develop a simple two-term analytic formula which interpolates between the Levinson limit and the Born limit and fits the rectified set of Schroedinger phase shifts quite well. In effect, they arrive at a four-dimensional function delta(E,l,Z,d) that should be useful as a convenient approximation in atomic, nuclear, particle, and plasma physics. 12 references, 9 figures, 2 tables.« less