Bessel Functions Michael Taylor. Lecture Notes for Math 524

Bessel Functions Michael Taylor Lecture Notes for Math 54 Contents 1. Introduction. Conversion to first order systems 3. The Bessel functions J ν 4. The Bessel functions Y ν 5. Relations between J ν and J ν±1 6. The functions J n+1/ 7. Integral formula for J ν 8. The Hankel functions H ν (k) 9. Large t behavior 1. Zeros of J ν and Y ν A. Where Bessel functions come from B. The Euler gamma function Γ(z) C. Derivative of 1/Γ(z) at integer points 1. Introduction Here we consider solutions to the Bessel equation (1.1) u (t) + 1 t u (t) + (1 ν t ) u(t) =, for ν R. This equation arises in the study of partial differential equations on regions with special shapes in n-dimensional Euclidean space, and in cases n = and n = 3 such PDE are very important in mathematical physics. See Appendix A of these notes for material on this. The equation (1.1) has a regular singular point at t =. In we convert (1.1) to a first order system, of the form (1.) t dx dt = A(t)x, with a regular singular point at t =. This allows us to take advantage of material developed in 11, Chapter 3 of [T1], though our presentation in is actually self contained. 1

In 3 we study a certain Bessel function J ν (t) solving (1.1). Results of guarantee that J ν (t) = t ν J ν (t) has a convergent power series a kt k, and we derive a recursion formula for the coefficients a k. We produce a solution to this recursion, and hence define J ν (t). The solution involves the gamma function Γ(z), and we make use of results on Γ(z) given in Appendix B. We see that, for ν R, both J ν (t) and J ν (t) solve (1.1). If ν is not an integer, we show that these solutions are linearly independent, and hence form a basis for solutions to (1.1) on t (, ). If ν = n is an integer, we show J n (t) = ( 1) n J n (t). Section 4 produces a solution Y ν (t) to (1.1), having the property that the Wronskian (1.3) W (J ν, Y ν )(t) = πt, for all ν R, t (, ), so {J ν, Y ν } forms a basis of solutions to (1.1) for each ν R. In 5 we give some identities connecting J ν (t) with J ν±1 (t). In 6 we discuss J n+1/ (t), which is an elementary function when n is an integer. In 7 we establish an integral formula for J ν (t). In 8 we introduce the Hankel functions H ν (1) and H ν (). In 9 we produce results on the asymptotic behavior of J ν (t), H ν (1) (t), and Y ν (t) as t +. In 1 we apply results of 9 to show that, for each ν R, real valued solutions to (1.1) have infinitely many zeros (say, ϑ kν ) on (, ), and give their approximate values, for large k. Material in 8 9 is less self-contained than in other sections. We refer to [W] for a derivation of the integral formula (8.1) for H ν (1) (t), and to Chapter 3 of [T3] for the methods of asymptotic analysis that apply to the integral formula (7.1) for J ν (t). As already indicated, Appendix A explains why Bessel functions are so important for certain PDEs, and Appendix B introduces the gamma function, which is a useful tool in the analysis of Bessel functions. Appendix C evaluates the derivative of 1/Γ(z) for integer z, of use in the formula for Y n (t) derived in 4.

3. Conversion to first order systems To begin our treatment of Bessel s equation (.1) u (t) + 1 t u (t) + we convert it to a first-order system for (.) x(t) = (1 ν t ) u(t) =, ( ) u(t), v(t) = tu (t), v(t) obtaining (.3) ( ) dx u dt = u + tu, hence (.4) t dx dt = This has the form ( ) 1 ν t x. (.5) t dx dt = A(t)x, A(t) = A + A t, with (.6) A = ( ) 1 ν, A = ( ). 1 Generally, (.5) has a solution given as a convergent power series (.7) x(t) = x + x 1 t + x t +, as long as A x =, provided also A has no eigenvalues that are positive integers. This is established in Lemma 11.1, Chapter 3, of [T1]. (A proof that works in the cases we consider here is given in (.16) (.4) below.) In case (.6), A has a nonzero null space only for ν =. To treat the more general cases, we set (.8) x(t) = t ν y(t).

4 This is motivated by the model Euler equation (.9) t dx dt = A x, whose solution for t > is (.1) x(t) = e (log t)a = t A, upon noting that A in (1.6) has eigenvalues ±ν. From (.8) we have x = t ν y + νt ν 1 y, so (.11) A(t)x = t dx dt = tν t dy dt + νtν y, so (.5) is equivalent to (.1) t dy dt = ( A(t) νi ) y(t). In case (.5) (.6), we have ( (.13) A(t) νi = (A νi) + A t ν 1, A νi = ν ν Since ( ) 1 (.14) (A νi) = ν ( ), Lemma 11.1 (Chapter 3) of [T1] (which, we recall, is proved in the case needed here below) implies (.1) has a solution given by a convergent power series (.15) y(t) = y + y 1 t + y t +, y = ( ) 1, ν provided A νi has no positive integers as eigenvalues. Since Spec(A νi) = {, ν}, this works as long as ν is not a negative integer. Actually, the fact that only even powers of t appear on the right side of (.13) allows us to improve this a bit. Let us rewrite (.1) (.13) as (.16) t dy dt = B(t )y, with (.17) B(s) = B + B 1 s, B = A νi, B 1 = A. ).

5 Set (.18) y(t) = z(t ), t >. Thus ty (t) = t z (t ) = B(t )y(t), so (.19) s dz ds = 1 B(s)z(s), and since (.) 1 B = 1 ( ) ν 1 ν, Spec 1 ν B = {, ν}, Lemma 11.1 mentioned above applies, to yield a solution to (.19) (.1) z(s) = z k s k, z = ( ) 1, ν as long as ν is not a negative integer. In fact, plugging this power series into (.19) gives the recursion (.) kz k = 1 B z k + 1 B 1z k 1, i.e., (.3) z k = (ki B ) 1 B 1 z k 1, starting with z = (1, ν) t, and then inductively specifying z k for k {1,, 3,... }, as long as k / Spec B for each such k, i.e., as long as Spec B contains no positive even integer, hence as long as ν is not a negative integer. Note that (.3) implies (.4) z k C k B 1 z k 1, for k sufficiently large, and this implies convergence of (.1) for all s R. This argument proves the special case of Lemma 11.1 (Chapter 3) of [T1] that we have been citing.

6 3. The Bessel functions J ν Here we construct a solution to (1.1), for t >, called the Bessel function J ν (t). In this case, (.) becomes (3.1) x(t) = and (.8) becomes (3.) y(t) = Thus we are motivated to set ( ) Jν (t) tj ν(t), ( t ν ) J ν (t) t 1 ν J ν(t). (3.3) J ν (t) = t ν J ν (t), so J ν (t) = t ν J ν (t), hence (3.4) J ν(t) = t ν J ν(t) + νt ν 1 J ν (t), so ( ) J ν (t) (3.5) y(t) = tj ν(t). + νj ν (t) The equation (.1), i.e., (3.6) t dy dt = ( A(t) νi ) y(t), becomes (3.7) t d ( ) J ν dt tj ν = + νj ν [ ( ) ν 1 ν + ν Comparing the first components, we get ( ) 1 tj ν = νj ν + tj ν + νj ν, J ν t ]( tj ν + νj ν which is a tautology, and comparing the second components, we get (3.8) t ( tj ν + (ν + 1)J ν) = ν J ν ν(tj ν + νj ν ) t J ν, ).

7 which we rewrite as (3.9) J ν (t) + ν + 1 J t ν(t) + J ν (t) =. The material developed in (.16) (.4) applies to (3.7), and guarantees that (3.9) has a solution given as a convergent power series of the form (3.1) J ν (t) = a k t k, t R, as long as ν is not a negative integer. Here we give an iterative construction of such coefficients. Given (3.1), the left side of (3.9) is (3.11) } {(k + )(k + ν + )a k+1 + a k t k. As long as ν / { 1,, 3,... }, one can fix a = a (ν) and solve recursively for a k+1, for each k : (3.1) a k+1 = 1 4 a k (k + 1)(k + ν + 1). We now establish an explicit solution to (3.1), making use of the Euler gamma function Γ(z). See Appendix B of these notes for a definition of Γ(z) and a derivation of some of its basic properties. For our present purposes, we mention that Γ(z) is well defined for z / {, 1,, 3,... }, and satisfies (3.13) Γ(z + 1) = zγ(z). Also, Γ(1) = 1, so, if k is a positive integer, (3.14) Γ(k + 1) = k! Moreover, (3.15) 1 is well defined for all z R, Γ(z) 1 = for z {, 1,, 3,... }. Γ(z) In light of these results, we see that a solution to (3.1) is given by (3.16) a k = ( 1 ) k 1 4 Γ(k + 1)Γ(k + ν + 1).

8 Furthermore, though (3.1) seems to break down for ν { 1,, 3,... }, (3.16) is well defined for all ν R, and all k Z +, and we have (3.17) 4(k + 1)(k + ν + 1)a k+1 + a k =, so (3.11) vanishes. Thus a solution to (3.9) is given by (3.18) J ν (t) = 1 ν ( 1) k t ) k. Γ(k + 1)Γ(k + ν + 1)( (We have multiplied the coefficients a k in (3.16) by ν, for consistency with standard conventions.) This is valid for all ν, t R. Convergence of (3.18) follows from the result 1 Γ(k + 1)Γ(k + ν + 1) 1 k!, for k, k + ν sufficiently large, which follows from the fact that Γ(z) + as z +, itself a direct consequence of the definition (B.1). Returning to (3.3), we have the Bessel function (3.19) J ν (t) = t ν J ν (t) = ( 1) k t ) k+ν, Γ(k + 1)Γ(k + ν + 1)( solving (1.1), for ν R, t (, ). Note that (1.1) for ν coincides with (1.1) for ν. Thus, both J ν (t) and J ν (t) solve (1.1). Now the leading term in (3.19) is 1 t ) ν (3.) J ν (t) = +. Γ(ν + 1)( As long as ν is not a negative integer, this leading term is nonvanishing, so (3.1) ν / Z = J ν (t) and J ν (t) are linearly independent solutions of (1.1). However, if ν is a negative integer, 1/Γ(ν +1) =, so (3.) does not establish such linear independence. Also, clearly J ν = J ν when ν =. Moreover, when ν = n is a negative integer, we have (3.) J n (t) = = k=n l= ( 1) k t ) k n Γ(k + 1)Γ(k n + 1)( ( 1) l+n Γ(l + 1)Γ(l + n + 1)( t = ( 1) n J n (t). ) l+n Thus we have not constructed two linearly independent solutions to (1.1) when ν is an integer. We turn to this task in 4.

9 4. The Bessel functions Y ν We have seen that J ν (t) and J ν (t) are linearly independent solutions to Bessel s equation (4.1) u (t) + 1 t u (t) + (1 ν t ) u(t) =, as long as ν is not an integer, while if n Z, then J n (t) = ( 1) n J n (t). We want to construct a basis of solutions to (4.1), uniformly good for all ν. This construction can be motivated by a calculation of the Wronskian. Generally, for a pair of solutions u 1 and u to a second order ODE (4.) a(t)u + b(t)u + c(t)u =, u 1 and u are linearly independent if and only if the Wronskian ( ) u1 u (4.3) W (t) = W (u 1, u )(t) = det u 1 u is nonvanishing. Now W (t) satisfies the first order ODE (4.4) W (t) = b(t) a(t) W (t). In the case of Bessel s equation (4.1), this becomes (4.5) W (t) = W (t), t so (4.6) W (t) = K t, for some K, independent of t, but perhaps depending on ν. To be specific, (4.7) W (J ν, J ν )(t) = K(ν). t We turn to the task of evaluating K(ν). Using the relations J ν (t) = t ν J ν (t) and J ν (t) = t ν J ν (t), we have (4.8) W (J ν, J ν )(t) = W (J ν, J ν )(t) ν t J ν(t)j ν (t).

1 Since J ν and J ν are smooth on R, we get (4.9) W (J ν, J ν )(t) = ν J ν()j ν () t + g(t), where g(t) is smooth on [, ). Comparison with (4.7) gives (4.1) W (J ν, J ν )(t) = νj ν()j ν (), t so, in (4.7), K(ν) = νj ν ()J ν (). Now (3.18) yields (4.11) J ν () = so 1 ν Γ(ν + 1), J ν() = ν (4.1) K(ν) = Γ(ν + 1)Γ(1 ν). 1 ν Γ(1 ν), By (3.13) (cf. (B.3)), Γ(ν + 1) = νγ(ν). Furthermore (cf. (B.8)), (4.13) Γ(ν)Γ(1 ν) = π sin πν. Hence (4.14) W (J ν, J ν )(t) = π sin πν. t In particular, this recovers our previous observation that J ν and J ν are linearly independent if and only if ν / Z. Taking a cue from (4.14), we set (4.15) Y ν (t) = J ν(t) cos πν J ν (t), sin πν when ν is not an integer. Then, for n Z, we define (4.16) Y n (t) = lim Y ν (t) = 1 [ ν n π ν J ν(t) ( 1) n ν(t)] ν J ν=n. Then we have (4.17) W (J ν, Y ν )(t) = πt, for all ν.

We now produce a series representation for Y n (t), when n Z +. Recall that, for ν R, ( t ) k+ν, ( 1) k (4.18) J ν (t) = α k (ν) αk (ν) = k!γ(k + ν + 1). Hence (4.19) ν J ν(t) = ( t k+ν ( α k(ν) ) + log t ) ( t ) k+ν, α k (ν) Note that the last sum is equal to J ν (t). Hence (4.) ν J ( t k+n ( ν(t) = α ν=n k(n) + log ) t ) J n (t). Similarly, (4.1) ν J ν(t) = ( t k ν ( α k( ν) ) log t ) J ν (t). Recalling from (3.13) that α k ( n) = for k < n, we have (4.) ν J ( t k n ( ν(t) = α ν=n k( n) log ) t ) J n (t). Plugging (4.) and (4.) into (4.16) then yields the series for Y n (t), convergent for t (, ). In detail, if n is a positive integer, (4.3) and Y n (t) = π ( log t (4.4) Y (t) = π )J n (t) + ( 1)n π + 1 π ( t α k(n) ( log t ) J (t) + π ( t k n α k( n) ) ) k+n, ( t ) k. α k() The evaluation of α k (±n) is discussed in Appendix C. We see that (4.5) Y (t) log t, as t, π while, if n is a positive integer, ( (4.6) Y n (t) ( 1)n t ) n, α π ( n) as t. As seen in Appendix C, when n is a positive integer, (4.7) α ( n) = β ( n + 1) = ( 1) n 1 (n 1)!. 11

1 5. Relations between J ν and J ν±1 The Bessel functions have the following remarkable relations: (5.1) (5.) J ν+1 (t) = J ν(t) + ν t J ν(t), J ν 1 (t) = J ν(t) + ν t J ν(t). Putting these two identities together yields (5.3) ( d dt ν 1 )( d t dt + ν ) J ν (t) = J ν (t), t which is equivalent to the Bessel equation (1.1). The identities (5.1) and (5.) are equivalent, respectively, to (5.4) J ν+1 (t) = 1 t J ν(t), and (5.5) J ν 1 (t) = tj ν(t) + νj ν (t). These follow readily from the power series (5.6) J ν (t) = 1 ν and its corollary (5.7) J ν(t) = 1 ν k=1 ( 1) k t ) k. k!γ(k + ν + 1)( ( 1) k k t ) k 1. k!γ(k + ν + 1)( The reader is invited to verify (5.4) (5.5), and hence (5.1) (5.). We mention that to establish (5.4) one uses Γ(k + 1 + 1) = (k + 1)Γ(k + 1), and to establish (5.5) one uses Γ(k + ν + 1) = (k + ν)γ(k + ν). We also have analogues of (5.1) (5.) for the functions Y ν : (5.8) Y ν+1 (t) = Y ν(t) + ν t Y ν(t), Y ν 1 (t) = Y ν(t) + ν t Y ν(t). We leave it to the reader to establish this, first for ν / Z, via (4.15), and then for ν Z, via a limiting argument.

13 6. The functions J n+1/ Looking at (3.9), we see that J 1/ (t) solves (6.1) u (t) + u(t) =, whose solutions are C 1 cos t + C sin t. Meanwhile, the power series (6.) J ν (t) = 1 ν yields (6.3) J 1/ (t) = ( 1) k t ) k k!γ(k + ν + 1)( ( 1) k t ) k, k!γ(k + 1/)( which has only even powers of t, so J 1/ (t) = C 1 cos t. Furthermore, (6.3) implies (6.4) J 1/ () = the latter identity by (B.9). Hence (6.5) J 1/ (t) = and consequently (6.6) J 1/ (t) = Applying (5.1) with ν = 1/ yields (6.7) J 1/ (t) = Γ(1/) = π, π πt πt cos t, cos t. sin t. Then repeated applications of (5.1) (5.) give, for n Z +, (6.8) J n+1/ (t) = ( 1) k{ n ( d dt j 1/ )} sin t, t πt j=1

14 and { n ( d (6.9) J n 1/ (t) = dt j 1/ )} cos t. t πt j=1 Returning to the identity (6.5), we recall that (6.1) cos t = Comparison with (6.3) gives (6.11) ( 1) k (k)! tk. π(k)! = k!γ( k k + 1 ), which can also be deduced from the identity Γ(1/) = π and, via Γ(z+1) = zγ(z), (6.1) ( Γ k + 1 ) ( = k 1 ) ( Γ k 1 ) = = k 1 k 3 1 ( 1 ) Γ.

15 7. Integral formula for J ν Our goal in this section is to establish the following integral formula: (7.1) J ν (t) = (t/) ν Γ(1/)Γ(ν + 1/) 1 1 (1 s ) ν 1/ e ist ds, for ν > 1/. Motivation for this formula can be found in Chapter 3, 6, of [T3]. To verify (7.1), we replace e ist by its power series, integrate term by term, and use an identity from Appendix B. To begin, the integral on the right side of (7.1) is equal to 1 1 (7.) (ist) k (1 s ) ν 1/ ds. (k)! 1 The identity (B.14) implies (7.3) 1 1 so the right side of (7.1) equals (7.4) s k (1 s ) ν 1/ ds = (t/) ν Γ(1/)Γ(ν + 1/) Γ(k + 1/)Γ(ν + 1/), Γ(k + ν + 1) 1 Γ(k + 1/)Γ(ν + 1/) (it)k. (k)! Γ(k + ν + 1) As seen in (6.11), we have ( 1 ( (7.5) Γ (k)! = ) k k! Γ k + 1 ), so the right side of (7.1) is equal to (7.6) ( t ) ν ( 1) k t ) k, k!γ(k + ν + 1)( which agrees with our formula (3.19) for J ν (t). This proves (7.1). Note that the integral in (7.1) is easy to evaluate when ν = 1/. We have (7.7) recovering the formula (6.7). J 1/ (t) = 1 ( t ) 1/ 1 e ist ds π 1 = 1 ( t ) 1/ sin t π t = sin t, πt

16 8. The Hankel functions H (k) ν The Hankel functions H ν (1) and H ν () are defined by (8.1) H ν (1) (t) = J ν (t) + iy ν (t), H ν () (t) = J ν (t) iy ν (t). The Wronskian relation (4.17) readily yields (8.) W (H ν (1), H ν () )(t) = 4i πt. Series expansions for H (1) ν (t) and H () ν (t) follow from (3.19) and (4.15) for ν / Z, and (4.3) (4.4) for ν = n Z +. Using (4.15), we have (8.3) H (1) ν H () ν (t) = J ν(t) e πiν J ν (t), i sin πν (t) = eπiν J ν (t) J ν (t), i sin πν evaluated as usual by l Hospital s rule for ν Z. These formulas imply (8.4) H (1) ν (t) = e πiν H ν (1) (t), H () ν (t) = e πiν H ν () (t). We mention the following useful integral formula: ( 1/ e (8.5) H ν (1) i(t πν/ π/4) (t) = e πt) s s ν 1/( 1 s ) ν 1/ ds, Γ(ν + 1/) it valid for t >, ν > 1/. There is a similar formula for H () ν (t). We will not prove (8.5) here. See [W], p. 168, for a derivation. See also [L], p. 139.

17 9. Large t behavior The integral formula (7.1) for J ν (t), valid for ν > 1/, t >, represents J ν (t) as (t/) ν times the Fourier transform of a function supported on the line segment I = [ 1, 1], and smooth on the interior of this segment. Classical methods of Fourier analysis (cf. [T3], Chapter 3) read off the asymptotic behavior of such a Fourier integral from the nature of the singularities of χ I (s)(1 s ) ν 1/ at the endpoints s = ±1. The following results: ( 1/ ( (9.1) J ν (t) = cos t πt) νπ π ) + O(t 3/ ), t +, 4 given ν > 1/. We also have ( 1/ ( (9.) J ν(t) = sin t πt) νπ π ) + O(t 3/ ), t +. 4 In fact, (9.) follows from the identity (5.1), i.e., (9.3) J ν(t) = J ν+1 (t) + ν t J ν(t), and (9.1), applied both to J ν (t) and to J ν+1 (t), and noting that ( (9.4) cos x π ) = sin x. The integral formula (8.5) for H ν (1) (t), valid for ν > 1/, t >, can be used to obtain an analogous result on the asymptotic behavior of H ν (1) (t). In fact, one can show that, if ν > 1/, (9.5) e s s ν 1/( 1 s ) ν 1/ ( ds = Γ ν + 1 ) + O(t 1 ), it as t +. It follows that (9.6) H (1) ν (t) = ( πt ) 1/e i(t νπ/ π/4) + O(t 3/ ), t +, given ν > 1/. Since H ν (1) (t) = J ν (t) + iy ν (t), comparison with (9.1) gives ( 1/ ( (9.7) Y ν (t) = sin t πt) νπ π ) + O(t 3/ ), t +, 4 for ν > 1/. An argument parallel to that yielding (9.), with (9.3) replaced by (5.8), yields ( 1/ ( (9.8) Y ν(t) = cos t πt) νπ π ) + O(t 3/ ), t +, 4 for ν > 1/. Note that a comparison with (9.1) (9.) gives (9.9) Y ν (t) = J ν(t) + O(t 3/ ), Y ν(t) = J ν (t) + O(t 3/ ), as t +, for ν > 1/.

18 1. Zeros of J ν and Y ν The functions J ν (t) and Y ν (t) each have an infinite number of zeros on the half line (, ). We show how this follows readily from the asymptotic formulas (9.1) and (9.7), and consider how to approximate these zeros, at least for large t. More generally, any real-valued solution u ν to Bessel s equation (1.1) is a linear combination u ν = aj ν + by ν, with a, b R, and hence there exists ψ [, π] and A R such that ( 1/ (1.1) u ν (t) = A cos (t πt) νπ π ) 4 ψ + O(t 3/ ), as t +. We assume u ν is not identically zero, so A. We can restrict attention to ν [, ), and note that (1.1) applies also to u ν = J ν and u ν = Y ν. We also have ( 1/ (1.) u ν(t) = A sin (t πt) νπ π ) 4 ψ + O(t 3/ ). Now set (1.3) α kν = kπ + νπ + π 4 + ψ, ( β kν = k + 1 ) π + νπ + π 4 + ψ = α kν + π. We have, with t = α kν, (1.4) ( ) 1/ u ν (α kν ) = ( 1) k A + O(t 3/ ), πt u ν (β kν ) = O(t 3/ ), and (still with t = α kν ) (1.5) u ν(α kν ) = O(t 3/ ), ( ) 1/ u ν(β kν ) = ( 1) k+1 A + O(t 3/ ). πt The first part of (1.4) guarantees that, for each ν [, ), there exists T ν < such that, whenever α kν T ν, u ν (α kν ) and u(α (k+1)ν ) = u ν (α kν +π) have opposite signs, hence there exists (1.6) ϑ kν (α kν, α kν + π) such that u ν (ϑ kν ) =.

The second part of (1.4) suggests that ϑ kν is close to β kν = α kν + π/, and indeed the second part of (1.5) yields (1.7) ϑ kν = β kν + O(k 1 ). 19 Remark. The Hankel functions have no zeros on the half line (, ). In fact, since J ν (t) and Y ν (t) are real valued for t (, ), (1.8) t (, ), H (1) ν (t ) = = J ν (t ) = Y ν (t ) =, contradicting the fact that {J ν, Y ν } is a basis of solutions to (1.1).

A. Where Bessel functions come from (A.1) Bessel functions arise in the natural generalization of the equation d u dx + k u =, with solutions sin kx and cos kx, to partial differential equations (A.) u + k u =, where is the Laplace operator, acting on a function u on a domain Ω R n by (A.3) u = u x 1 + + u x. n We can eliminate k from (A.) by scaling. Set u(x) = v(kx). Then equation (A.) becomes (A.4) ( + 1)v =. We specialize to the case n = and write (A.5) u = u x + u y. For a number of special domains Ω R, such as circular domains, annular domains, angular sectors, and pie-shaped domains, it is convenient to switch to polar coordinates (r, θ), related to (x, y)-coordinates by (A.6) x = r cos θ, y = r sin θ. In such coordinates, ( (A.7) v = r + 1 r r + 1 ) r θ v. A special class of solutions to (A.4) has the form (A.8) v = w(r)e iνθ. By (A.7), for such v, (A.9) ( + 1)v = [ d w dr + 1 dw ( ) ] r dr + 1 ν r w e iνθ,

1 so (A.4) holds if and only if (A.1) d w dr + 1 dw ( ) r dr + 1 ν r w =. This is Bessel s equation (1.1) (with different variables). Note that if v solves (A.4) on Ω R and if Ω is a circular domain or an annular domain, centered at the origin, then ν must be an integer. However, if Ω is an angular sector or a pie-shaped domain, with vertex at the origin, ν need not be an integer. In n dimensions, the Laplace operator (A.3) can be written ( (A.11) v = r + n 1 r r + 1 ) r S v, where S is a second-order differential operator acting on functions on the unit sphere S n 1 R n, called the Laplace-Beltrami operator. Generalizing (A.8), one looks for solutions to (A.4) of the form (A.1) v(x) = w(r)ψ(ω), where x = rω, r (, ), ω S n 1. Parallel to (A.9), for such v, (A.13) ( + 1)v = provided [ d w dr + n 1 dw ( ) ] r dr + 1 ν r w ψ(ω), (A.14) S ψ = ν ψ. The equation (A.15) d w dr + n 1 dw ( ) r dr + 1 ν r w = is a variant of Bessel s equation. If we set (A.16) ϕ(r) = r n/ 1 w(r), then (A.15) is converted into the Bessel equation (A.17) d ϕ dr + 1 dϕ ( ) r dr + 1 µ r ϕ =, µ = ν + ( n ). The study of solutions to (A.14) gives rise to the study of spherical harmonics, and from there to other special functions, such as Legendre functions.

The search for solutions of the form (A.1) is a key example of the method of separation of variables for partial differential equations. It arises in numerous other contexts. Here are a couple of other examples: (A.18) ( x + k )u =, and (A.19) ( + K x + k) u =. The first describes the n-dimensional quantum harmonic oscillator. The second (for n = 3) describes the quantum mechanical model of a hydrogen atom, according to Schrödinger. Study of these equations leads to other special functions defined by differential equations, such as Hermite functions and Whittaker functions. Much further material on these topics can be found in books on partial differential equations, such as [T3] (particularly Chapters 3 and 8).

3 B. The Euler gamma function Γ(z) We define the gamma function by (B.1) Γ(z) = e t t z 1 dt, for z >. Actually, we can take z to be complex, and (B.1) works for Re z >. In such a case, (B.1) is absolutely convergent. There is the crucial identity (B.) Γ(z + 1) = = = = zγ(z), e t t z dt d dt (e t ) t z dt e t d dt tz dt for Re z >, where we use integration by parts plus the identity (d/dt)t z = zt z 1. In other words, (B.3) Γ(z + 1) = zγ(z), for Re z >. The definition (B.1) readily gives Γ(1) = 1, so, for each positive integer k, the relation (B.3) gives inductively (B.4) Γ(k) = (k 1)! While Γ(z) is defined in (B.1) for Re z >, note that the left side of (B.3) is well defined for Re z > 1, so this identity extends Γ(z) to {z C : Re z > 1}, with a pole at z =. Iterating this argument, we extend Γ(z) to {z C : z, 1,,... }. Similarly, (B.5) β(z) = 1 Γ(z), initially defined for Re z >, satisfies (B.6) β(z) = zβ(z + 1), and then β extends to be well defined for all z C, with (B.7) β(z) = for z {, 1,, 3,... }.

4 For use in 4 (cf. (4.13)) we record the following identity: (B.8) Γ(z)Γ(1 z) = π sin πz. A demonstration of this would take us far afield, so we refer to 18 of [T] for a proof, and for further material on the gamma function. Note that setting z = 1/ in (B.8) yields (B.9) ( 1 Γ = ) π. A direct proof of (B.9) can be given as follows. Taking t = s in (B.1) gives (B.1) Γ(z) = e s s z 1 ds, hence (B.11) Thus (B.1) ( 1 Γ = e ) s ds = e s ds. ( 1 ( Γ = ) )( ) e x dx e y dy = = = π π R e r r dr dθ = π e r r dr e s ds = π. e (x +y ) dx dy The following identity will be useful in 7. Namely, the beta function, defined for x, y > by (B.13) B(x, y) = 1 (with u = s/(1 s)), satisfies s x 1 (1 s) y 1 ds = (1 + u) x y u x 1 du (B.14) B(x, y) = Γ(x)Γ(y) Γ(x + y). To prove this, note that since (B.15) Γ(z)p z = e pt t z 1 dt,

5 we have (B.16) (1 + u) x y = so 1 e (1+u)t t x+y 1 dt, Γ(x + y) (B.17) 1 B(x, y) = e t t x+y 1 e ut u x 1 du dt Γ(x + y) = Γ(x) e t t y 1 dt Γ(x + y) = Γ(x)Γ(y) Γ(x + y), as asserted. For closer contact with (7.3), note that setting s = t in (B.13) gives 1 (B.18) B(x, y) = t x 1 (1 t ) y 1 dt, so, if k Z + and ν > 1/, (B.19) ( B k + 1, ν + 1 ) 1 = t k (1 t ) ν 1/ dt. 1

6 C. Derivative of 1/Γ(z) at integer points α k The formula for Y n (t), given by (4.16), (4.), and (4.), involves α k (n) and ( n), where (C.1) α k (ν) = Thus ( 1) k k!γ(k + ν + 1). (C.) α k(ν) = ( 1)k β (k + ν + 1), k! where, as in (B.5), (C.3) β(z) = 1 Γ(z). Thus we want to evaluate β (k ± n + 1), or equivalently, evaluate β (l) for all l Z. To start, we have from (B.1) that, for z > 1, (C.4) Γ (z) = In particular, (C.5) Γ (1) = where γ is Euler s constant, e t t z 1 (log t) dt. e t log t dt = γ, (C.6) γ.577156649. See 18 and Appendix J of [T] for more on this constant. We deduce that (C.7) β (1) = γ. To evaluate β (l) for other l Z, we can use (C.8) β(z) = zβ(z + 1), which implies (C.9) β (z) = β(z + 1) + zβ (z + 1).

7 Thus (C.1) β () = β(1) = 1. For l = m, a negative integer, we have β( m + 1) =, hence (C.11) β ( m) = mβ ( m + 1), hence β ( 1) = β () = 1, β ( ) = β ( 1) =, and, inductively, (C.1) β ( m) = ( 1) m m!, when m is a positive integer. To evaluate β(l) for an integer l, we can turn (C.9) around: (C.13) β (z + 1) = β (z) β(z + 1). z Hence (C.14) β () = β (1) β() = γ 1, and generally (C.15) β (l + 1) = β (l) l 1 l l!.

8 References [L] N. Lebedev, Special Functions and Their Applications, Dover, New York, 197. [O] F. Olver, Asymptotics and Special Functions, Academic Press, New York, 1974. [T1] M. Taylor, Introduction to Differential Equations, American Mathematical Society, Providence, RI, 11. [T] M. Taylor, Introduction to Complex Analysis, Lecture Notes, available at http://www.unc.edu/math/faculty/met/complex.html [T3] M. Taylor, Partial Differential Equations, Vols. 1 3, Springer-Verlag, New York, 1996 (nd ed. 11). [W] G. Watson, A Treatise on the Theory of Bessel Functions, Cambridge Univ. Press, Cambridge UK, 1944 (Library Ed. 1996).