[go: up one dir, main page]

14 Legendre and Related FunctionsReal Arguments

§14.20 Conical (or Mehler) Functions

Contents
  1. §14.20(i) Definitions and Wronskians
  2. §14.20(ii) Graphics
  3. §14.20(iii) Behavior as x→1
  4. §14.20(iv) Integral Representation
  5. §14.20(v) Trigonometric Expansion
  6. §14.20(vi) Generalized Mehler–Fock Transformation
  7. §14.20(vii) Asymptotic Approximations: Large τ, Fixed μ
  8. §14.20(viii) Asymptotic Approximations: Large τ, 0≤μ≤A⁢τ
  9. §14.20(ix) Asymptotic Approximations: Large μ, 0≤τ≤A⁢μ
  10. §14.20(x) Zeros and Integrals

§14.20(i) Definitions and Wronskians

Throughout §14.20 we assume that ν=−12+i⁢τ, with μ≥0 and τ≥0. (14.2.2) takes the form

14.20.1 (1−x2)⁢d2wdx2−2⁢x⁢dwdx−(τ2+14+μ21−x2)⁢w=0.

Solutions are known as conical or Mehler functions. For −1<x<1 and τ>0, a numerically satisfactory pair of real conical functions is 𝖯−12+i⁢τ−μ⁡(x) and 𝖯−12+i⁢τ−μ⁡(−x).

Another real-valued solution 𝖰^−12+i⁢τ−μ⁡(x) of (14.20.1) was introduced in Dunster (1991). This is defined by

14.20.2 𝖰^−12+i⁢τ−μ⁡(x)=ℜ⁡(eμ⁢π⁢i⁢𝖰−12+i⁢τ−μ⁡(x))−12⁢π⁢sin⁡(μ⁢π)⁢𝖯−12+i⁢τ−μ⁡(x).

Equivalently,

14.20.3 𝖰^−12+i⁢τ−μ⁡(x)=π⁢e−τ⁢π⁢sin⁡(μ⁢π)⁢sinh⁡(τ⁢π)2⁢(cosh2⁡(τ⁢π)−sin2⁡(μ⁢π))⁢𝖯−12+i⁢τ−μ⁡(x)+π⁢(e−τ⁢π⁢cos2⁡(μ⁢π)+sinh⁡(τ⁢π))2⁢(cosh2⁡(τ⁢π)−sin2⁡(μ⁢π))⁢𝖯−12+i⁢τ−μ⁡(−x).

𝖰^−12+i⁢τ−μ⁡(x) exists except when μ=12,32,… and τ=0; compare §14.3(i). It is an important companion solution to 𝖯−12+i⁢τ−μ⁡(x) when τ is large; compare §§14.20(vii), 14.20(viii), and 10.25(iii).

14.20.4 𝒲⁡{𝖯−12+i⁢τ−μ⁡(x),𝖯−12+i⁢τ−μ⁡(−x)}=2|Γ⁡(μ+12+i⁢τ)|2⁢(1−x2).
14.20.5 𝒲⁡{𝖯−12+i⁢τ−μ⁡(x),𝖰^−12+i⁢τ−μ⁡(x)}=π⁢(e−τ⁢π⁢cos2⁡(μ⁢π)+sinh⁡(τ⁢π))|Γ⁡(μ+12+i⁢τ)|2⁢(cosh2⁡(τ⁢π)−sin2⁡(μ⁢π))⁢(1−x2),

provided that 𝖰^−12+i⁢τ−μ⁡(x) exists.

Lastly, for the range 1<x<∞, P−12+i⁢τ−μ⁡(x) is a real-valued solution of (14.20.1); in terms of Q−12±i⁢τμ⁡(x) (which are complex-valued in general):

§14.20(ii) Graphics

See accompanying text
Figure 14.20.1: 𝖯−12+i⁢τ0⁡(x), τ=0,1,2,4,8. Magnify
See accompanying text
Figure 14.20.2: 𝖰^−12+i⁢τ0⁡(x), τ=0,12,1,2,4. Magnify
See accompanying text
Figure 14.20.3: 𝖯−12+i⁢τ−1/2⁡(x), τ=0,1,2,4,8. Magnify
See accompanying text
Figure 14.20.4: 𝖰^−12+i⁢τ−1/2⁡(x), τ=12,1,2,4. (This function does not exist when τ=0.) Magnify
See accompanying text
Figure 14.20.5: 𝖯−12+i⁢τ−1⁡(x), τ=0,1,2,4,8. Magnify
See accompanying text
Figure 14.20.6: 𝖰^−12+i⁢τ−1⁡(x), τ=0,12,1,2,4. Magnify
See accompanying text
Figure 14.20.7: 𝖯−12+i⁢τ−2⁡(x),τ=0,1,2,4,8. Magnify
See accompanying text
Figure 14.20.8: 𝖰^−12+i⁢τ−2⁡(x), τ=0,12,1,2,4. Magnify

§14.20(iii) Behavior as x→1

The behavior of 𝖯−12+i⁢τ−μ⁡(±x) as x→1− is given in §14.8(i). For μ>0 and x→1−,

14.20.7 𝖰^−12+i⁢τμ⁡(x)∼12⁢Γ⁡(μ)⁢(21−x)μ/2,
14.20.8 𝖰^−12+i⁢τ−μ⁡(x)∼π⁢Γ⁡(μ)⁢(e−τ⁢π⁢cos2⁡(μ⁢π)+sinh⁡(τ⁢π))2⁢(cosh2⁡(τ⁢π)−sin2⁡(μ⁢π))⁢|Γ⁡(μ+12+i⁢τ)|2⁢(21−x)μ/2.

§14.20(iv) Integral Representation

When 0<θ<π,

§14.20(v) Trigonometric Expansion

14.20.10 𝖯−12+i⁢τ⁡(cos⁡θ)=1+4⁢τ2+1222⁢sin2⁡(12⁢θ)+(4⁢τ2+12)⁢(4⁢τ2+32)22⋅42⁢sin4⁡(12⁢θ)+⋯,
0≤θ≤π.

From (14.20.9) or (14.20.10) it is evident that 𝖯−12+i⁢τ⁡(cos⁡θ) is positive for real θ.

§14.20(vi) Generalized Mehler–Fock Transformation

14.20.11 f⁡(τ)=τπ⁢sinh⁡(τ⁢π)⁢Γ⁡(12−μ+i⁢τ)⁢Γ⁡(12−μ−i⁢τ)⁢∫1∞P−12+i⁢τμ⁡(x)⁢g⁡(x)⁢dx,

where

14.20.12 g⁡(x)=∫0∞P−12+i⁢τμ⁡(x)⁢f⁡(τ)⁢dτ.

Special cases:

§14.20(vii) Asymptotic Approximations: Large τ, Fixed μ

For τ→∞ and fixed μ,

14.20.15 𝖯−12+i⁢τ−μ⁡(cos⁡θ) =1τμ⁢(θsin⁡θ)1/2⁢Iμ⁡(τ⁢θ)⁢(1+O⁡(1/τ)),
14.20.16 𝖰^−12+i⁢τ−μ⁡(cos⁡θ) =1τμ⁢(θsin⁡θ)1/2⁢Kμ⁡(τ⁢θ)⁢(1+O⁡(1/τ)),

uniformly for θ∈(0,π−δ], where I and K are the modified Bessel functions (§10.25(ii)) and δ is an arbitrary constant such that 0<δ<π. For asymptotic expansions and explicit error bounds, see Olver (1997b, pp. 473–474). See also Žurina and Karmazina (1966).

§14.20(viii) Asymptotic Approximations: Large τ, 0≤μ≤A⁢τ

In this subsection and §14.20(ix), A and δ denote arbitrary constants such that A>0 and 0<δ<2.

As τ→∞,

14.20.17 𝖯−12+i⁢τ−μ⁡(x)=σ⁡(μ,τ)⁢(α2+η1+α2−x2)1/4⁢Iμ⁡(τ⁢η1/2)⁢(1+O⁡(1/τ)),
14.20.18 𝖰^−12+i⁢τ−μ⁡(x)=σ⁡(μ,τ)⁢(α2+η1+α2−x2)1/4⁢Kμ⁡(τ⁢η1/2)⁢(1+O⁡(1/τ)),

uniformly for x∈[−1+δ,1) and μ∈[0,A⁢τ]. Here

14.20.19 α=μ/τ,
14.20.20 σ⁡(μ,τ)=exp⁡(μ−τ⁢arctan⁡α)(μ2+τ2)μ/2.

The variable η is defined implicitly by

14.20.21 (α2+η)1/2+12⁢α⁢ln⁡η−α⁢ln⁡((α2+η)1/2+α)=arccos⁡(x(1+α2)1/2)+α2⁢ln⁡(1+α2+(α2−1)⁢x2−2⁢α⁢x⁢(1+α2−x2)1/2(1+α2)⁢(1−x2)),

where the inverse trigonometric functions take their principal values. The interval −1<x<1 is mapped one-to-one to the interval 0<η<∞, with the points x=−1 and x=1 corresponding to η=∞ and η=0, respectively.

§14.20(ix) Asymptotic Approximations: Large μ, 0≤τ≤A⁢μ

As μ→∞,

14.20.22 𝖯−12+i⁢τ−μ⁡(x)=exp⁡(μ⁢β⁢arctan⁡β)Γ⁡(μ+1)⁢(1+β2)μ/2⁢e−μ⁢ρ(1+β2−x2⁢β2)1/4⁢(1+O⁡(1μ)),

uniformly for x∈(−1,1) and τ∈[0,A⁢μ]. Here

14.20.23 β=τ/μ,

and the variable ρ is defined by

14.20.24 ρ=12⁢ln⁡((1−β2)⁢x2+1+β2+2⁢x⁢(1+β2−β2⁢x2)1/21−x2)+β⁢arctan⁡(β⁢x1+β2−β2⁢x2)−12⁢ln⁡(1+β2),

with the inverse tangent taking its principal value. The interval −1<x<1 is mapped one-to-one to the interval −∞<ρ<∞, with the points x=−1, x=0, and x=1 corresponding to ρ=−∞, ρ=0, and ρ=∞, respectively.

With the same conditions, the corresponding approximation for 𝖯−12+i⁢τ−μ⁡(−x) is obtainable by replacing e−μ⁢ρ by eμ⁢ρ on the right-hand side of (14.20.22). Approximations for 𝖯−12+i⁢τμ⁡(x) and 𝖰^−12+i⁢τ−μ⁡(x) can then be achieved via (14.9.7) and (14.20.3).

For extensions to complex arguments (including the range 1<x<∞), asymptotic expansions, and explicit error bounds, see Dunster (1991). For the case of purely imaginary order and argument see Dunster (2013).

§14.20(x) Zeros and Integrals

For zeros of 𝖯−12+i⁢τ⁡(x) see Hobson (1931, §237).

For integrals with respect to τ involving 𝖯−12+i⁢τ⁡(x), see Prudnikov et al. (1990, pp. 218–228).