ASYMPTOTIC EXPANSIONS FOR THE INCOMPLETE GAMMA FUNCTION IN THE TRANSITION REGIONS GERGO} NEMES AND ADRI B. OLDE DAALHUIS Abstract. We construct asymptotic expansions for the normalised incom- plete gamma function Q(a; z) = Γ(a; z)=Γ(a) that are valid in the transition regions, including the case z ≈ a, and have simple polynomial coefficients. For Bessel functions, these type of expansions are well known, but for the nor- malised incomplete gamma function they were missing from the literature. A detailed historical overview is included. We also derive an asymptotic expan- sion for the corresponding inverse problem, which has importance in proba- bility theory and mathematical statistics. The coefficients in this expansion are again simple polynomials, and therefore its implementation is straightfor- ward. As a byproduct, we give the first complete asymptotic expansion as a ! −∞ of the unique negative zero of the regularised incomplete gamma function γ∗(a; x). 1. Introduction and main results The normalised incomplete gamma function Q(a; z) = Γ(a; z)=Γ(a) is one of the most widely used special functions of two variables. It is used in constructing gamma distributions [9, Ch. 17], which appear naturally in the theory associated with normally distributed random variables. In many applications, Poisson random variables are used in which the Poisson rate is not fixed, and the function Q(a; z) plays an important role in these cases (see, for example, [6]). Consequently, efficient and accurate approximations for the incomplete gamma function are essential. In many applications the inverse problem also arises naturally, that is, the problem of determining the solution x, real and positive, of the equation Q(a; x) = q for given 0 < q < 1 and a > 0. In the important papers [27, 31], Temme established asymptotic expansions for the normalised incomplete gamma function as a ! 1, that are valid uniformly for jzj 2 [0; 1). The asymptotic expansions of the corresponding inverse functions are discussed in [29]. The large region of validity is important in applications, and many papers in physics, applied statistics, network and control theory, and so on, refer to these uniform asymptotic expansions (see, e.g., [11, 14, 21, 22, 26, 33]). However, applying Temme's expansions in the neighbourhood of the (most in- teresting) point z = a can be difficult because its coefficients possess a removable arXiv:1803.07841v2 [math.CA] 5 Jul 2018 singularity at this point. To overcome this difficulty, Temme [27] gave power series expansions for these coefficients. Having to use local power series expansions for the coefficients in uniform asymptotic expansions, is not a very elegant solution. In this paper, we propose an alternative approach to the problem. We will construct what we call transition region expansions, and provide full details of the inversion of these new expansions. These are expansions that are valid in the regions in which Q(a; z) changes dramatically, and their coefficients are 2010 Mathematics Subject Classification. 33B20, 41A60. Key words and phrases. incomplete gamma functions, asymptotic expansions, transition regions. The authors' research was supported by a research grant (GRANT11863412/70NANB15H221) from the National Institute of Standards and Technology. 1 2 G. NEMES AND A. B. OLDE DAALHUIS polynomials satisfying simple recurrence relations. The region of validity overlaps with those of the non-uniform \outer" expansions. Furthermore, the coefficients of their inversions are simple polynomials, whose computation and implementation are straightforward. It is surprising to us that these transition region expansions for the normalised incomplete gamma function have not yet been discussed in the literature, given the fact that expansions of similar type for Bessel functions are well known (see [18, x10.19(iii)]). What our new expansions and the transition region expansions for the Bessel functions have in common is that both mimic the corresponding uniform expansions; compare the similarities between (1.1) and (2.6), and between [18, Eq. 10.19.8] and [18, Eq. 10.20.4]. The ideas in this paper should also be applicable to other cumulative distribution functions using the results of [28]. The proofs of our new results rely heavily on the uniform asymptotic expansions by Temme. Thus, we shall provide full details of his expansions, and even include new recurrence relations for the coefficients in the local power series expansions mentioned above (see AppendixB). For more information about existing results in the literature and their connection with those we prove in this paper, the reader is referred to the historical overview in Section2. The following theorem is the main result of the paper. Theorem 1.1. The normalised incomplete gamma function admits the asymptotic expansions 2 1 1 1 − 1 1 τ X Cn(τ) (1.1) Q a; a + τa 2 ∼ erfc 2 2 τ + p exp − 2 2 an=2 2πa n=0 and ±πia e 1 ±πi 1 − 1 Q −a; (a + τa 2 )e ∼ ± erfc ±2 2 iτ 2i sin(πa) 2 (1.2) 2 1 n i τ X (−i) Cn(iτ) − p exp ; 2 an=2 2πa n=0 as a ! 1 in the sector j arg aj ≤ π − δ < π, uniformly with respect to bounded complex values of τ. Here, erfc denotes the complementary error function [18, Eq. 7.2.2]. The coefficients Cn(τ) are polynomials in τ of degree 3n + 2 and satisfy 1 1 C (τ) = τ 2 − ; 0 3 3 0 00 2 2 0 00 Cn(τ) + τCn(τ) − Cn(τ) = τ(τ − 2)Cn−1(τ) − (2τ − 1)Cn−1(τ) + τCn−1(τ) for n ≥ 1. In addition, the even- and odd-order polynomials are even and odd functions, respectively. The explicit forms of the next two coefficients in (1.1) are as follows: 1 11 1 1 29 133 23 1 C (τ) = τ 5 − τ 3 + τ; C (τ) = τ 8 − τ 6 + τ 4 − τ 2 − : 1 18 36 12 2 162 324 540 540 540 For higher coefficients and for an algorithm to generate them, see the end of Section 3. The most common case in applications is when both variables of Q(a; z) are real and positive. It is possible to combine our transition region expansion (1.1), with the outer expansions that exist in the literature and cover the complete range τ > −a1=2 (or equivalently, the range z > 0). First we note that in Theorem 1.1 we assume that τ is bounded. We can weaken this restriction as follows. Since µ Cn(τ) is a polynomial of degree 3n + 2 it follows that we can allow τ = O (jaj ) as 1 a ! 1, as long as µ < 6 . See the proposition below. In Section2, we will give ASYMPTOTICS FOR THE INCOMPLETE GAMMA FUNCTION 3 the expansions (2.1) and (2.2). These are also asymptotic expansions with simple coefficients, and with the notation of this section, these two outer expansions are valid for jτj ≥ aµ as a ! +1, with any µ > 0 (see AppendixA). Therefore, by combining these three expansions, we can cover the whole range τ > −a1=2. Proposition 1.2. The expansions (1.1) and (1.2) are still valid in the sense of gen- eralised asymptotic expansions (see, e.g., [18, x2.1(v)]) provided that τ = O (jajµ) 1 as a ! 1 with some constant µ < 6 . In the inverse problem we are trying to find x = x(a; q) such that (1.3) Q(a; x) = q holds, where 0 < q < 1 and a > 0 are given. The inversion of the normalised incom- plete gamma function is an important topic in mathematical statistics and probabil- ity theory, in particular for computing percentage points of the χ2-distribution (see, for example, [8]). A suitable inversion of the expansion (1.1) yields the following result. Theorem 1.3. For any 0 < q < 1, define τ0 to be the unique real root of the equation 1 1 − 2 (1.4) q = 2 erfc 2 τ0 : Then the inverse function x = x(a; q) that satisfies the equation (1.3) has the as- ymptotic expansion 1 1 X dn(τ0) (1.5) x(a; q) ∼ a + τ a 2 + ; 0 an=2 n=0 as a ! +1, provided that q is bounded away from 0 and 1. The coefficients dn(τ0) are polynomials in τ0 of degree n + 2. In addition, the even- and odd-order polynomials are even and odd functions, respectively. The first few coefficients in (1.5) are as follows: 1 1 1 7 1 7 8 d (τ ) = τ 2 − ; d (τ ) = τ 3 − τ ; d (τ ) = − τ 4 − τ 2 + : 0 0 3 0 3 1 0 36 0 36 0 2 0 270 0 810 0 405 For higher coefficients and for an algorithm to generate them, see Section5. In Theorem 1.3, we assume that q is bounded away from 0 and 1. In practice, these bounds can be very small. It follows from (1.4) that τ0 ! +1 as q ! 0 and τ0 ! −∞ as q ! 1. Since 2 1 1 τ 1 − 2 0 2 erfc 2 τ0 ∼ p exp − 2πτ0 2 as τ0 ! +1, and 2 1 1 τ 1 − 2 0 2 erfc 2 τ0 ∼ 1 + p exp − 2πτ0 2 p as τ0 ! −∞, we can infer that jτ0j ∼ −2 log(q(1 − q)) as q ! 0 or 1. Since the dn(τ0)'s are polynomials in τ0 of degree n + 2, for (1.5) to be an asymptotic 1=2 expansion we need τ0 = o jaj as a ! 1. Hence, we need log(q(1−q)) = o(jaj), i.e., q can be sub-exponentially close in jaj to 0 or 1.