The Zeros of Special Functions from a Fixed Point Method

A scheme for the computation of the zeros of special functions and orthogonal polynomials is developed. We study the structure of the first order difference-differential equations (DDEs) satisfied by two fundamental sets of solutions of second order ODEs $y_n^{\prime\prime}(x)+A_{n}(x)y_n (x)=0$, n being the order of the solutions and An (x) a family of continuous functions. It is proved that, with a convenient normalization of the solutions, $T_{\pm1}(z)=z\pm \mbox{\rm sign}(d)\arctan(y_n (x(z))/y_{n\pm 1}(x(z)))$ are globally convergent iterations with fixed points $z(x_n^{(i)})$, $x_n^{(i)}$ being the zeros of yn(x); d is one of the coefficients in the DDEs and z(x) is a primitive of d. The structure of the DDEs is also used to set global bounds on the differences between adjacent zeros of functions of consecutive orders and to find iteration steps which guarantee that all the zeros inside a given interval can be found with certainty. As an illustration, we describe how to implement this scheme for the calculation of the zeros of arbitrary solutions of the Bessel, Coulomb, Legendre, Hermite, and Laguerre equations.

[1]  Issei Fujishiro,et al.  Computing zeros and orders of Bessel functions , 1991 .

[2]  Michael N. Vrahatis,et al.  RFSFNS: A portable package for the numerical determination of the number and the calculation of roots of Bessel functions , 1995 .

[3]  Javier Segura,et al.  A global Newton method for the zeros of cylinder functions , 2004, Numerical Algorithms.

[4]  Yasuhiko Ikebe,et al.  The zeros of regular Coulomb wave functions and of their derivatives , 1975 .

[5]  R. Piessens,et al.  ON THE COMPUTATION OF ZEROS AND TURNING POINTS OF BESSEL FUNCTIONS , 1993 .

[6]  Michael N. Vrahatis,et al.  On the Localization and Computation of Zeros of Bessel Functions , 1997 .

[7]  DongSheng Cai,et al.  Error analysis for the computation of zeros of regular Coulomb wave function and its first derivative , 2001, Math. Comput..

[8]  Paolo Ricci,et al.  On asymptotic formulas for the first zero of the Bessel function J v , 1996 .

[9]  Michael N. Vrahatis,et al.  Locating and Computing Zeros of Airy Functions , 1996 .

[10]  Árpád Elbert,et al.  An upper bound for the zeros of the cylinder function C_v(x) , 1998 .

[11]  Milton Abramowitz,et al.  Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables , 1964 .

[12]  N. Temme.  Special Functions: An Introduction to the Classical Functions of Mathematical Physics , 1996 .

[13]  Martin E. Muldoon,et al.  Convexity Properties of Special Functions and Their Zeros , 1998 .

[14]  J. Gard,et al.  Method for Evaluation of Zeros of Bessel Functions , 1973 .

[15]  Amparo Gil,et al.  ELF and GNOME: Two tiny codes to evaluate the real zeros of the Bessel functions of the first kind for real orders , 1999 .

[16]  F. W. J. Olver,et al.  On the Reversion of an Asymptotic Expansion and the Zeros of the Airy Functions , 1999, SIAM Rev..

[17]  Nico M. Temme,et al.  An algorithm with ALGOL 60 program for the computation of the zeros of ordinary bessel functions and those of their derivatives , 1979 .

[18]  Michael N. Vrahatis,et al.  Locating and Computing All the Simple Roots and Extrema of a Function , 1996, SIAM J. Sci. Comput..

[19]  D. Owen.  Handbook of Mathematical Functions with Formulas , 1965 .

[20]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[21]  James S. Ball.  Automatic Computation of Zeros of Bessel Functions and Other Special Functions , 1999, SIAM J. Sci. Comput..

[22]  Javier Segura,et al.  Bounds on differences of adjacent zeros of Bessel functions and iterative relations between consecutive zeros , 2001, Math. Comput..

[23]  P. Siafarikas,et al.  A differential inequality for the positive zeros of Bessel functions , 1992 .

[24]  Imanuel Marx.  On the structure of recurrence relations. , 1953 .

[25]  Robert Piessens,et al.  Chebyshev series approximations for the zeros of the Bessel functions , 1984 .