arxiv:math/ v1 [math.ca] 9 Jul 1993

Transcription

1 Some Integrals Involving Bessel Functions arxiv:math/9373v [math.ca] 9 Jul 993 M.L. Glasser Department of Physics and Department of Mathematics and Computer Science, Clarkson University, Potsdam, N.Y (USA) E. Montaldi Dipartimento di Fisica, Universita di Milano, Via Celoria 6, 33 Milano, Italy Abstract A number of new definite integrals involving Bessel functions are presented. These have been derived by finding new integral representations for the product of two Bessel functions of different order and argument in terms of the generalized hypergeometric function with subsequent reduction to special cases. Connection is made with Weber s second exponential integral and Laplace transforms of products of three Bessel functions. Keywords: Bessel Function, Hypergeometric Function, Integral Representation AMS Classification No. 33A4, 33A3

2 Some Integrals Involving Bessel Functions M.L. Glasser* and E. Montaldi** *Department of Physics and Department of Mathematics and Computer Science, Clarkson University, Potsdam, N.Y (USA) **Dipartimento di Fisica, Universita di Milano, Via Celoria 6, 33 Milano, Italy. Introduction The aim of this work is to derive a number of infinite integrals involving Bessel functions which appear to be new. The approach is, beginning with an expression for the product of two Bessel functions as a sum of Gauss functions, to integrate and perform a resummation to obtain other hypergeometric functions, and then to reduce these to more familiar form. In particular we exploit the relation between Bessel functions and the function F 3 []. Finally, we generalize Weber s second exponential integral by expressing the Laplace transform of a product of three Bessel functions as an infinite series of products of modified Bessel functions. We have not aimed at complete rigor or generality; operations such as interchange of limits and Hankel inversion are carried out formally, i.e. without verifying that the conditions stated ensure their validity. In the majority of cases approriate conditions can be supplied by appeal to convergence and analytic continuation.. On the product J µ (ax)j ν (bx) Our starting point is the familiar expansion [] Γ(µ + )Γ(ν + )J µ (ax)j ν (bx) = ( ax)µ ( bx)ν m= ( ) m ( ax)m m!(µ + ) m F ( m, µ m; ν + ; b a) (.)

3 By making use of the standard transformation[3], we have F ( m, µ m; ν + ; b a ) = = ( b F a )µ+ν+m+ (ν + m +, µ + ν + m + ; ν + ; b a) (.) Next, we recall that [4] F (ν + m +, ν + µ + m + ; ν + ; b a ) = = 4(a/b) ν Γ(µ + ν + )(ν + ) m (µ + ν + ) m t µ+ν+m+ K µ (t)i ν ( b t)dt (.3) a m =,,,..., Re ν >, Re (µ + ν) >, Re( b ) <. a By inserting (.) and (.3) into (.), we get Γ(ν + )Γ(µ + )Γ(µ + ν + )J µ (ax)j ν (bx) = 4( ax)µ ( bx)ν ( b a )µ+ν+ ( a b )ν t µ+ν+ F 3 (µ +, ν +, µ + ν + ; z t )K µ (t)i ν ( b t)dt (.4) a Re ν >, Re (µ + ν) >, Re( b ) <, a where z = b ax( ). a It is interesting to observe that eq. (.4) enables us to derive the Mellin transform of J µ (ax)j ν (bx), i.e. the Weber-Schafheitlin integral [5], in a simple way. Indeed from [6], x µ+ν s F 3 (µ +, ν +, µ + ν + ; z t )dx = = [ b a( and [4] )Γ(µ + )Γ(ν + )Γ(µ + ν + ), (.5) )t]s µ ν Γ(µ+ν s+ a Γ( µ ν+s+ )Γ( ν µ+s+ )Γ( µ+ν+s+ ) Re (µ + ν s) >, 3

4 t s K µ (t)i ν ( b a t)dt = we immediately obtain F ( ν + µ + s + (b/a)ν 4Γ(ν + ) Γ(µ + ν + s + )Γ( ν µ + s + ), ν µ + s + Re(ν ± µ + s) >, Re( b a ) < x s J µ (ax)j ν (bx)dx = s b ν a s ν ; ν + ; b a) (.6) Γ( µ+ν s+ ) Γ(ν + )Γ( µ ν+s+ ) ( b F a )s ( µ + ν + s +, ν µ + s + ; ν + ; b a ) = s b ν a s ν Γ( µ+ν s+ ) Γ(ν + )Γ( µ ν+s+ ) F ( ν µ s + Re(µ + ν s) >, < b < a With b ib (b > ),and x positive, eq. (.4) becomes Γ(µ + )Γ(ν + )Γ(µ + ν + )J µ (ax)i ν (bx) = ( 4 ax)µ ( 4 bx)ν ( + b a )µ+ν+ ( a b )ν, ν + µ s + ; ν+; b a ) (.7) t µ+ν+ F 3 (µ+, ν+, µ+ν+; 6 a x (+ b t )K a ) µ (t)j ν ( b t)dt (.8) a Re ν >, Re(µ + ν) > or, by writing b=ay and x=a/4 Γ(µ + )Γ(ν + )Γ(µ + ν + )J µ ( 4 a )I ν ( 4 a y) = ( a 6 )µ+ν ( + y ) µ+ν+ t µ+ν+ F 3 (µ +, ν +, µ + ν + ; a4 56 ( + y ) t )K µ (t)j ν (yt)dt (.9) Re ν >, Re(µ + ν) >. Also, with a 6a +y, 4

5 4a Γ(µ + )Γ(ν + )Γ(µ + ν + )J µ ( + y )I ν( 4ay + y ) = ( + y )a µ+ν t µ+ν+ F 3 (µ +, ν +, µ + ν + ; a t )K µ (t)j ν (yt)dt (.) Re ν >, Re(µ + ν) >, a, y >. For Reν / we may use the Hankel inversion formula [7], and eq. (.) gives (at) µ+ν F 3 (µ+, ν+, µ+ν+; a t )K µ (t) = Γ(µ+)Γ(ν+)Γ(µ+ν+) y 4a + y J ν(ty)j µ ( + y )I ν( 4ay + y)dy (.) This is an interesting addition to the class of Sonine-Gegenbauer integrals [8]. In particular, by dividing by a µ+ν and taking the limit a, which may be taken under the integral sign, we have ( t y ν+ )µ+ν K µ (t) = Γ(µ + ν + ) ( + y ) µ+ν+j ν(ty)dy (.) The four particular cases µ =, ν = ; µ =, ν = ; ν =, µ = ; ν =, µ = of eq.(.9) are of some interest. By using and x F 3 (, x4, ; ) = ber(x) (.3) 56 4 F 3 ( 3, 3, : x4 ) = bei(x) (.4) 56 we obtain the four integral representations π ( + y ) / J ( a 4 ) cosh(a y 4 ) = 5 K (t)ber(a ( + y )t) cos(yt)dt (.5)

6 π (+y ) / J ( a 4 ) sinh(a y ( + y ) / I ( a y 4 ) cos(a 4 ) = 4 ) = ( + y ) / I ( a y 4 ) sin(a 4 ) = Therefore, K (t)ber(a t) = K (t)bei(a t) = t e t ber(a t) = t e t bei(a t) = K (t)bei(a ( + y )t) sin(yt)dt (.6) e t ber(a ( + y )t)j (yt)dt (.7) e t bei(a ( + y )t)j (yt)dt (.8). ( + y ) / J ( a 4 + y ) cosh( a y 4 + y) cos(ty)dy (.9) (+y ) / J ( a 4 + y ) sinh( a y 4 + y) sin(ty)dy (.) y( + y ) / I ( a y 4 + y ) cos( a 4 + y )J (ty)dy (.) y( + y ) / I ( a y 4 + y ) sin( a 4 + y )J (ty)dy (.) A curious formula arising from (.4) may be mentioned. We set µ =, replace I ν ( b t) by Poisson s integral [9] a I ν ( b a t) = Γ(ν + )(bt a )ν e b a tu ( u ) ν du (.3) Re ν > interchange the order of integration, and then use the Laplace transform e βt t ν+ F 3 ( 3, ν +, ν + 3 ; Γ(ν + ) α αt )dt = 4 α β (ν+) sin(4 β ) (.4) Re ν >, Re β > which is readily established by expressing the F 3 as its power series and integrating term-by-term. This leads to sin(ax)j ν (bx) = 6

7 = πγ(ν + )( bx)ν (a b ) ν+ Re ν >, a > b >. ( u ) ν b (a + bu) ν+ sin(a x)du (.5) a + bu Therefore, with x = π a and b = ay J ν ( πy ) = πγ(ν + )( 4 πy)ν ( y ) ν+ ( u ) ν / ( + uy) sin(π y )du (.6) ν+ + uy Re ν >, y <. By taking advantage of other known reductions for the F 3 functions in (.4) many additional new integrals can be derived. For example from [3] we obtain e x I ν (x sin θ)[sin( 3πν )ber ν( cos θ ux)+cos( 3πν )bei ν( cos θ ux)]dx and in particular = sec θ sin u J ν (u sin θ), e x I n (x sin θ)bei 4n ( cos θ ux)dx = ( ) n sec θ sin u J n (u sin θ) e x I n+ (x sin θ)ber 4n+ ( cos θ ux)dx = ( ) n+ sec θ sin u J n+ (u sin θ). We conclude this section by deriving a further integral representation for the product J ν (ax)j ν (bx) in terms of a F 3, but different from (.4). Let us consider eq.(.) with µ = ν. According to the quadratic transformation [] F ( m, ν m; ν +; b b a) = (+ F a )m ( m, m ab ; ν +; ( ) = a + b ) = ( + b m! a )m [ m ] r= ab ( (.7) (m r)!r!(ν + ) r a + b )r 7

8 and (.) becomes By using [] it follows that r= J ν (ax)j ν (bx) = ( 4 abx ) ν [ m ] r= m= ( ) m ( x a + b ) m Γ(ν + m + ) ( ab ) r a +b (m r)!r!γ(ν + r + ). (.8) [ m ] c(m, r) = c(m + r, r) (.9) m= r= m,r= J ν (ax)j ν (bx) = r!γ(ν + r + ) ( abx a + b )ν+r m= ( ) m ( x a + b ) ν+r+m m!γ(ν + r + m + ) = ( abx abx ( ) r a +b a + b )ν r!γ(ν + r + ) J ν+r(x a + b ) (.3) r= We note, in passing, that (.3) provides a quick derivation of Weber s second integral []. Indeed, by using we have e px x ν+r+ J ν+r (x a + b )dx = (a + b ) ν/+r xe px J ν (ax)j ν (bx)dx = a +b p e 4p (p) ν+r+ r= = e a +b 4p (.3) (ab/4p) ν+r r!γ(ν + r + ) = = a +b p e 4p I ν ( ab p ) (.3) We also observe that (.3) can be obtained from the formula [3] = r= F (a, b; c; x) F (a, b; c; y) = (a) r (b) r (c a) r (c b) r r!(c) r (c) r (xy) r F (a+r, b+r; c+r; x+y xy) (.33) 8

9 by applying the confluence principle [4] twice. We first replace x by x/b, y by y/b, and let b. This gives F (a; c; x) F (a; c; y) = r= (a) r (c a) r r!(c) r (c) r ( xy) r F (a + r; c + r; x + y) Next, we replace x by x/a, y by y/a and let a ; then F (c; x) F (c; x) = r= (.34). (xy) r r!(c) r (c) r F (c + r; x + y) (.35) which is equivalent to the desired result. If on the right hand side of (.3) we write [9] J ν+r (x a + b ) = πγ(ν + we get (Re ν > ) J ν (ax)j ν (bx) = ((/)x a + b ) ν+r ) ( ν + ) 4 r( ν + 3 ) 4 r r ( t ) ν +r cos(xt a + b )dt, πγ(ν + ) (abx ) ν ( t ) ν cos(xt a + b ) (.36) F 3 (ν +, ν + 4, ν ; 64 a b x 4 ( t ) )dt (.37) and with b =, x = a /, and a = 4 ( u + + u ), u + + u J ν ( u + u )J ν ( ) = ( t ) ν / cos(ut) πγ(ν + ) F 3 (ν +, ν + 4, ν ; 64 ( t ) )dt (.38) Finally, with ν = / (.37) yields [3] a + b cos(ut )[I (u t ab ) + J (u dt t )] t 9

10 = (/u) sin(u a/b) sin(u b/a). Further, more complex, evaluations are possible by the same procedure. 3. The integral e αx J (β x)j (β x)j (β 3 x)dx. The integral in this section heading, which we denote I(β, β, β 3 ) is an extension of the ν = case of Weber s second integral, eq. (.3), to which it reduces when one of the three parameters vanishes. By writing (see eq.(.)) and [5] ( ) n ( J (β x)j (β β 3 x) = 3 x) n F n! ( n, n; ; β ) (3.) β3 n= e αx x n n! J (β x)dx = /4α L n ( β ), (3.) α n+e β 4α L n (x) being a Laguerre polynomial, we first have I(β, β, β 3 ) = α e β Now [6] 4α F ( n, n; ; β β 3 ( ) n ( β 3 n= n! 4α )n F ( n, n; ; β β3 ) = ( β ) n F β3 ( n, n + ; ; ) β β3 )L n ( β 4α ) (3.3) = ( β ) n P β3 n ( β 3 + β ) = π (β β3 β πβ3 n + β 3 β β 3 cos θ) n dθ. (3.4) Therefore, in terms of a well known generating function[7], (3.3) becomes β I(β, β, β 3 ) = e 4α πα π (/4α) n (β + β3 β β 3 cos θ) n L n ( β n! 4α )dθ n= = e β /4α πα π β e R /4α I ( β R )dθ (3.5) α where R = (β +β 3 β β 3 cos θ) /. Next, by using Graf s addition theorem [8] I ( β R α ) = ( δ no )I n ( β β α )I n( β β ) cos(nθ) (3.6) α n=

11 and the familiar integral representation we finally have π e β β α cos θ cos(nθ)dθ = πi n ( β β ), (3.7) α where = α e e αx J (β x)j (β x)j (β 3 x)dx 4α (β +β +β 3 ) n= ( δ n )F n (β, β, β 3 ) (3.8) F n (β, β, β 3 ) = I n ( β β α )I n( β β 3 α )I n( β β 3 α ) (3.9) We conclude by sketching the evaluation of the more general integral I m (β, β, β 3 ) = e αx J (β x)jm (β x)jm (β 3 x)dx (3.) where m is a positive integer. In the first place, (3.4) takes the form I m (β, β, β 3 ) = m!α (β β 3 4α )m e β ( ) n ( β 3 n= n! 4α )n F ( n, n m; m + ; β β3 By observing that [9] F ( n, n m; m + ; β β 3 eq. (3.) becomes 4α ) = ( β ) n n!(m)! β3 (m + n)! Cm+/ )L m+n ( β 4α ) (3.) n ( β 3 + β ) β3 β = m m! π sin m θ (β (m)!πβ3 n + β3 β β 3 cos θ) n dθ (3.) I m (β, β, β 3 ) = m! πα (m)! (β β 3 α )m e β /4α π sin m θ ( /4α) n n= n!

12 Now, from [], (β + β3 β β 3 cos θ) n L m+n ( β )dθ. (3.3) 4α L m+n (x) = n! (m + n)! ex ( d dx )m e x x m L m n (x) (3.4) and [7] n= t n (m + n)! Lm n (x) = (xt) m/ e t J m ( xt) (3.5) it follows that (where R has the same meaning as before) n= n! ( R 4α )n L m+n ( β 4α ) = e(β R )/4α m ( d dx )m e x x m By inserting (3.6) into (3.3) and using [] and [] (R x α ) m I m (R x α ) x=β /4α. (3.6) x x (R α ) m I m (R α ) = (β β 3 x α ) m (m )! we finally obtain n= π (m + n)cn m (cos θ) n= I m+n (β x α )I m+n(β 3 x α ) (3.8) sin m θ e β β 3 α cos θ C m n (cos θ)dθ = π m (m + n )! ( α ) m I m+n ( β β 3 ), (3.8) n!(m )! β β 3 α = α e αx J (β x)jm (β x)jm (β 3 x)dx m+ m! ( α ) m e β +β 3 d 4α ( (m)! β β 3 dx )m e x (3.9) (n + m)(m + n )! x x I m+n (β n! α )I m+n(β 3 α )I m+n( β β 3 α ) x=β /4α.

13 In particular, by letting β 3, after dividing by β m 3, e αx J (β x)jm (β x)x m/ dx = α e β /4α ( d dx )m e x ( x x α )m/ I m (β α ) x=β /4α (3.) = α ( β ) α )m e β +β m 4α n=( ) n( m ( β ) n I n ( β β n β α ). which is of interest in connection with formulas (39)-(4) on page 86 of reference [5]. In conclusion, we point out that the derivation of (.4) can be extended to give e βx F 3 (µ, ν, ν + ; a x )dx = (a) µ Γ(µ)Γ(ν)β µ ν J µ (4a/β) (3.) Therefore, the results of section, for example, are capable of extension in a variety of directions. We leave this for the future and merely quote one example, interesting because it contains each of the four types of Bessel functions: xj (ax)i (ax)y (x)k (x) dx = (πa ) ln( a 4 ) (3.) where < a <. References [] T.J. Lardner, Siam Review, 69 (969). [] Bateman Manuscript Project, Higher Transcendental Functions (McGraw- Hill, 953), Vol., p., eq. (47). [3] Ibid. Vol., p.64, eq. (3). [4] Ibid. Vol., p.5, eq. (3). [5] Ibid. Vol., p.5, eq. (9). [6] L.J. Slater, Generalized Hypergeometric Functions (Cambridge 966), Sec [7] I.N. Sneddon, The Use of Integral Transforms (McGraw-Hill 97), Chap. 5, p.39. 3

14 [8] See Ref.[] Vol., p.96, eq. (59). [9] Ibid. Vol., p.8, eq. (). [] See Ref. [], p., eq. (). [] H.M. Srivastava and H.L. Manocha, A treatise on Generating Functions (John Wiley 983), p., eq. (6). [] See Ref. [], Vol., p.5, eq. (5). [3] J.L. Burchnall and T.W. Chaundy, Quart. J. Math. Oxford Ser., 49-7 (94). [4] See Ref. [], p.36. [5] Bateman Manuscript Project, Tables of Integral Transforms (McGraw- Hill 953), Vol., p. 3, eq. (3). [6] See Ref. [], p. 57, eq. (5). [7] See Ref. [], p. 89, eq. (8). [8] Ibid. Vol., p.44, eq.(5). [9] Ibid. Vol., p.77, eq.(3). [] Ibid. Vol., p.9, eq.(8). [] Ibid. Vol., p., eq.(3). [] See Ref. [5], p.8, eq.(7). [3] A.P. Prudnikov, Ya. A. Brychkov and O.I. Marichev, Integrals and Series (Gordon and Breach, N.Y. 989) Vol. 3 Sec Eq. (6). [4] Ibid. Sec Eq.(5). 4