Arxiv:1708.09715V1 [Math.CA] 31 Aug 2017
Total Page:16
File Type:pdf, Size:1020Kb
INTEGRAL REPRESENTATION FOR BESSEL'S FUNCTIONS OF THE FIRST KIND AND NEUMANN SERIES ENRICO DE MICHELI Abstract. A Fourier-type integral representation for Bessel's function of the first kind and complex order is obtained by using the Gegenbuaer extension of Poisson's integral representation for the Bessel function along with a trigono- metric integral representation of Gegenbauer's polynomials. This representa- tion lets us express various functions related to the incomplete gamma function in series of Bessel's functions. Neumann series of Bessel functions are also con- sidered and a new closed-form integral representation for this class of series is given. The density function of this representation is simply the analytic function on the unit circle associated with the sequence of coefficients of the Neumann series. Examples of new closed-form integral representations of spe- cial functions are also presented. 1. Introduction The Bessel function of the first kind and order ν is defined by the series [27, Eq. 8, p. 40] 1 : z ν X ( 1)k (z=2)2k (1.1) Jν (z) = − ; (z C; ν C; ν = 1; 2;:::); 2 k! Γ(ν + k + 1) 2 2 6 − − k=0 which is convergent absolutely and uniformly in any closed domain of z and in any bounded domain of ν. It is the solution of Bessel's equation d2y dy (1.2) z2 + z + (z2 ν2) = 0; dz2 dz − which is nonsingular at z = 0. The function Jν (z) is therefore an analytic function of z for any z, except for the branch point z = 0 if ν is not an integer. In this paper, a new integral representation for the Bessel functions of the first 1 kind and order ν C, Re ν > 2 , is obtained. This representation, which is given in Section 3, generalizes2 to complex− values of the order ν the classical Bessel's arXiv:1708.09715v1 [math.CA] 31 Aug 2017 integral n Z π ( i) iz cos θ (1.3) Jn(z) = − e cos nθ dθ (n Z); π 0 2 which holds only for integral order values (see, e.g., formula (3.21)). The Bessel functions of the first kind are expressed in terms of an integral of Fourier-type, which involves the regularized incomplete gamma function P ( ν ; z), with ν be- ing the (complex) fractional part of ν. This result is achievedf byg using the Poissonf g 1 integral representation of Jν+n(z) (Re ν > , ` = 0; 1; 2;:::) given in terms of − 2 2010 Mathematics Subject Classification. 33C10,40C10,33B20. Key words and phrases. Bessel function of the first kind, Neumann series of Bessel functions, Integral representations. 1 2 E. DE MICHELI Gegenbauer polynomials and then by exploiting the integral representation of the latter polynomials that we prove in Section 2. The Fourier form of the integral representation of Jν+n(z) motivates us to con- sider the Fourier inversion formula in Section 4. This analysis allows us to obtain, for either integer and half-integer orders ν, a trigonometric expansion of the lower incomplete gamma function, whose coefficients are related to (modified and unmod- ified) Bessel functions of the first kind. Thus, classical Bessel series expansions of various special functions which are connected with the incomplete gamma function, can be easily obtained in a unified form. A Neumann series of Bessel functions is an expansion of the type 1 X (1.4) an Jν+n(z)(z; ν C); n=0 2 1 where an n=0 are given coefficients. Many special functions of the mathematical physicsf enjoyg an expansion of this type, e.g., Kummer confluent hypergeometric's, Lommel's , Kelvin's, Whittaker's, and so on. In Section 5 we exploit the novel representation of Bessel's functions presented in Section 3 in order to derive a sim- ple closed-form integral representation of expansions (1.4). The set of coefficients 1 an n=0, which characterizes the series (1.4), comes into this integral representation inf ag very simple way through the associated analytic function A(z) in the unit disk P1 n D, i.e.: A(z) = n=0 anz with a weak restriction on the sequence of coefficients an , which is sufficient to guarantee the boundedness of A(z) on the unit circle. Examplesf g of special functions with this kind of closed-form integral representation are also given. 2. Connection between the coefficients of Gegenbauer and Fourier expansions ν 1 The Gegenbauer (ultraspherical) polynomials C` (x) of order ν (ν > 2 ; ν = 0) may be defined by means of the generating function [23, Eq. (4.7.23)] (see− also6 [23, Eq. (4.7.6)]): 1 X (2.1) Cν (x) t` = (1 2xt + t2)−ν (x [ 1; 1]; t ( 1; 1)): ` − 2 − 2 − `=0 For fixed ν, the Gegenbauer polynomials are orthogonal on the interval [ 1; 1] with respect to the weight function wν (x) = (1 x2)(ν−1=2), that is: − − Z 1 ν ν ν ν (2.2) C` (x) Cn(x) w (x) dx = h` δ`;n; −1 where: 1−2ν ν 1 2 π Γ(` + 2ν) (2.3) h` = : (` + ν) `! [Γ(ν)]2 Particularly important special cases of Gegenbauer polynomials are obtained for ν = 1=2, which gives the Legendre polynomials P`(x), for ν = 1 yielding the Cheby- shev polynomials of the second kind U`(x), and in a suitable limiting form for ν 0 ν ! the Chebyshev polynomials of the first kind T`(x) [23]: limν!0(`!=(2ν)`)C` (x) = T`(x), where ( )` denotes the Pochhammer symbol. · BESSEL FUNCTIONS OF THE FIRST KIND 3 eiu δ (u) + δ (u) − δ0(u) 1 0 1 − iu e− Figure 1. Integration path for the evaluation of the integral representation (2.5) of the Gegenbauer polynomials. For Re ν > 0 they admit the following integral representation [26, Eq. (1), p. 559]: ν Z π ν (` + ν) h` p 2 ` 2ν−1 (2.4) C` (x) = (x + i 1 x cos η) (sin η) dη: π 0 − Now, our aim in this section is to present an integral representation of the Gegen- bauer polynomials. To this end, we prove the following proposition (see also [6, 9]). Proposition 2.1. For Re ν > 0, the following integral representation for the ν Gegenbauer polynomials C` (cos u) (u [0; π], ` N0) holds: 2 2 (` + ν) hν 2ν−1e−iπν Z 2π−u (2.5) Cν (cos u) = ` ei(`+ν)t (cos u cos t)ν−1 dt: ` (2ν−1) π (sin u) u − Proof. In representation (2.4) we first introduce the variable u, defined by x = cos u (u [0; π]): 2 ν Z π ν (` + ν) h` ` 2ν−1 (2.6) C` (cos u) = (cos u + i sin u cos η) (sin η) dη: π 0 Then, in the integral in (2.6) we substitute to η the complex integration variable τ defined by eiτ = cos u + i sin u cos η: It can be checked that (2.7) 2eiτ (cos τ cos u) = (eiτ eiu)(eiτ e−iu) = sin2u sin2η: − − − 4 E. DE MICHELI Now, since eiτ dτ = sin u sin η dη, the integrand on the right-hand side of (2.6) can be written as follows:− (ν−1) (sin u)−(2ν−1) ei(`+1)τ (eiτ eiu)(eiτ e−iu) dτ (2.8) − − − = (sin u)−(2ν−1) ei(`+ν)τ [2(cos τ cos u)](ν−1) dτ: − − In order to determine the integration path, consider the intermediate step where eiτ is chosen as integration variable; the original path (corresponding to η [0; π]) iu −iu 2 is the (oriented) linear segment δ0(u) starting at e and ending at e . Since (as shown by (2.8)) the integrand is an analytic function of eiτ in the disk eiτ < 1 j j (since ` N0), the integration path δ0(u) can be replaced by the circular path 2 iτ δ+(u) = e : τ = t; u 6 t 6 2π u (see Fig. 1). Moreover, by using the fact that iτ iuf iτ −iu ν−1 − g iτ iτ iπ [(e e )(e e )] is positive for e δ0(u) R and therefore at e = e , we conclude− from− the left equality in (2.7) that2 in the[ r.h.s. of (2.8) the following specification holds (for τ = t; u t 2π u): 6 6 − [2(cos t cos u)]ν−1 = e−i(ν−1)π[2(cos u cos t)]ν−1: − − Finally, accounting for this latter expression, the integral representation (2.5) fol- lows directly from (2.6). Remark 1. In the proof of Proposition 2.1 we could even choose the integration iτ path δ−(u) = e : τ = t; u 6 t 6 u , oriented as shown in Fig. 1. Following the same argumentsf given above,− it is easilyg seen that this option yields the following alternative integral representation of the Gegenbauer polynomials [6, 8]: ν ν−1 Z u ν (` + ν) h` 2 i(`+ν)t ν−1 (2.9) C` (cos u) = 2ν−1 e (cos t cos u) dt: π (sin u) −u − 3. Fourier-type integral representation of Bessel functions of the first kind It is well-known that the Bessel functions of the first kind of integral order n can be represented by Bessel's trigonometric integral (1.3) (see also [11, Eq. n 10.9.2]). For n Z, Jn(x) satisfies the symmetry relation: J−n(x) = ( 1) Jn(x). This latter relation2 no longer holds when n is not an integral number− since, for ν Z, Jν (x) and J−ν (x) are linearly independent solutions of the Bessel equation of62 order ν.