Arxiv:1705.01190V1 [Math.CA]
Total Page:16
File Type:pdf, Size:1020Kb
Noname manuscript No. (will be inserted by the editor) Uniform asymptotic expansions for Laguerre polynomials and related confluent hypergeometric functions T. M. Dunster · A. Gil · J. Segura Received: date / Accepted: date Abstract Uniform asymptotic expansions involving exponential and Airy (α) functions are obtained for Laguerre polynomials Ln (x), as well as comple- mentary confluent hypergeometric functions. The expansions are valid for n large and α small or large, uniformly for unbounded real and complex values of (α) x. The new expansions extend the range of computability of Ln (x) compared to previous expansions, in particular with respect to higher terms and large values of α. Numerical evidence of their accuracy for real and complex values of x is provided. Keywords Asymptotic expansions Laguerre polynomials Confluent hypergeometric functions Turning point· theory, WKB methods· Numerical computation · · Mathematics Subject Classification (2000) MSC 34E05 33C45 33C15 34E20 33F05 · · · · The authors acknowledge support from Ministerio de Econom´ıa y Competitividad, project MTM2015-67142-P (MINECO/FEDER, UE). T. M. Dunster Department of Mathematics and Statistics San Diego State University. 5500 Campanile Drive San Diego, CA, USA. E-mail: [email protected] A. Gil Departamento de Matem´atica Aplicada y CC. de la Computaci´on. ETSI Caminos. Universidad de Cantabria. 39005-Santander, Spain. arXiv:1705.01190v1 [math.CA] 2 May 2017 E-mail: [email protected] J. Segura Departamento de Matem´aticas, Estad´ıstica y Computaci´on. Universidad de Cantabria. 39005-Santander, Spain. E-mail: [email protected] 2 T. M. Dunster et al. 1 Introduction In this paper we shall obtain computable asymptotic expansions for Laguerre (α) polynomials Ln (x), for the case n large and α small or large. These are defined by n n + α ( x)k L(α)(x)= − . (1.1) n n k k! Xk=0 − In terms of the confluent hypergeometric functions we have Γ (n + α + 1) ( 1)n L(α)(x)= M( n, α +1, x)= − U( n, α +1, x), (1.2) n n! − n! − where M denotes Olver’s scaled confluent hypergeometric function [9, Chap. 7, sec. 9]. We will use the techniques described in [3] for computing uniform asymp- totic expansions of turning point problems. These involve Airy function expan- sions for solutions of differential equations having a simple turning point, but in a form that differs from the classical Airy function expansions of [9, Chap. 11]. Specifically, the coefficients in these new expansions are significantly easier to compute than previously, and close to the turning point the method utilises Cauchy’s integral formula. We remark that the Cauchy integral method of [3] has potential applications to other forms of differential equations, and also the method was subsequently used in [4] to compute coefficients appearing in certain asymptotic expansions of integrals. In this paper we shall also obtain Liouville-Green (L-G) expansions for Laguerre polynomials and confluent hypergeometric functions in domains that do not contain turning points, and these involve the exponential function. The form of these expansions differs from the standard ones (see [9, Chap. 10]), insomuch the coefficients in the expansions appear in the argument of the exponential. The advantage of this form is that the coefficients are again generally easier to compute. Indeed, L-G expansions of this form were used in [3] to obtain the new form of the Airy expansions described above. A recent reference on the computation of Laguerre polynomials using asymp- totics is [6], where three different asymptotic approximations are used: two expansions in terms of Bessel functions (from [5] and [12, Sect. 10.3.4]) and the Airy expansion of Frenzen and Wong [5]. These expansions, and in par- ticular the Airy expansion, give accurate approximations for large n when α is not large. For large α some expansions are available, like for instance the expansions for Whittaker functions of [1] and [8], or some of the expansions in [11], but the coefficients of these expansions are very hard to compute. In this paper, we provide expansions that are computable and which can also be used for large α, and we provide numerical evidence of their accuracy for real and complex values of x. The Laguerre polynomials satisfy the following form of the confluent hy- pergeometric equation d2y dy x + (α +1 x) + ny =0. (1.3) dx2 − dx Uniform asymptotic expansions for Laguerre polynomials 3 (α+1)/2 1 By the transformationy ˜ = x exp 2 x y we can remove the first derivative in the usual manner, and we get − d2y˜ x2 2(2n + α + 1)x + α2 1 = − − y.˜ (1.4) dx2 4x2 1 On replacing x by uz, where u = n + 2 (a choice being explained after (2.9) below), this can then be re-written in the form d2w = u2f(a,z)+ g(z) w, (1.5) dz2 (α+1)/2 1 (α) with a solution w = z exp uz Ln (uz). On defining a by − 2 α = u a2 1 (a 0), (1.6) − ≥ we have in (1.5) (z z )(z z ) 1 f(z)= − 1 − 2 , g(z)= , (1.7) 4z2 −4z2 where z = (a 1)2, z = (a + 1)2. (1.8) 1 − 2 The latter are the turning points for large u, which is assumed here. These turning points coalesce with each other when a = 0 (α = u), and z1 coalesces with the pole at z = 0 when a = 1 (α = 0). − We shall consider the following cases separately. Case 1a: We obtain expansions for z lying in domains containing the turn- ing point z , which includes the interval 0 z z δ (here and throughout 1 ≤ ≤ 2 − δ > 0). The turning point z = z2 is excluded, and this is covered in Case 2 below. We assume 1+ δ a2 a < for fixed a (1, ). Thus z cannot ≤ ≤ 1 ∞ 1 ∈ ∞ 1 coalesce with the pole at z = 0 nor with the other turning point z = z2, and α is positive and large, satisfying uδ α u(a1 1). Case 1b: This is the same as Case≤ 1a,≤ except− now we consider α negative. 2 In particular, we assume 0 < a0 a 1 δ (for fixed a0 (0, 1)), and in this case u(1 a ) α uδ <≤ 0. Again≤ −z cannot coalesce∈ with the pole − − 0 ≤ ≤− 1 at z = 0 nor with the other turning point z = z2. Case 2: In this case expansions are derived for z lying in domains containing the turning point z , including the interval z + δ z < , but not z or 2 1 ≤ ∞ 1 0. Here a lies in a larger interval than Cases 1a and 1b, namely 0 < a0 a2 a < , and hence u(1 a ) α u(a 1). Note this can include≤ ≤ 1 ∞ − − 0 ≤ ≤ 1 − large negative values of α if a0 < 1. The turning points again cannot coalesce (a = 0) in this case, but z1 can coalesce with z = 0 (since both points are excluded). The more general case 0 a < with 0 z < will require an application of asymptotic expansions≤ valid∞ for a coalescing≤ ∞ turning point and double pole, and also for two coalescing turning points. This will be studied in a subsequent paper. 4 T. M. Dunster et al. The plan of this paper is as follows. In section 2 we consider Case 1a and construct L-G expansions for the functions. In this section a detailed description of the Liouvile transformation is provided. In section 3 we obtain Airy function expansions for Case 1a. In section 4 we obtain L-G and Airy expansions for Case 1b, and in section 5 we likewise do this for Case 2. Finally, in section 6 we present numerical results for the expansions of Cases 1a and 2 for Laguerre polynomials. 2 L-G expansions: Case 1a We make the Liouville transformation z1 ξ = f 1/2(t)dt. (2.1) Zz On explicit integration (and noting that √z z = a2 1 ) we obtain 1 2 − 1 1 ξ = a2 +1 ln a2 +1 z S (z) S(z) max a2, 1 ln(2a) 2 − − − 2 − 2 (2.2) 1 a2 1 + a2 1 S(z) a2 +1 z + a2 1 ln − − − . 2 − ( z ) In this S(z) is given by 2 1/2 S(z)= (z z)(z z) 1/2 = z2 2 a2 +1 z + a2 1 , (2.3) { 1 − 2 − } − − n o where the principal logarithms are taken, and the branches of the square roots are chosen so that S(z) > 0 for <z<z , and is continuous elsewhere −∞ 1 in the plane having a cut on the interval [z1,z2]. This means S(z) < 0 for z <z< . 2 ∞ Thus ξ is real and positive when z (0,z1), and varies continuously in the z plane having cuts along ( , 0] and∈ [z , ). Note that z = z is mapped to −∞ 1 ∞ 1 ξ = 0, and z = z2 i0 (respectively above and below the cut) is mapped to ξ = min a2, 1 πi±. We± find that ξ + as z 0+, such that → ∞ → ξ = 1 a2 1 2 ln a2 1 ln (z) 1 a2 ln(a) + (z). (2.4) 2 − − − − − | | O Furthermore, we have that ξ = 1 z 1 a2 +1 ln(z)+1 + a2 ln(a) + min a2, 1 πi + z−1 , (2.5) 2 − 2 { } O as z in upper half plane. Figs.→ ∞ 1 and 2 depict the z-ξ map for 0 arg(z) π, with corresponding points labeled with the same letters. For the≤ lower≤ half z plane we can use from the Schwarz reflection principle that ξ(z)= ξ (¯z). Uniform asymptotic expansions for Laguerre polynomials 5 H ()0 z ()1 ()z2 A B C D E F G Fig.