[go: up one dir, main page]

14 Legendre and Related FunctionsReal Arguments

§14.15 Uniform Asymptotic Approximations

Contents
  1. §14.15(i) Large μ, Fixed ν
  2. §14.15(ii) Large μ, 0≤ν+12≤(1−δ)⁢μ
  3. §14.15(iii) Large ν, Fixed μ
  4. §14.15(iv) Large ν, 0≤μ≤(1−δ)⁢(ν+12)
  5. §14.15(v) Large ν, (ν+12)⁢δ≤μ≤(ν+12)/δ

§14.15(i) Large μ, Fixed ν

For the interval −1<x<1 with fixed ν, real μ, and arbitrary fixed values of the nonnegative integer J,

14.15.1 𝖯ν−μ⁡(±x)=(1∓x1±x)μ/2⁢(∑j=0J−1(ν+1)j⁢(−ν)jj!⁢Γ⁡(j+1+μ)⁢(1∓x2)j+O⁡(1Γ⁡(J+1+μ)))

as μ→∞, uniformly with respect to x. In other words, the convergent hypergeometric series expansions of 𝖯ν−μ⁡(±x) are also generalized (and uniform) asymptotic expansions as μ→∞, with scale 1/Γ⁡(j+1+μ), j=0,1,2,…; compare §2.1(v).

Provided that μ−ν∉ℤ the corresponding expansions for 𝖯νμ⁡(x) and 𝖰ν∓μ⁡(x) can be obtained from the connection formulas (14.9.7), (14.9.9), and (14.9.10).

For the interval 1<x<∞ the following asymptotic approximations hold when μ→∞, with ν (≥−12) fixed, uniformly with respect to x:

where u is given by (14.12.10). Here I and K are the modified Bessel functions (§10.25(ii)).

For asymptotic expansions and explicit error bounds, see Dunster (2003b) and Gil et al. (2000). See also Temme (2015, Chapter 29).

§14.15(ii) Large μ, 0≤ν+12≤(1−δ)⁢μ

In this and subsequent subsections δ denotes an arbitrary constant such that 0<δ<1.

As μ→∞,

14.15.4 𝖯ν−μ⁡(x)=1Γ⁡(μ+1)⁢(1−α2)−μ/2⁢(1−α1+α)(ν/2)+(1/4)⁢(px)1/2⁢e−μ⁢ρ⁢(1+O⁡(1μ)),

uniformly with respect to x∈(−1,1) and ν+12∈[0,(1−δ)⁢μ], where

14.15.5 α=ν+12μ(<1),
14.15.6 p=x(α2⁢x2+1−α2)1/2,

and

14.15.7 ρ=12⁢ln⁡(1+p1−p)+12⁢α⁢ln⁡(1−α⁢p1+α⁢p).

With the same conditions, the corresponding approximation for 𝖯ν−μ⁡(−x) is obtained by replacing e−μ⁢ρ by eμ⁢ρ on the right-hand side of (14.15.4). Approximations for 𝖯νμ⁡(x) and 𝖰ν∓μ⁡(x) can then be achieved via (14.9.7), (14.9.9), and (14.9.10).

Next,

14.15.8 Pν−μ⁡(x)=(2⁢μπ)1/2⁢1Γ⁡(μ+1)⁢(1−α1+α)(ν/2)+(1/4)⁢(1−α2)−μ/2×(α2+η2α2⁢(x2−1)+1)1/4⁢Kν+12⁡(μ⁢η)⁢(1+O⁡(1μ)),
14.15.9 𝑸νμ⁡(x)=(π2)1/2⁢(eμ)ν+(1/2)⁢(1−α1+α)μ/2⁢(1−α2)−(ν/2)−(1/4)×(α2+η2α2⁢(x2−1)+1)1/4⁢Iν+12⁡(μ⁢η)⁢(1+O⁡(1μ)),

uniformly with respect to x∈(1,∞) and ν+12∈[0,(1−δ)⁢μ]. Here α is again given by (14.15.5), and η is defined implicitly by

14.15.10 α⁢ln⁡((α2+η2)1/2+α)−α⁢ln⁡η−(α2+η2)1/2=12⁢ln⁡((1+α2)⁢x2+1−α2−2⁢x⁢(α2⁢x2−α2+1)1/2(x2−1)⁢(1−α2))+12⁢α⁢ln⁡(α2⁢(2⁢x2−1)+1+2⁢α⁢x⁢(α2⁢x2−α2+1)1/21−α2).

The interval 1<x<∞ is mapped one-to-one to the interval 0<η<∞, with the points x=1 and x=∞ corresponding to η=∞ and η=0, respectively. For asymptotic expansions and explicit error bounds, see Dunster (2003b).

§14.15(iii) Large ν, Fixed μ

For ν→∞ and fixed μ (≥0),

uniformly for θ∈(0,π−δ]. For the Bessel functions J and Y see §10.2(ii), and for the env functions associated with J and Y see §2.8(iv).

For asymptotic expansions and explicit error bounds, see Olver (1997b, Chapter 12, §§12, 13) and Jones (2001). For convergent series expansions see Dunster (2004). See also Temme (2015, Chapter 29).

See also Olver (1997b, pp. 311–313) and §18.15(iii) for a generalized asymptotic expansion in terms of elementary functions for Legendre polynomials Pn⁡(cos⁡θ) as n→∞ with θ fixed.

§14.15(iv) Large ν, 0≤μ≤(1−δ)⁢(ν+12)

As ν→∞,

14.15.15 𝖯ν−μ⁡(x)=β⁢(y−α21−α2−x2)1/4⁢(Jμ⁡((ν+12)⁢y1/2)+O⁡(1ν)⁢env⁡Jμ⁡((ν+12)⁢y1/2)),

uniformly with respect to x∈[0,1) and μ∈[0,(1−δ)⁢(ν+12)]. For α, β, and y see below.

Next,

14.15.17 Pν−μ⁡(x)=β⁢(α2−yx2−1+α2)1/4⁢Iμ⁡((ν+12)⁢|y|1/2)⁢(1+O⁡(1ν)),
14.15.18 𝑸νμ⁡(x)=1β⁢Γ⁡(ν+μ+1)⁢(α2−yx2−1+α2)1/4⁢Kμ⁡((ν+12)⁢|y|1/2)⁢(1+O⁡(1ν)),

uniformly with respect to x∈(1,∞) and μ∈[0,(1−δ)⁢(ν+12)]. In (14.15.15)–(14.15.18)

14.15.19 α=μν+12(<1),
14.15.20 β=eμ⁢(ν−μ+12ν+μ+12)(ν/2)+(1/4)⁢((ν+12)2−μ2)−μ/2,

and the variable y is defined implicitly by

14.15.21 (y−α2)1/2−α⁢arctan⁡((y−α2)1/2α)=arccos⁡(x(1−α2)1/2)−α2⁢arccos⁡((1+α2)⁢x2−1+α2(1−α2)⁢(1−x2)),
x≤(1−α2)1/2, y≥α2,

and

14.15.22 (α2−y)1/2+12⁢α⁢ln⁡|y|−α⁢ln⁡((α2−y)1/2+α)=ln⁡(x+(x2−1+α2)1/2(1−α2)1/2)+α2⁢ln⁡((1−α2)⁢|1−x2|(1+α2)⁢x2−1+α2+2⁢α⁢x⁢(x2−1+α2)1/2),
x≥(1−α2)1/2, y≤α2,

where the inverse trigonometric functions take their principal values (§4.23(ii)). The points x=(1−α2)1/2, x=1, and x=∞ are mapped to y=α2, y=0, and y=−∞, respectively. The interval 0≤x<∞ is mapped one-to-one to the interval −∞<y≤y0, where y=y0 is the (positive) solution of (14.15.21) when x=0.

For asymptotic expansions and explicit error bounds, see Boyd and Dunster (1986).

§14.15(v) Large ν, (ν+12)⁢δ≤μ≤(ν+12)/δ

Here we introduce the envelopes of the parabolic cylinder functions U⁡(−c,x), U¯⁡(−c,x), which are defined in §12.2. For U⁡(−c,x) or U¯⁡(−c,x), with c and x nonnegative,

14.15.23 env⁢U⁡(−c,x) ={(U2⁡(−c,x)+U¯2⁡(−c,x))1/2,0≤x≤Xc,2⁢U⁡(−c,x),Xc≤x<∞,
env⁢U¯⁡(−c,x) ={(U2⁡(−c,x)+U¯2⁡(−c,x))1/2,0≤x≤Xc,2⁢U¯⁡(−c,x),Xc≤x<∞,

where x=Xc denotes the largest positive root of the equation U⁡(−c,x)=U¯⁡(−c,x).

As ν→∞,

14.15.24 𝖯ν−μ⁡(x)=1(ν+12)1/4⁢2(ν+μ)/2⁢Γ⁡(12⁢ν+12⁢μ+34)⁢(ζ2−α2x2−a2)1/4×(U⁡(μ−ν−12,(2⁢ν+1)1/2⁢ζ)+O⁡(ν−2/3)⁢env⁢U⁡(μ−ν−12,(2⁢ν+1)1/2⁢ζ)),
14.15.25 𝖰ν−μ⁡(x)=π(ν+12)1/4⁢2(ν+μ+2)/2⁢Γ⁡(12⁢ν+12⁢μ+34)⁢(ζ2−α2x2−a2)1/4×(U¯⁡(μ−ν−12,(2⁢ν+1)1/2⁢ζ)+O⁡(ν−2/3)⁢env⁢U¯⁡(μ−ν−12,(2⁢ν+1)1/2⁢ζ)),

uniformly with respect to x∈[0,1) and μ∈[δ⁢(ν+12),ν+12]. Here

14.15.26 a =((ν+μ+12)⁢|ν−μ+12|)1/2ν+12,
α =(2⁢|ν−μ+12|ν+12)1/2,

and the variable ζ is defined implicitly by

14.15.27 12⁢ζ⁢(ζ2−α2)1/2−12⁢α2⁢arccosh⁡(ζα)=(1−a2)1/2⁢arctanh⁡(1x⁢(x2−a21−a2)1/2)−arccosh⁡(xa),
a≤x<1, α≤ζ<∞,

and

14.15.28 12⁢α2⁢arcsin⁡(ζα)+12⁢ζ⁢(α2−ζ2)1/2=arcsin⁡(xa)−(1−a2)1/2⁢arctan⁡(x⁢(1−a2a2−x2)1/2),
−a≤x≤a, −α≤ζ≤α,

when a>0, and

14.15.29 ζ2=−ln⁡(1−x2),
−1<x<1,

when a=0. The inverse hyperbolic and trigonometric functions take their principal values (§§4.23(ii), 4.37(ii)).

When a>0 the interval −a≤x<1 is mapped one-to-one to the interval −α≤ζ<∞, with the points x=−a, x=a, and x=1 corresponding to ζ=−α, ζ=α, and ζ=∞, respectively. When a=0 the interval −1<x<1 is mapped one-to-one to the interval −∞<ζ<∞, with the points x=−1, 0, and 1 corresponding to ζ=−∞, 0, and ∞, respectively.

Next, as ν→∞,

14.15.30 𝖯ν−μ⁡(x)=1(ν+12)1/4⁢2(ν+μ)/2⁢Γ⁡(12⁢ν+12⁢μ+34)⁢(ζ2+α2x2+a2)1/4×U⁡(μ−ν−12,(2⁢ν+1)1/2⁢ζ)⁢(1+O⁡(ν−1⁢ln⁡ν)),

uniformly with respect to x∈(−1,1) and μ∈[ν+12,(1/δ)⁢(ν+12)]. Here ζ is defined implicitly by

14.15.31 12⁢ζ⁢(ζ2+α2)1/2+12⁢α2⁢arcsinh⁡(ζα)=(1+a2)1/2⁢arctanh⁡(x⁢(1+a2x2+a2)1/2)−arcsinh⁡(xa),
−1<x<1, −∞<ζ<∞,

when a>0, which maps the interval −1<x<1 one-to-one to the interval −∞<ζ<∞: the points x=−1 and x=1 correspond to ζ=−∞ and ζ=∞, respectively. When a=0 (14.15.29) again applies. (The inverse hyperbolic functions again take their principal values.)

Since (14.15.30) holds for negative x, corresponding approximations for 𝖰ν∓μ⁡(x), uniformly valid in the interval −1<x<1, can be obtained from (14.9.9) and (14.9.10).

For error bounds and other extensions see Olver (1975b).