Bessel Functions Michael Taylor. Lecture Notes for Math 524
|
|
- Christina Kelley
- 6 years ago
- Views:
Transcription
1 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
2 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 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 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 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 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 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 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 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).
10 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 ν.
11 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
12 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 ) = (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 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 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 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 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 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 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ν ) =.
19 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).
20 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 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νθ,
21 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.
22 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).
23 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,... }.
24 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,
25 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
26 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) γ 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).
27 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!.
28 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, [T1] M. Taylor, Introduction to Differential Equations, American Mathematical Society, Providence, RI, 11. [T] M. Taylor, Introduction to Complex Analysis, Lecture Notes, available at [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).
S. Ghorai 1. Lecture XV Bessel s equation, Bessel s function. e t t p 1 dt, p > 0. (1)
S Ghorai 1 1 Gamma function Gamma function is defined by Lecture XV Bessel s equation, Bessel s function Γp) = e t t p 1 dt, p > 1) The integral in 1) is convergent that can be proved easily Some special
More informationSurvival Guide to Bessel Functions
Survival Guide to Bessel Functions December, 13 1 The Problem (Original by Mike Herman; edits and additions by Paul A.) For cylindrical boundary conditions, Laplace s equation is: [ 1 s Φ ] + 1 Φ s s s
More informationENGI 9420 Lecture Notes 1 - ODEs Page 1.01
ENGI 940 Lecture Notes - ODEs Page.0. Ordinary Differential Equations An equation involving a function of one independent variable and the derivative(s) of that function is an ordinary differential equation
More informationThis ODE arises in many physical systems that we shall investigate. + ( + 1)u = 0. (λ + s)x λ + s + ( + 1) a λ. (s + 1)(s + 2) a 0
Legendre equation This ODE arises in many physical systems that we shall investigate We choose We then have Substitution gives ( x 2 ) d 2 u du 2x 2 dx dx + ( + )u u x s a λ x λ a du dx λ a λ (λ + s)x
More informationIntroduction to Spherical Harmonics
Introduction to Spherical Harmonics Lawrence Liu 3 June 4 Possibly useful information. Legendre polynomials. Rodrigues formula:. Generating function: d n P n x = x n n! dx n n. wx, t = xt t = P n xt n,
More informationSolutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions. ρ + (1/ρ) 2 V
Solutions to Laplace s Equation in Cylindrical Coordinates and Numerical solutions Lecture 8 1 Introduction Solutions to Laplace s equation can be obtained using separation of variables in Cartesian and
More informationStrauss PDEs 2e: Section Exercise 4 Page 1 of 6
Strauss PDEs 2e: Section 5.3 - Exercise 4 Page of 6 Exercise 4 Consider the problem u t = ku xx for < x < l, with the boundary conditions u(, t) = U, u x (l, t) =, and the initial condition u(x, ) =, where
More information7.3 Singular points and the method of Frobenius
284 CHAPTER 7. POWER SERIES METHODS 7.3 Singular points and the method of Frobenius Note: or.5 lectures, 8.4 and 8.5 in [EP], 5.4 5.7 in [BD] While behaviour of ODEs at singular points is more complicated,
More informationLecture 4b. Bessel functions. Introduction. Generalized factorial function. 4b.1. Using integration by parts it is easy to show that
4b. Lecture 4b Using integration by parts it is easy to show that Bessel functions Introduction In the previous lecture the separation of variables method led to Bessel's equation y' ' y ' 2 y= () 2 Here
More informationAppendix A Vector Analysis
Appendix A Vector Analysis A.1 Orthogonal Coordinate Systems A.1.1 Cartesian (Rectangular Coordinate System The unit vectors are denoted by x, ŷ, ẑ in the Cartesian system. By convention, ( x, ŷ, ẑ triplet
More informationNotes on Special Functions
Spring 25 1 Notes on Special Functions Francis J. Narcowich Department of Mathematics Texas A&M University College Station, TX 77843-3368 Introduction These notes are for our classes on special functions.
More informationPhysics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I
Physics 342 Lecture 23 Radial Separation Lecture 23 Physics 342 Quantum Mechanics I Friday, March 26th, 2010 We begin our spherical solutions with the simplest possible case zero potential. Aside from
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationPHYS 404 Lecture 1: Legendre Functions
PHYS 404 Lecture 1: Legendre Functions Dr. Vasileios Lempesis PHYS 404 - LECTURE 1 DR. V. LEMPESIS 1 Legendre Functions physical justification Legendre functions or Legendre polynomials are the solutions
More informationLecture 31. Basic Theory of First Order Linear Systems
Math 245 - Mathematics of Physics and Engineering I Lecture 31. Basic Theory of First Order Linear Systems April 4, 2012 Konstantin Zuev (USC) Math 245, Lecture 31 April 4, 2012 1 / 10 Agenda Existence
More informationIn this chapter we study several functions that are useful in calculus and other areas of mathematics.
Calculus 5 7 Special functions In this chapter we study several functions that are useful in calculus and other areas of mathematics. 7. Hyperbolic trigonometric functions The functions we study in this
More informationConformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder.
Lent 29 COMPLEX METHODS G. Taylor A star means optional and not necessarily harder. Conformal maps. (i) Let f(z) = az + b, with ad bc. Where in C is f conformal? cz + d (ii) Let f(z) = z +. What are the
More informationODE Homework Solutions of Linear Homogeneous Equations; the Wronskian
ODE Homework 3 3.. Solutions of Linear Homogeneous Equations; the Wronskian 1. Verify that the functions y 1 (t = e t and y (t = te t are solutions of the differential equation y y + y = 0 Do they constitute
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More informationOn an Eigenvalue Problem Involving Legendre Functions by Nicholas J. Rose North Carolina State University
On an Eigenvalue Problem Involving Legendre Functions by Nicholas J. Rose North Carolina State University njrose@math.ncsu.edu 1. INTRODUCTION. The classical eigenvalue problem for the Legendre Polynomials
More informationBessel s Equation. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Bessel s Equation MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Background Bessel s equation of order ν has the form where ν is a constant. x 2 y + xy
More informationODE Homework Series Solutions Near an Ordinary Point, Part I 1. Seek power series solution of the equation. n(n 1)a n x n 2 = n=0
ODE Homework 6 5.2. Series Solutions Near an Ordinary Point, Part I 1. Seek power series solution of the equation y + k 2 x 2 y = 0, k a constant about the the point x 0 = 0. Find the recurrence relation;
More information4 Power Series Solutions: Frobenius Method
4 Power Series Solutions: Frobenius Method Now the ODE adventure takes us to series solutions for ODEs, a technique A & W, that is often viable, valuable and informative. These can be readily applied Sec.
More information23 Elements of analytic ODE theory. Bessel s functions
23 Elements of analytic ODE theory. Bessel s functions Recall I am changing the variables) that we need to solve the so-called Bessel s equation 23. Elements of analytic ODE theory Let x 2 u + xu + x 2
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationSeries Solutions of ODEs. Special Functions
C05.tex 6/4/0 3: 5 Page 65 Chap. 5 Series Solutions of ODEs. Special Functions We continue our studies of ODEs with Legendre s, Bessel s, and the hypergeometric equations. These ODEs have variable coefficients
More informationAsymptotic of Enumerative Invariants in CP 2
Peking Mathematical Journal https://doi.org/.7/s4543-8-4-4 ORIGINAL ARTICLE Asymptotic of Enumerative Invariants in CP Gang Tian Dongyi Wei Received: 8 March 8 / Revised: 3 July 8 / Accepted: 5 July 8
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 21 Quantum Mechanics in Three Dimensions Lecture 21 Physics 342 Quantum Mechanics I Monday, March 22nd, 21 We are used to the temporal separation that gives, for example, the timeindependent
More informationRiemann integral and volume are generalized to unbounded functions and sets. is an admissible set, and its volume is a Riemann integral, 1l E,
Tel Aviv University, 26 Analysis-III 9 9 Improper integral 9a Introduction....................... 9 9b Positive integrands................... 9c Special functions gamma and beta......... 4 9d Change of
More informationDifferential Equations
Electricity and Magnetism I (P331) M. R. Shepherd October 14, 2008 Differential Equations The purpose of this note is to provide some supplementary background on differential equations. The problems discussed
More informationLaplace s Equation in Cylindrical Coordinates and Bessel s Equation (I)
Laplace s Equation in Cylindrical Coordinates and Bessel s Equation I) 1 Solution by separation of variables Laplace s equation is a key equation in Mathematical Physics. Several phenomena involving scalar
More informationORDINARY DIFFERENTIAL EQUATIONS
ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 4884 NOVEMBER 9, 7 Summary This is an introduction to ordinary differential equations We
More informationLaplace s Method for Ordinary Differential Equations
Physics 2400 Spring 2016 Laplace s Method for Ordinary Differential Equations Lecture notes by M. G. Rozman Last modified: March 31, 2016 We can [1, pp. 124-8], [2, h. VIII], [3, Apps. a-b], [4, h. 18],
More informationAnalogues for Bessel Functions of the Christoffel-Darboux Identity
Analogues for Bessel Functions of the Christoffel-Darboux Identity Mark Tygert Research Report YALEU/DCS/RR-1351 March 30, 2006 Abstract We derive analogues for Bessel functions of what is known as the
More informationASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS. Tian Ma. Shouhong Wang
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 11, Number 1, July 004 pp. 189 04 ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS Tian Ma Department of
More informationAnswers to Problem Set Number MIT (Fall 2005).
Answers to Problem Set Number 5. 18.305 MIT Fall 2005). D. Margetis and R. Rosales MIT, Math. Dept., Cambridge, MA 02139). November 23, 2005. Course TA: Nikos Savva, MIT, Dept. of Mathematics, Cambridge,
More informationMATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012
MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to
More informationRecurrence Relations and Fast Algorithms
Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations
More informationQualification Exam: Mathematical Methods
Qualification Exam: Mathematical Methods Name:, QEID#41534189: August, 218 Qualification Exam QEID#41534189 2 1 Mathematical Methods I Problem 1. ID:MM-1-2 Solve the differential equation dy + y = sin
More informationMath 300 Winter 2011 Advanced Boundary Value Problems I. Bessel s Equation and Bessel Functions
Math 3 Winter 2 Avance Bounary Value Problems I Bessel s Equation an Bessel Functions Department of Mathematical an Statistical Sciences University of Alberta Bessel s Equation an Bessel Functions We use
More informationNotes on the CEV model
Notes on the CEV model Andrew Lesniewski Ellington Management Group 53 Forest Avenue Old Greenwich, CT 687 First draft of March 2, 24 This draft of November 6, 29 Statement of the problem We construct
More informationMath 259: Introduction to Analytic Number Theory More about the Gamma function
Math 59: Introduction to Analytic Number Theory More about the Gamma function We collect some more facts about Γs as a function of a complex variable that will figure in our treatment of ζs and Ls, χ.
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationFIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS II: Homogeneous Linear Systems with Constant Coefficients
FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS II: Homogeneous Linear Systems with Constant Coefficients David Levermore Department of Mathematics University of Maryland 28 January 2012 Because
More informationSeries Solutions. 8.1 Taylor Polynomials
8 Series Solutions 8.1 Taylor Polynomials Polynomial functions, as we have seen, are well behaved. They are continuous everywhere, and have continuous derivatives of all orders everywhere. It also turns
More informationMATH5685 Assignment 3
MATH5685 Assignment 3 Due: Wednesday 3 October 1. The open unit disk is denoted D. Q1. Suppose that a n for all n. Show that (1 + a n) converges if and only if a n converges. [Hint: prove that ( N (1 +
More informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georgia Tech PHYS 612 Mathematical Methods of Physics I Instructor: Predrag Cvitanović Fall semester 2012 Homework Set #5 due October 2, 2012 == show all your work for maximum credit, == put labels, title,
More informationTopic 13 Notes Jeremy Orloff
Topic 13 Notes Jeremy Orloff 13 Analytic continuation and the Gamma function 13.1 Introduction In this topic we will look at the Gamma function. This is an important and fascinating function that generalizes
More informationx exp ( x 2 y 2) f (x) dx. (1.2) The following Laplace-type transforms which are the L 2n transform and the L 4n transform,
Journal of Inequalities and Special Functions ISSN: 17-433, URL: www.ilirias.com/jiasf Volume 8 Issue 517, Pages 75-85. IDENTITIES ON THE E,1 INTEGRAL TRANSFORM NEŞE DERNEK*, GÜLSHAN MANSIMLI, EZGI ERDOĞAN
More informationThe Laplacian and Mean and Extreme Values
Mathematical Assoc. of America American Mathematical Monthly 2: September 29, 205 2:9 p.m. LaplaceMeanValueFinal.tex page The Laplacian and Mean and Extreme Values Jeffrey S. Ovall Abstract. The Laplace
More informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationζ (s) = s 1 s {u} [u] ζ (s) = s 0 u 1+sdu, {u} Note how the integral runs from 0 and not 1.
Problem Sheet 3. From Theorem 3. we have ζ (s) = + s s {u} u+sdu, (45) valid for Res > 0. i) Deduce that for Res >. [u] ζ (s) = s u +sdu ote the integral contains [u] in place of {u}. ii) Deduce that for
More informationUNIFORM BOUNDS FOR BESSEL FUNCTIONS
Journal of Applied Analysis Vol. 1, No. 1 (006), pp. 83 91 UNIFORM BOUNDS FOR BESSEL FUNCTIONS I. KRASIKOV Received October 8, 001 and, in revised form, July 6, 004 Abstract. For ν > 1/ and x real we shall
More informationHIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS. David Levermore Department of Mathematics University of Maryland.
HIGHER-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATIONS David Levermore Department of Mathematics University of Maryland 28 March 2008 The following is a review of some of the material that we covered on higher-order
More information221A Lecture Notes Steepest Descent Method
Gamma Function A Lecture Notes Steepest Descent Method The best way to introduce the steepest descent method is to see an example. The Stirling s formula for the behavior of the factorial n! for large
More informationThe integrating factor method (Sect. 1.1)
The integrating factor method (Sect. 1.1) Overview of differential equations. Linear Ordinary Differential Equations. The integrating factor method. Constant coefficients. The Initial Value Problem. Overview
More informationBoundary Value Problems in Cylindrical Coordinates
Boundary Value Problems in Cylindrical Coordinates 29 Outline Differential Operators in Various Coordinate Systems Laplace Equation in Cylindrical Coordinates Systems Bessel Functions Wave Equation the
More informationRelevant sections from AMATH 351 Course Notes (Wainwright): 1.3 Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls): 1.1.
Lecture 8 Qualitative Behaviour of Solutions to ODEs Relevant sections from AMATH 351 Course Notes (Wainwright): 1.3 Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls): 1.1.1 The last few
More informationRelevant sections from AMATH 351 Course Notes (Wainwright): Relevant sections from AMATH 351 Course Notes (Poulin and Ingalls):
Lecture 5 Series solutions to DEs Relevant sections from AMATH 35 Course Notes (Wainwright):.4. Relevant sections from AMATH 35 Course Notes (Poulin and Ingalls): 2.-2.3 As mentioned earlier in this course,
More informationMATH 311 Topics in Applied Mathematics Lecture 25: Bessel functions (continued).
MATH 311 Topics in Applied Mathematics Lecture 25: Bessel functions (continued). Bessel s differential equation of order m 0: z 2 d2 f dz 2 + z df dz + (z2 m 2 )f = 0 The equation is considered on the
More informationOn the Weyl symbol of the resolvent of the harmonic oscillator
On the Weyl symbol of the resolvent of the harmonic oscillator Jan Dereziński, Maciej Karczmarczyk Department of Mathematical Methods in Physics, Faculty of Physics University of Warsaw, Pasteura 5, -93,
More informationWe denote the derivative at x by DF (x) = L. With respect to the standard bases of R n and R m, DF (x) is simply the matrix of partial derivatives,
The derivative Let O be an open subset of R n, and F : O R m a continuous function We say F is differentiable at a point x O, with derivative L, if L : R n R m is a linear transformation such that, for
More informationMath 1B Final Exam, Solution. Prof. Mina Aganagic Lecture 2, Spring (6 points) Use substitution and integration by parts to find:
Math B Final Eam, Solution Prof. Mina Aganagic Lecture 2, Spring 20 The eam is closed book, apart from a sheet of notes 8. Calculators are not allowed. It is your responsibility to write your answers clearly..
More informationTaylor and Laurent Series
Chapter 4 Taylor and Laurent Series 4.. Taylor Series 4... Taylor Series for Holomorphic Functions. In Real Analysis, the Taylor series of a given function f : R R is given by: f (x + f (x (x x + f (x
More informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationChapter 4. Series Solutions. 4.1 Introduction to Power Series
Series Solutions Chapter 4 In most sciences one generation tears down what another has built and what one has established another undoes. In mathematics alone each generation adds a new story to the old
More informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
More informationORDINARY DIFFERENTIAL EQUATIONS
ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 48824. JANUARY 3, 25 Summary. This is an introduction to ordinary differential equations.
More informationMATH 32A: MIDTERM 2 REVIEW. sin 2 u du z(t) = sin 2 t + cos 2 2
MATH 3A: MIDTERM REVIEW JOE HUGHES 1. Curvature 1. Consider the curve r(t) = x(t), y(t), z(t), where x(t) = t Find the curvature κ(t). 0 cos(u) sin(u) du y(t) = Solution: The formula for curvature is t
More informationLecture 4.6: Some special orthogonal functions
Lecture 4.6: Some special orthogonal functions Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4340, Advanced Engineering Mathematics
More informationPractice Problems For Test 3
Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)
More informationCOMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE
BIT 6-85//41-11 $16., Vol. 4, No. 1, pp. 11 118 c Swets & Zeitlinger COMPUTATION OF BESSEL AND AIRY FUNCTIONS AND OF RELATED GAUSSIAN QUADRATURE FORMULAE WALTER GAUTSCHI Department of Computer Sciences,
More informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationAsymptotics of Integrals of. Hermite Polynomials
Applied Mathematical Sciences, Vol. 4, 010, no. 61, 04-056 Asymptotics of Integrals of Hermite Polynomials R. B. Paris Division of Complex Systems University of Abertay Dundee Dundee DD1 1HG, UK R.Paris@abertay.ac.uk
More informationarxiv:math/ v2 [math.ap] 3 Oct 2006
THE TAYLOR SERIES OF THE GAUSSIAN KERNEL arxiv:math/0606035v2 [math.ap] 3 Oct 2006 L. ESCAURIAZA From some people one can learn more than mathematics Abstract. We describe a formula for the Taylor series
More informationA review: The Laplacian and the d Alembertian. j=1
Chapter One A review: The Laplacian and the d Alembertian 1.1 THE LAPLACIAN One of the main goals of this course is to understand well the solution of wave equation both in Euclidean space and on manifolds
More informationMATH 425, HOMEWORK 3 SOLUTIONS
MATH 425, HOMEWORK 3 SOLUTIONS Exercise. (The differentiation property of the heat equation In this exercise, we will use the fact that the derivative of a solution to the heat equation again solves the
More informationExercises for Multivariable Differential Calculus XM521
This document lists all the exercises for XM521. The Type I (True/False) exercises will be given, and should be answered, online immediately following each lecture. The Type III exercises are to be done
More information5.4 Bessel s Equation. Bessel Functions
SEC 54 Bessel s Equation Bessel Functions J (x) 87 # with y dy>dt, etc, constant A, B, C, D, K, and t 5 HYPERGEOMETRIC ODE At B (t t )(t t ), t t, can be reduced to the hypergeometric equation with independent
More informationSeries solutions to a second order linear differential equation with regular singular points
Physics 6C Fall 0 Series solutions to a second order linear differential equation with regular singular points Consider the second-order linear differential equation, d y dx + p(x) dy x dx + q(x) y = 0,
More informationFinal Exam. Monday March 19, 3:30-5:30pm MAT 21D, Temple, Winter 2018
Name: Student ID#: Section: Final Exam Monday March 19, 3:30-5:30pm MAT 21D, Temple, Winter 2018 Show your work on every problem. orrect answers with no supporting work will not receive full credit. Be
More informationELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS EXERCISES I (HARMONIC FUNCTIONS)
ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS EXERCISES I (HARMONIC FUNCTIONS) MATANIA BEN-ARTZI. BOOKS [CH] R. Courant and D. Hilbert, Methods of Mathematical Physics, Vol. Interscience Publ. 962. II, [E] L.
More informationA Note on the Differential Equations with Distributional Coefficients
MATEMATIKA, 24, Jilid 2, Bil. 2, hlm. 115 124 c Jabatan Matematik, UTM. A Note on the Differential Equations with Distributional Coefficients Adem Kilicman Department of Mathematics, Institute for Mathematical
More information1 Solutions in cylindrical coordinates: Bessel functions
1 Solutions in cylindrical coordinates: Bessel functions 1.1 Bessel functions Bessel functions arise as solutions of potential problems in cylindrical coordinates. Laplace s equation in cylindrical coordinates
More informationMathematics of Physics and Engineering II: Homework answers You are encouraged to disagree with everything that follows. P R and i j k 2 1 1
Mathematics of Physics and Engineering II: Homework answers You are encouraged to disagree with everything that follows Homework () P Q = OQ OP =,,,, =,,, P R =,,, P S = a,, a () The vertex of the angle
More informationExistence Theory: Green s Functions
Chapter 5 Existence Theory: Green s Functions In this chapter we describe a method for constructing a Green s Function The method outlined is formal (not rigorous) When we find a solution to a PDE by constructing
More informationA FEW REMARKS ON THE SUPREMUM OF STABLE PROCESSES P. PATIE
A FEW REMARKS ON THE SUPREMUM OF STABLE PROCESSES P. PATIE Abstract. In [1], Bernyk et al. offer a power series and an integral representation for the density of S 1, the maximum up to time 1, of a regular
More informationWaves and the Schroedinger Equation
Waves and the Schroedinger Equation 5 april 010 1 The Wave Equation We have seen from previous discussions that the wave-particle duality of matter requires we describe entities through some wave-form
More informationNotes on uniform convergence
Notes on uniform convergence Erik Wahlén erik.wahlen@math.lu.se January 17, 2012 1 Numerical sequences We begin by recalling some properties of numerical sequences. By a numerical sequence we simply mean
More informationProblem 1 In each of the following problems find the general solution of the given differential
VI Problem 1 dt + 2dy 3y = 0; dt 9dy + 9y = 0. Problem 2 dt + dy 2y = 0, y(0) = 1, y (0) = 1; dt 2 y = 0, y( 2) = 1, y ( 2) = Problem 3 Find the solution of the initial value problem 2 d2 y dt 2 3dy dt
More information1+t 2 (l) y = 2xy 3 (m) x = 2tx + 1 (n) x = 2tx + t (o) y = 1 + y (p) y = ty (q) y =
DIFFERENTIAL EQUATIONS. Solved exercises.. Find the set of all solutions of the following first order differential equations: (a) x = t (b) y = xy (c) x = x (d) x = (e) x = t (f) x = x t (g) x = x log
More information1 Bessel, Neumann, and Hankel Functions: J ν (x), N ν (x), H (1)
Bessel, Neumann, and Hankel Functions: J ν x), N ν x), H ) ν x), H ) ν x) Bessel functions are solutions of the following differential equation: x y +xy +x ν )y =.) Any two of the following functions are
More information16. Local theory of regular singular points and applications
16. Local theory of regular singular points and applications 265 16. Local theory of regular singular points and applications In this section we consider linear systems defined by the germs of meromorphic
More informationSpecial Functions. SMS 2308: Mathematical Methods
Special s Special s SMS 2308: Mathematical Methods Department of Computational and Theoretical Sciences, Kulliyyah of Science, International Islamic University Malaysia. Sem 1 2014/2015 Factorial Bessel
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationChapter 3a Topics in differentiation. Problems in differentiation. Problems in differentiation. LC Abueg: mathematical economics
Chapter 3a Topics in differentiation Lectures in Mathematical Economics L Cagandahan Abueg De La Salle University School of Economics Problems in differentiation Problems in differentiation Problem 1.
More informationd Wave Equation. Rectangular membrane.
1 ecture1 1.1 2-d Wave Equation. Rectangular membrane. The first problem is for the wave equation on a rectangular domain. You can interpret this as a problem for determining the displacement of a flexible
More informationUNIVERSITY OF MANITOBA
Question Points Score INSTRUCTIONS TO STUDENTS: This is a 6 hour examination. No extra time will be given. No texts, notes, or other aids are permitted. There are no calculators, cellphones or electronic
More informationPhysics 200 Lecture 4. Integration. Lecture 4. Physics 200 Laboratory
Physics 2 Lecture 4 Integration Lecture 4 Physics 2 Laboratory Monday, February 21st, 211 Integration is the flip-side of differentiation in fact, it is often possible to write a differential equation
More information