[go: up one dir, main page]

8 Incomplete Gamma and Related FunctionsIncomplete Gamma Functions

§8.11 Asymptotic Approximations and Expansions

Contents
  1. §8.11(i) Large z, Fixed a
  2. §8.11(ii) Large a, Fixed z
  3. §8.11(iii) Large a, Fixed z/a
  4. §8.11(iv) Large a, Bounded (x−a)/(2⁢a)12
  5. §8.11(v) Other Approximations

§8.11(i) Large z, Fixed a

Define

8.11.1 uk=(−1)k⁢(1−a)k=(a−1)⁢(a−2)⁢⋯⁢(a−k),
8.11.2 Γ⁡(a,z)=za−1⁢e−z⁢(∑k=0n−1ukzk+Rn⁡(a,z)),
n=1,2,….

Then as z→∞ with a and n fixed

where δ denotes an arbitrary small positive constant.

If a is real and z (=x) is positive, then Rn⁡(a,x) is bounded in absolute value by the first neglected term un/xn and has the same sign provided that n≥a−1. For bounds on Rn⁡(a,z) when a is real and z is complex see Olver (1997b, pp. 109–112). For an exponentially-improved asymptotic expansion (§2.11(iii)) see Olver (1991a).

§8.11(ii) Large a, Fixed z

8.11.4 γ⁡(a,z)=za⁢e−z⁢∑k=0∞zk(a)k+1,
a≠0,−1,−2,….

This expansion is absolutely convergent for all finite z, and it can also be regarded as a generalized asymptotic expansion (§2.1(v)) of γ⁡(a,z) as a→∞ in |ph⁡a|≤π−δ.

Also,

8.11.5 P⁡(a,z)∼za⁢e−zΓ⁡(1+a)∼(2⁢π⁢a)−12⁢ea−z⁢(z/a)a,
a→∞, |ph⁡a|≤π−δ.

§8.11(iii) Large a, Fixed z/a

If z=λ⁢a, with λ fixed, then as a→∞

8.11.6 γ⁡(a,z)∼−za⁢e−z⁢∑k=0∞(−a)k⁢bk⁡(λ)(z−a)2⁢k+1,
0<λ<1, |ph⁡a|≤π2−δ.
8.11.7 Γ⁡(a,z)∼za⁢e−z⁢∑k=0∞(−a)k⁢bk⁡(λ)(z−a)2⁢k+1,
λ>1, |ph⁡a|≤3⁢π2−δ.

where

8.11.8 b0⁡(λ) =1,
b1⁡(λ) =λ,
b2⁡(λ) =λ⁢(2⁢λ+1),

and for k=1,2,…,

8.11.9 bk⁡(λ)=λ⁢(1−λ)⁢bk−1′⁡(λ)+(2⁢k−1)⁢λ⁢bk−1⁡(λ).

Sharp error bounds and an exponentially-improved extension for (8.11.7) can be found in Nemes (2016). This reference also contains explicit formulas for bk⁡(λ) in terms of Stirling numbers and for the case λ>1 an asymptotic expansion for bk⁡(λ) as k→∞.

The expansion (8.11.7) also applies when a is replaced by −a, λ<0 and |ph⁢a|≤3⁢π2−ω−δ with ω=ph⁢(−λ+ln⁡(−λ)+π⁢i), 0<ω<π. For error bounds and an exponentially-improved extension for this later expansion, see Nemes (2015c). In the case that a=n, a positive integer, the z-region of validity of (8.11.7) is discussed in Ameur and Cronvall (2023).

§8.11(iv) Large a, Bounded (x−a)/(2⁢a)12

If x=a+(2⁢a)12⁢y and a→+∞, then

8.11.11 γ∗⁡(1−a,−x)=xa−1⁢(−cos⁡(π⁢a)+sin⁡(π⁢a)π⁢(2⁢π⁢F⁡(y)+23⁢2⁢πa⁢(1−y2))⁢ey2+O⁡(a−1)),

in both cases uniformly with respect to bounded real values of y. For Dawson’s integral F⁡(y) see §7.2(ii). See Tricomi (1950b) for these approximations, together with higher terms and extensions to complex variables. For related expansions involving Hermite polynomials see Pagurova (1965).

§8.11(v) Other Approximations

As z→∞,

8.11.12 Γ⁡(z,z)∼zz−1⁢e−z⁢(π2⁢z12−13+2⁢π24⁢z12−4135⁢z+2⁢π576⁢z32+82835⁢z2+…),
|ph⁡z|≤2⁢π−δ.

For sharp error bounds and an exponentially-improved extension, see Nemes (2016). This reference also contains explicit formulas for the coefficients in terms of Stirling numbers.

For the function en⁡(z) defined by (8.4.11),

8.11.13 limn→∞en⁡(n⁢x)en⁢x={0,x>1,12,x=1,1,0≤x<1.

With x=1, an asymptotic expansion of en⁡(n⁢x)/en⁢x follows from (8.11.14) and (8.11.16).

If Sn⁡(x) is defined by

8.11.14 en⁢x=en⁡(n⁢x)+(n⁢x)nn!⁢Sn⁡(x),

then

8.11.15 Sn⁡(x)=γ⁡(n+1,n⁢x)(n⁢x)n⁢e−n⁢x.

As n→∞

8.11.16 Sn⁡(1)−12⁢n!⁢ennn∼−23+4135⁢n−1−82835⁢n−2−168505⁢n−3+…,
8.11.17 Sn⁡(−1)∼−12+18⁢n−1+132⁢n−2−1128⁢n−3−13512⁢n−4+….

Also,

8.11.18 Sn⁡(x)∼∑k=0∞dk⁡(x)⁢n−k,
n→∞,

uniformly for x∈(−∞,1−δ], with

8.11.19 d0⁡(x) =x/(1−x),
dk⁡(x) =(−1)k⁢bk⁡(x)(1−x)2⁢k+1,
k=1,2,3,…,

and bk⁡(x) as in §8.11(iii).

For (8.11.18) and extensions to complex values of x see Buckholtz (1963). For a uniformly valid expansion for n→∞ and x∈[δ,1], see Wong (1973b).