International Journal of Mathematics and Mathematical Sciences

Volume 2012 (2012), Article ID 873078, 10 pages

doi:10.1155/2012/873078

Research Article

Approximate Closed-Form Formulas for the Zeros of the Bessel Polynomials

Rafael G. Campos and Marisol L. Calderón

Facultad de Ciencias Físico-Matemáticas, Universidad Michoacana, 58060 Morelia, MN, Mexico

Received 11 June 2012; Revised 10 September 2012; Accepted 23 September 2012

Academic Editor: Stefan Samko

Copyright © 2012 Rafael G. Campos and Marisol L. Calderón. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We find approximate expressions and for the real and imaginary parts of the th zero of the Bessel polynomial . To obtain these closed-form formulas we use the fact that the points of well-defined curves in the complex plane are limit points of the zeros of the normalized Bessel polynomials. Thus, these zeros are first computed numerically through an implementation of the electrostatic interpretation formulas and then, a fit to the real and imaginary parts as functions of , and is obtained. It is shown that the resulting complex number is -convergent to for fixed .

1. Introduction

The polynomial solutions of the differential equation

were studied systematically in [1] by the first time. They are named (generalized) Bessel polynomials and are given explicitly by

as it can be shown in [2]. Here, is the Pochhammer symbol and . Many properties as well as applications are associated to this equation; the traveling waves in the radial direction which are solutions of the wave equation in spherical coordinates can be written in terms of the polynomial solutions of (1.1). Also, this equation has application in network and filter design, isotropic turbulence fields, and more (see the monograph [2] or [3–14] and references therein for some other results). Among these, several results about the important problem concerning the location of its zeros have been obtained [8–11] and in [12], explicit expressions for sum rules and for the homogeneous product sum symmetric functions of zeros of these polynomials are given. On the other hand, the electrostatic interpretation of these zeros as the equilibrium configuration in the complex plane with a logarithmic electric potential and a dipole at the origin has been given in [13], and in [14] it is shown that this equilibrium configuration is not stable. Thus, these cases show that it is desirable to acquire new analytical knowledge about the location of the zeros of the Bessel polynomials.

In this paper we give approximate explicit formulas for both the real and imaginary parts of the th zero of and show that the approximation order of these new formulas to the exact zeros of the Bessel polynomials is for fixed .

The approach followed in this paper is simple and based on three items. The first is the electrostatic interpretation of the zeros of polynomials satisfying second-order differential equations [15–17], the second is a simple curve fitting of numerical data, and the third is the known fact that the points of well-defined curves in the complex plane are limit points of the zeros of the normalized Bessel polynomials [8–11]. The formulas yielded by the electrostatic interpretation of the zeros of Bessel polynomials are used to find them numerically as it has been done previously with these and other sets of points [7–19]. Several sets of zeros are computed in this way and the sets of real and imaginary values are fitted by polynomials depending on the index whose coefficients depend on and .

2. Asymptotic Expressions for the Zeros

Let , , be the zeros of the Bessel polynomial , ordered according to the imaginary part. Then, from (1.1) it follows that (A procedure for obtaining this kind of nonlinear equations for the zeros of a polynomial satisfying second and higher order differential equations is given in [19].)

where , that is, the real and imaginary parts of the zeros should satisfy the electrostatic equations

This set of nonlinear equations can be solved by standard methods. We have used a Newton method to solve them up to and .

Let be the th zero of the normalized Bessel polynomials , that is, and . As it is shown in Figure 1, the piecewise linear interpolation of the real and imaginary parts of can be fitted by polynomials of the second and third degree in the index .

(a)

(b)

(c)

(d)

(e)

(f)

Figure 1: Real and imaginary parts of the zeros of the normalized Bessel polynomials for . They are plotted as functions of for , in gray-level intensity, from lower to higher, according to the value of .

Thus, we propose the following expressions

for the approximate zero to fit our data. To find the dependence of the coefficients of these polynomials on and , we take into account the numerical behavior of the data at the middle and end points.

We begin by finding the coefficients of the second-order polynomial giving the real part by fitting the values of at and . In Figure 2 we show the dependence of and on for some values of .

(a)

(b)

Figure 2: Dependence of and on for . Plots are shown in gray-level intensity, from lower to higher, according to the value of .

A fit of these data to the models and yields

These conditions and the symmetry of with respect to the middle point lead to the following coefficients:

where

Now, to find the coefficients of the third-order polynomial giving the imaginary part we follow a similar procedure. The dependence of the real part of and on for some values of is shown in Figure 3.

(a)

(b)

Figure 3: Dependence of and on for . Plots are shown in gray-level intensity, from lower to higher, according to the value of .

Again, a fit of these data to the models and yields

In addition, we have that

therefore, the coefficients are given by

where

Thus, the substitution of (2.10), (2.9), (2.6), and (2.5), respectively, in (2.3) yields the approximate closed-form expressions

where , for the zeros of the unnormalized Bessel polynomial . The expressions given in (2.11) converge to the zeros of these polynomials, as we will show in the following.

3. Convergence

Following [9], we define

and denote by the curve defined by

which contains the limit points of the zeros of the normalized Bessel polynomial . Then, it has been proved in [8] that the zero of approaches to order the limit value , that is,

as .

Thus, if we show that , we will have proved that

and therefore, taking into account that , the explicit expression (2.11) approaches to order the zero of the Bessel polynomial .

To this purpose, we simply substitute the expression for given by (2.3) in (3.1) to obtain, after a lengthy calculation, that the expansion of in terms of is

where

This implies that

for fixed and . Thus, approaches to order the curve and (3.4) follows. From here we have that

as . Numerical calculations confirm and extend this result. Figure 4 shows the behavior of the maxima of over as they depend on for the particular case of . The numbers computed by (2.2) are taken as the exact zeros . A fit of these data gives with .

Figure 4: Plot of the values of against for the case of .

4. Some Few Tests

Just to give examples of the application of the approximate expression (2.11), we consider the following cases.

4.1. The Real Zero

A closed-form formula for the unique real zero of the Bessel polynomial can be obtained by the substitution of in the real part of (2.11), . This gives

as our new result. In [2, 11] very accurate expressions for are given. Particularly, the following formula:

is given in [11]. Expanding in powers of we find that

indicating good relative agreement between the two results.

4.2. Power Sums

Here we carry out the corresponding multiplications and use some cases of Faulhaber’s formula. Then we compare our results with the exact ones.

(1)  Sum of the Zeros. The simple sum of (cf. (2.11)) gives a complicated expression for the real part. However, expanding both the real and imaginary parts of this sum gives

The exact result is , as can be seen from (1.2).

(2)  Sum of the Squares of the Zeros. In this case we take the particular case of . For this value we obtain

The exact sum

can be found elsewhere [12].

The use of the approximate formula (2.11) for obtaining sums of higher powers of these zeros is not expected to give satisfactory results, since the powers and the sum magnify the total error.

5. Final Comment

The approximate formula for given above is far from being unique. There exist many other functions to fit the zeros obtained through the electrostatic equations (2.2), and there are other conditions to impose at the extreme and middle points of the fitting interval. For instance, the imaginary part can be fitted by a polynomial of degree 5, but this does not improve the rate of convergence and, on the other hand, the calculations become more complicated.

Acknowledgment

The authors thank Consejo Nacional de Ciencia y Tecnología for the financial support given to this project.

References

H. L. Krall and O. Frink, “A new class of orthogonal polynomials: the Bessel polynomials,” Transactions of the American Mathematical Society, vol. 65, pp. 100–115, 1949.

E. Grosswald, Bessel Polynomials, vol. 698 of Lecture Notes in Mathematics, Springer, Berlin, Germany, 1978.

H. M. Srivastava, “Some orthogonal polynomials representing the energy spectral functions for a family of isotropic turbulence fields,” Zeitschrift für Angewandte Mathematik und Mechanik, vol. 64, no. 6, pp. 255–257, 1984.

J. L. López and N. M. Temme, “Large degree asymptotics of generalized Bessel polynomials,” Journal of Mathematical Analysis and Applications, vol. 377, no. 1, pp. 30–42, 2011.

Ö. Eğecioğlu, “Bessel polynomials and the partial sums of the exponential series,” SIAM Journal on Discrete Mathematics, vol. 24, no. 4, pp. 1753–1762, 2010.

C. Berg and C. Vignat, “Linearization coefficients of Bessel polynomials and properties of Student -distributions,” Constructive Approximation, vol. 27, no. 1, pp. 15–32, 2008.

L. Pasquini, “Accurate computation of the zeros of the generalized Bessel polynomials,” Numerische Mathematik, vol. 86, no. 3, pp. 507–538, 2000.

A. J. Carpenter, “Asymptotics for the zeros of the generalized Bessel polynomials,” Numerische Mathematik, vol. 62, no. 4, pp. 465–482, 1992.

F. W. J. Olver, “The asymptotic expansion of Bessel functions of large order,” Philosophical Transactions of the Royal Society of London Series A, vol. 247, pp. 328–368, 1954.

H.-J. Runckel, “Zero-free parabolic regions for polynomials with complex coefficients,” Proceedings of the American Mathematical Society, vol. 88, no. 2, pp. 299–304, 1983.

M. G. de Bruin, E. B. Saff, and R. S. Varga, “On the zeros of generalized Besselpolynomials I, II,” Indagationes Mathematicae, vol. 84, pp. 1–25, 1981.

F. Gálvez and J. S. Dehesa, “Some open problems of generalised Bessel polynomials,” Journal of Physics A, vol. 17, no. 14, pp. 2759–2766, 1984.

E. Hendriksen and H. van Rossum, “Electrostatic interpretation of zeros,” in Orthogonal Polynomials and Their Applications, M. Alfaro, J. S. Dehesa, F. J. Marcellan, J. L. Rubio de Francia, and J. Vinuesa, Eds., vol. 1329 of Lecture Notes in Mathematics, pp. 241–250, Springer, Berlin, Germany, 1988.

G. Valent and W. Van Assche, “The impact of Stieltjes' work on continued fractions and orthogonal polynomials: additional material,” Journal of Computational and Applied Mathematics, vol. 65, no. 1–3, pp. 419–447, 1995.

G. Szegő, Orthogonal Polynomials, Colloquium Publications, American Mathematical Society, Providence, RI, USA, 1975.

F. Marcellán, A. Martínez-Finkelshtein, and P. Martinez-Gonzalez, “Electrostatic models for zeros of polynomials: old, new, and some open problems,” Journal of Computational and Applied Mathematics, vol. 207, pp. 258–272, 2007.

R. G. Campos, “Perturbed zeros of classical orthogonal polynomials,” Boletin de la Sociedad Matematica Mexicana, vol. 5, no. 1, pp. 143–153, 1999.

R. G. Campos, “Solving singular nonlinear two-point boundary value problems,” Boletin de la Sociedad Matematica Mexicana, vol. 3, no. 2, pp. 279–297, 1997.

R. G. Campos and L. A. Avila, “Some properties of orthogonal polynomials satisfying fourth order differential equations,” Glasgow Mathematical Journal, vol. 37, no. 1, pp. 105–113, 1995.