Reliable Computation of the Zeros of Solutions of Second Order Linear ODEs Using a Fourth Order Method

A fourth order fixed point method to compute the zeros of solutions of second order homogeneous linear ODEs is obtained from the approximate integration of the Riccati equation associated with the ODE. The method requires the evaluation of the logarithmic derivative of the function and also uses the coefficients of the ODE. An algorithm to compute with certainty all the zeros in an interval is given which provides a fast, reliable, and accurate method of computation. The method is illustrated by the computation of the zeros of Gauss hypergeometric functions (including Jacobi polynomials) and confluent hypergeometric functions (Laguerre polynomials, Hermite polynomials, and Bessel functions included) among others. The examples show that typically 4 or 5 iterations per root are enough to provide more than 100 digits of accuracy, without requiring a priori estimations of the roots.

[1]  David M. Miller,et al.  Handbook of Mathematical Functions With Formulas, Graphs and Mathematical Tables (National Bureau of Standards Applied Mathematics Series No. 55) , 1965 .

[2]  Sadayuki Murashima,et al.  Analysis of Potential Probleme and the Modified Bessel Functions of Purely Imaginary Order , 1974 .

[3]  M. Abramowitz,et al.  Handbook of Mathematical Functions With Formulas, Graphs and Mathematical Tables (National Bureau of Standards Applied Mathematics Series No. 55) , 1965 .

[4]  Ilia Krasikov.  On the zeros of polynomials and allied functions satisfying second order differential equations , 2002 .

[5]  Árpád Elbert,et al.  Some recent results on the zeros of Bessel functions and orthogonal polynomials , 2001 .

[6]  Javier Segura,et al.  Global Sturm inequalities for the real zeros of the solutions of the Gauss hypergeometric differential equation , 2007, J. Approx. Theory.

[7]  Árpád Elbert,et al.  On the Square of the Zeros of Bessel Functions , 1984 .

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

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

[10]  Xin Li,et al.  Bound on the extreme zeros of orthogonal polynomials , 1992 .

[11]  F. Olver.  Asymptotics and Special Functions , 1974 .

[12]  Knut Petras,et al.  On the computation of the Gauss-Legendre quadrature formula with a given precision , 1999 .

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

[14]  Gene H. Golub,et al.  Calculation of Gauss quadrature rules , 1967, Milestones in Matrix Computation.

[15]  T. M. Dunster.  Bessel functions of purely imaginary order, with an application to second-order linear differential equations having a large parameter , 1990 .

[16]  INEQUALITIES FOR THE SMALLEST ZEROS OF LAGUERRE POLYNOMIALS AND THEIR q-ANALOGUES , 2007 .

[17]  Nico M. Temme,et al.  Numerical methods for special functions , 2007 .

[18]  E. Yakimiw,et al.  Accurate Computation of Weights in Classical Gauss-Christoffel Quadrature Rules , 1996 .

[19]  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 .

[20]  Martin E. Muldoon,et al.  Continuous ranking of zeros of special functions , 2008 .

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

[22]  Javier Segura,et al.  The Zeros of Special Functions from a Fixed Point Method , 2002, SIAM J. Numer. Anal..

[23]  John B. Shoven,et al.  I , Edinburgh Medical and Surgical Journal.

[24]  Wolfram Koepf,et al.  Computing the Real Zeros of Hypergeometric Functions , 2004, Numerical Algorithms.

[25]  Amparo Gil,et al.  Computing the Zeros and Turning Points of Solutions of Second Order Homogeneous Linear ODEs , 2003, SIAM J. Numer. Anal..

[26]  L. Milne‐Thomson.  A Treatise on the Theory of Bessel Functions , 1945, Nature.

[27]  Amparo Gil,et al.  New inequalities from classical Sturm theorems , 2004, J. Approx. Theory.