The Bessel Functions - Brown University

1 April 9, 201318:41 World Scientific Book - 9in x 6inMAINARDI BOOK-FINALA ppendix BThe Bessel FunctionsAs Rainville pointed out in his classic booklet[Rainville (1960)],no other special Functions have received such detailed treatment inreadily available treatises as the Bessel Functions . Consequently, wehere present only a brief introduction to the subject including therelated Laplace transform pairs used in this The standard Bessel functionsThe Bessel Functions of the first and second kind:J ,Y.

2 The Bessel Functions of the first kindJ (z) are defined from theirpower series representation:J (z) := k=0( 1)k (k+ 1) (k+ + 1)(z2)2k+ ,( )wherezis a complex variable and is a parameter which can takearbitrary real or complex values. When is integer it turns out asan entire function; in this caseJ n(z) = ( 1)nJn(z), n= 1,2,..( )In factJn(z) = k=0( 1)kk!(k+n)!(z2)2k+n,J n(z) = k=n( 1)kk!(k n)!(z2)2k n= s=0( 1)n+s(n+s)!s!(z2)2s+ 9, 201318:41 World Scientific Book - 9in x 6inMAINARDI BOOK-FINAL174 Fractional Calculus and Waves in linear ViscoelasticyWhen is not integer the Bessel Functions exhibit a branch pointatz= 0 because of the factor (z/2) , sozis intended with|arg(z)|< that is in the complex plane cut along the negative real a suggestion by Tricomi, see[Gatteschi (1973)], we canextract from the series in ( ) that singular factor and set:JT (z) := (z/2) J (z) = k=0( 1)kk!

3 (k+ + 1)(z2)2k.( )The entire functionJT (z) was referred to by Tricomi as theuniformBessel function. In some textbooks on special Functions , see [Kiryakova (1994)], p. 336, the related entire functionJC (z) :=z /2J (2z1/2) = k=0( 1)kzkk! (k+ + 1)( )is introduced and named theBessel-Clifford for fixedzin the cut plane the terms of the series ( )are analytic function of the variable , the fact that the series isuniformly convergent implies that the Bessel function of the firstkindJ (z) is an entire function of order.

4 The Bessel Functions are usually introduced in the framework ofthe Fucks Frobenius theory of the second order differential equationsof the formd2dz2u(z) +p(z)ddzu(z) +q(z)u(z) = 0,( )wherep(z) andq(z) are assigned analytic Functions . If we chose in( )p(z) =1z, q(z) = 1 2z2,( )and solve by power series, we would just obtain the series in ( ).As a consequence, we say that the Bessel function of the first kindsatisfies the equationu (z) +1zu (z) +(1 2z2)u(z) = 0,( )where, for shortness we have used the apices to denote differentiationwith respect toz.

5 It is customary to refer to Eq. ( ) as theBesseldifferential 9, 201318:41 World Scientific Book - 9in x 6inMAINARDI BOOK-FINALA ppendix B: The Bessel Functions175 When is not integer the general integral of the Bessel equationisu(z) = 1J (z) + 2J (z), 1, 2 C,( )sinceJ (z) andJ (z) are in this case linearly independent withWronskianW{J (z),J (z)}= 2 zsin( ).( )We have used the notationW{f(z),g(z)}:=f(z)g (z) f (z)g(z).In order to get a solution of Eq. ( ) that is linearly independentfromJ also when =n(n= 0, 1, ) we introduce the Besselfunction of the second kindY (z) :=J (z) cos( ) J (z)sin( ).

6 ( )For integer the of ( ) becomes indeterminate so in thiscase we defineYn(z) as the limitYn(z) := lim nY (z) =1 [ J (z) =n ( 1)n J (z) =n].( )We also note that ( ) impliesY n(z) = ( 1)nYn(z).( )Then, when is an arbitrary real number, the general integral of Eq.( ) isu(z) = 1J (z) + 2Y (z), 1, 2 C,( )and the corresponding Wronskian turns out to beW{J (z),Y (z)}=2 z.( )The Bessel Functions of the third kind:H(1) ,H(2) .In ad-dition to the Bessel Functions of the first and second kind it is cus-tomary to consider the Bessel function of the third kind, or Hankelfunctions, defined asH(1) (z) :=J (z) +iY (z), H(2) (z) :=J (z) iY (z).

7 ( )These Functions turn to be linearly independent with WronskianW{H(1) (z),H(2) (z)}= 4i z.( )April 9, 201318:41 World Scientific Book - 9in x 6inMAINARDI BOOK-FINAL176 Fractional Calculus and Waves in linear ViscoelasticyUsing ( ) to eliminateYn(z) from ( ), we obtain H(1) (z) :=J (z) e i J (z)isin( ),H(2) (z) :=e+i J (z) J (z)isin( ),( )which imply the important formulasH(1) (z) = e+i H(1) (z), H(2) (z) = e i H(2) (z).( )The recurrence relations for the Bessel func-tionsJ (z),Y (z),H(1) (z),H(2) (z) satisfy simplerecurrence rela-tions.

8 Denoting any one of them byC (z) we have: C (z) =z2 [C 1(z) +C +1(z)],C (z) =12[C 1(z) C +1(z)].( )In particular we noteJ 0(z) = J1(z), Y 0(z) = Y1(z).We note thatC stands forcylinder function, as it is usual to call thedifferent kinds of Bessel Functions . The origin of the termcylinder isdue to the fact that these Functions are encountered in studying theboundary value problems of potential theory for cylindrical more general differential equation for the Bessel differential equation ( ) can be generalized by intro-ducing three additional complex parameters ,p,qin such a wayz2w (z) + (1 2p)zw (z) +( 2q2z2q+p2 2q2)w(z) = 0.

9 ( )A particular integral of this equation is provided byw(z) =zpC ( zq).( )We see that for = 1,p= 0,q= 1 we recover Eq. ( ).April 9, 201318:41 World Scientific Book - 9in x 6inMAINARDI BOOK-FINALA ppendix B: The Bessel Functions177 The asymptotic representations for the Bessel asymptotic representations of the standard Bessel Functions forz 0 andz are provided by the first term of the convergentseries expansion aroundz= 0 and by the first term of the asymptoticseries expansion forz , 0 (with|arg(z)|< if is not integer) we have: J n(z) ( 1)n(z/2)nn!

10 , n= 0,1,..,J (z) (z/2) ( + 1), 6= 1, ( ) Y0(z) iH(1)0(z) iH(2)0(z) 2 log (z),Y (z) iH(1) (z) iH(2) (z) 1 ( )(z/2) , >0.( )Forz with|arg(z)|< and for any we have: J (z) 2 zcos(z 2 4),Y (z) 2 zsin(z 2 4),H(1) (z) 2 ze+i(z 2 4),H(2) (z) 2 ze i(z 2 4).( )The generating function of the Bessel Functions of Bessel Functions of the first kindJn(z) are simply re-lated to the coefficients of the Laurent expansion of the functionw(z,t) = ez(t 1/t)/2=+ n= cn(z)tn,0<|t|< .( )To this aim we multiply the power series of ezt/2, e z/(2t), and, aftersome manipulation, we getw(z,t) = ez(t 1/t)/2=+ n= Jn(z)tn,0<|t|<.