Elementary approximations for zeros of Bessel functions

Let j_vk, y_vk and c_vk denote the kth positive zeros of the Bessel functions J_v(x), Y_v(x) and of the general cylinder function C_v(x) = cos αJ_v(x)?sin αY_v(x), 0 ? α < π, respectively. In this paper we extend to c_vk, k = 2, 3,..., some linear inequalities presently known only for j_vk. In the case of the zeros y_vk we are able to extend these inequalities also to k = 1. Finally in the case of the first positive zero j_v1 we compare the linear enequalities given in [9] with some other known inequalities.     

Polynomial series expansions for confluent and Gaussian hypergeometric functions

Based on the Hadamard product of power series, polynomial series expansions for confluent hypergeometric functions and for Gaussian hypergeometric functions are introduced and studied. It turns out that the partial sums provide an interesting alternative for the numerical evaluation of the functions and , in particular, if the parameters are also viewed as variables.

     

The evaluation of Bessel functions via exp-arc integrals

A standard method for computing values of Bessel functions has been to use the well-known ascending series for small argument, and to use an asymptotic series for large argument; with the choice of the series changing at some appropriate argument magnitude, depending on the number of digits required. In a recent paper, D. Borwein, J. Borwein, and R. Crandall [D. Borwein, J.M. Borwein, R. Crandall, Effective Laguerre asymptotics, preprint at http://locutus.cs.dal.ca:8088/archive/00000334/] derived a series for an “exp-arc” integral which gave rise to an absolutely convergent series for the J and I Bessel functions with integral order. Such series can be rapidly evaluated via recursion and elementary operations, and provide a viable alternative to the conventional ascending-asymptotic switching. In the present work, we extend the method to deal with Bessel functions of general (non-integral) order, as well as to deal with the Y and K Bessel functions.     

Asymptotic expansions of integrals of two Bessel functions via the generalized hypergeometric and Meijer functions

Asymptotic expansions of certain finite and infinite integrals involving products of two Bessel functions of the first kind are obtained by using the generalized hypergeometric and Meijer functions. The Bessel functions involved are of arbitrary (generally different) orders, but of the same argument containing a parameter which tends to infinity. These types of integrals arise in various contexts, including wave scattering and crystallography, and are of general mathematical interest being related to the Riemann—Liouville and Hankel integrals. The results complete the asymptotic expansions derived previously by two different methods — a straightforward approach and the Mellin-transform technique. These asymptotic expansions supply practical algorithms for computing the integrals. The leading terms explicitly provide valuable analytical insight into the high-frequency behavior of the solutions to the wave-scattering problems.     

On an identity for zeros of Bessel functions

In this paper our aim is to present an elementary proof of an identity of Calogero concerning the zeros of Bessel functions of the first kind. Moreover, by using our elementary approach we present a new identity for the zeros of Bessel functions of the first kind, which in particular reduces to some other new identities. We also show that our method can be applied for the zeros of other special functions, like Struve functions of the first kind, and modified Bessel functions of the second kind.     

Numerical algorithms for uniform Airy-type asymptotic expansions

Airy-type asymptotic representations of a class of special functions are considered from a numerical point of view. It is well known that the evaluation of the coefficients of the asymptotic series near the transition point is a difficult problem. We discuss two methods for computing the asymptotic series. One method is based on expanding the coefficients of the asymptotic series in Maclaurin series. In the second method we consider auxiliary functions that can be computed more efficiently than the coefficients in the first method, and we do not need the tabulation of many coefficients. The methods are quite general, but the paper concentrates on Bessel functions, in particular on the differential equation of the Bessel functions, which has a turning point character when order and argument of the Bessel functions are equal.     

Maximal order for quadratures usingn evaluations

The problem is to find approximationsI (f; h) to the integralI(f; h)= ₀ ^h f. Such an approximation has local orderp ifI(f; h)–I (f; h)=O(h ^p ) ash0. Let(n) denote the maximal local order possible for a method usingn evaluations of a function or its derivatives. We show that (n)=2n+1 if the information used is Hermitian. This is conjectured to be true in general. The conjecture is established for all methods using three or fewer evaluations.This research was supported in part by the National Science Foundation under Grant MCS75-222-55 and the Office of Naval Research under Contract N00014-76-C-0370, NR 044-422. Numerical results reported in this paper were obtained through the computing facilities of the University of Maryland.     

Comments on the theory of Bessel functions with more than one index

We present an outline of the theory of Bessel-like functions with more than one index and one or more variables. Their link with other types of functions is discussed and their use in applications is touched on.     

A weighted uniform -estimate of Bessel functions: A note on a paper of Guo

An improved Guo's uniform estimate of Bessel functions is shown by using a uniform pointwise bound of Barceló and Córdoba.

     

Automatic,guaranteed integration of analytic functions

Davis introduced a method for estimating linear functionals of analytic functions by using Cauchy's Integral Formula. This is used to construct methods for numerical integration which give rigorous error bounds. By combining these bounds with strategies for order and subinterval adaptation, a program is developed for automatic integration of analytic functions. Interval analysis is used to validate the bounds.     

On sums of roots of infinite order entire functions

A class of infinite order entire functions is considered. Estimates for sums of the roots are derived. These estimates supplement the Hadamard theorem. Moreover, we establish a new estimate for the counting function of the roots, which in appropriate situations can be more useful than the Jensen inequality.     

Note on the Dočev-Grosswald asymptotic series for generalized Bessel polynomials

The Do?ev-Grosswald asymptotic series for the generalized Bessel polynomials y_n(z; a, b) is extended to O(1/n⁴) relative accuracy. The differential equation of the asymptotic factor, derived from the differential equation for y_n(z; a, b), is the basis of a different and easier method that employs simple recurrence relations and much less algebra for obtaining the same series. This is applied to the important special case of a = 1 to obtain the asymptotic series to O(1/n¹¹) relative accuracy.     

On a fast and accurate method for computing Fourier transforms

In this paper, a variable order method for the fast and accurate computation of the Fourier transform is presented. The increase in accuracy is achieved by applying corrections to the trapezoidal sum approximations obtained by the FFT method. It is shown that the additional computational work involved is of orderK(2m+2), wherem is a small integer andKn. Analytical expressions for the associated error is also given.     

Inequalities for formal power series and entire functions

We present several integral and exponential inequalities for formal power series and for both arbitrary entire functions of exponential type and generalized Borel transforms. They are obtained through certain limit procedures which involve the multiparameter binomial inequalities, integral inequalities for continuous functions, and weighted norm inequalities for analytic functions. Some applications to the confluent hypergeometric functions, Bessel functions, Laguerre polynomials, and trigonometric functions are discussed. Also some generalizations are given.     

Bounds on differences of adjacent zeros of Bessel functions and iterative relations between consecutive zeros

Bounds for the distance between adjacent zeros of cylinder functions are given; and are such that ; stands for the th positive zero of the cylinder (Bessel) function , , .

These bounds, together with the application of modified (global) Newton methods based on the monotonic functions and , give rise to forward ( ) and backward ( ) iterative relations between consecutive zeros of cylinder functions.

The problem of finding all the positive real zeros of Bessel functions for any real and inside an interval , 0$">, is solved in a simple way.

     

Asymptotic expansions and acceleration of convergence for higher order iteration processes

Summary We analyze the convergence behavior of sequences of real numbers {x _n }, which are defined through an iterative process of the formx _n :=T(x _n –1), whereT is a suitable real function. It will be proved that under certain mild assumptions onT, these numbersx _n possess an asymptotic (error) expansion, where the type of this expansion depends on the derivative ofT in the limit point ; this generalizes a result of G. Meinardus [6].It is well-known that the convergence of sequences, which possess an asymptotic expansion, can be accelerated significantly by application of a suitable extrapolation process. We introduce two types of such processes and study their main properties in some detail. In addition, we analyze practical aspects of the extrapolation and present the results of some numerical tests. As we shall see, even the convergence of Newton's method can be accelerated using the very simple linear extrapolation process.Dedicated to Professor Dr. Günter Meinardus on the occasion of his 65^th birthday