Digital Library of Mathematical Functions
About the Project
NIST
14 Legendre and Related FunctionsReal Arguments

§14.15 Uniform Asymptotic Approximations

Contents

§14.15(i) Large \mu, Fixed \nu

For the interval -1<x<1 with fixed \nu, real \mu, and arbitrary fixed values of the nonnegative integer J,

as \mu\to\infty, uniformly with respect to x. In other words, the convergent hypergeometric series expansions of \mathop{\mathsf{P}^{{-\mu}}_{{\nu}}\/}\nolimits\!\left(\pm x\right) are also generalized (and uniform) asymptotic expansions as \mu\to\infty, with scale \ifrac{1}{\mathop{\Gamma\/}\nolimits\!\left(j+1+\mu\right)}, j=0,1,2,\dots; compare §2.1(v).

Provided that \mu-\nu\notin\Integer the corresponding expansions for \mathop{\mathsf{P}^{{\mu}}_{{\nu}}\/}\nolimits\!\left(x\right) and \mathop{\mathsf{Q}^{{\mp\mu}}_{{\nu}}\/}\nolimits\!\left(x\right) can be obtained from the connection formulas (14.9.7), (14.9.9), and (14.9.10).

For the interval 1<x<\infty the following asymptotic approximations hold when \mu\to\infty, with \nu (\geq-\frac{1}{2}) fixed, uniformly with respect to x:

where u is given by (14.12.10). Here \mathop{I\/}\nolimits and \mathop{K\/}\nolimits are the modified Bessel functions (§10.25(ii)).

For asymptotic expansions and explicit error bounds, see Dunster (2003b) and Gil et al. (2000).

§14.15(ii) Large \mu, 0\leq\nu+\frac{1}{2}\leq(1-\delta)\mu

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

As \mu\to\infty,

uniformly with respect to x\in(-1,1) and \nu+\tfrac{1}{2}\in[0,(1-\delta)\mu], where

14.15.5\alpha=\frac{\nu+\frac{1}{2}}{\mu}\,(<1),
14.15.6p=\frac{x}{\left(\alpha^{2}x^{2}+1-\alpha^{2}\right)^{{1/2}}},

and

14.15.7\rho=\frac{1}{2}\mathop{\ln\/}\nolimits\!\left(\frac{1+p}{1-p}\right)+\frac{1}%
{2}\alpha\mathop{\ln\/}\nolimits\!\left(\frac{1-\alpha p}{1+\alpha p}\right).

With the same conditions, the corresponding approximation for \mathop{\mathsf{P}^{{-\mu}}_{{\nu}}\/}\nolimits\!\left(-x\right) is obtained by replacing e^{{-\mu\rho}} by e^{{\mu\rho}} on the right-hand side of (14.15.4). Approximations for \mathop{\mathsf{P}^{{\mu}}_{{\nu}}\/}\nolimits\!\left(x\right) and \mathop{\mathsf{Q}^{{\mp\mu}}_{{\nu}}\/}\nolimits\!\left(x\right) can then be achieved via (14.9.7), (14.9.9), and (14.9.10).

Next,

uniformly with respect to x\in(1,\infty) and \nu+\tfrac{1}{2}\in[0,(1-\delta)\mu]. Here \alpha is again given by (14.15.5), and \eta is defined implicitly by

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

§14.15(iii) Large \nu, Fixed \mu

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 Olver (1997b, pp. 311–313) and §18.15(iii) for a generalized asymptotic expansion in terms of elementary functions for Legendre polynomials \mathop{P_{{n}}\/}\nolimits\!\left(\mathop{\cos\/}\nolimits\theta\right) as n\to\infty with \theta fixed.

§14.15(iv) Large \nu, 0\leq\mu\leq(1-\delta)(\nu+\frac{1}{2})

Next,

uniformly with respect to x\in(1,\infty) and \mu\in[0,(1-\delta)(\nu+\frac{1}{2})]. In (14.15.15)–(14.15.18)

14.15.19\alpha=\frac{\mu}{\nu+\frac{1}{2}}\,(<1),

and the variable y is defined implicitly by

14.15.21\left(y-\alpha^{2}\right)^{{1/2}}-\alpha\mathop{\mathrm{arctan}\/}\nolimits\!%
\left(\frac{\left(y-\alpha^{2}\right)^{{1/2}}}{\alpha}\right)=\mathop{\mathrm{%
arccos}\/}\nolimits\!\left(\frac{x}{\left(1-\alpha^{2}\right)^{{1/2}}}\right)-%
\frac{\alpha}{2}\mathop{\mathrm{arccos}\/}\nolimits\!\left(\frac{\left(1+%
\alpha^{2}\right)x^{2}-1+\alpha^{2}}{\left(1-\alpha^{2}\right)\left(1-x^{2}%
\right)}\right),x\leq\left(1-\alpha^{2}\right)^{{1/2}}, y\geq\alpha^{2},

and

where the inverse trigonometric functions take their principal values (§4.23(ii)). The points x=\left(1-\alpha^{2}\right)^{{1/2}}, x=1, and x=\infty are mapped to y=\alpha^{2}, y=0, and y=-\infty, respectively. The interval 0\leq x<\infty is mapped one-to-one to the interval -\infty<y\leq y_{0}, where y=y_{0} 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 \nu, (\nu+\frac{1}{2})\delta\leq\mu\leq(\nu+\frac{1}{2})/\delta

Here we introduce the envelopes of the parabolic cylinder functions \mathop{U\/}\nolimits\!\left(-c,x\right), \mathop{\overline{U}\/}\nolimits\!\left(-c,x\right), which are defined in §12.2. For f(x)=\mathop{U\/}\nolimits\!\left(-c,x\right) or \mathop{\overline{U}\/}\nolimits\!\left(-c,x\right), with c and x nonnegative,

where x=X_{c} denotes the largest positive root of the equation \mathop{U\/}\nolimits\!\left(-c,x\right)=\mathop{\overline{U}\/}\nolimits\!%
\left(-c,x\right).

As \nu\to\infty,

uniformly with respect to x\in[0,1) and \mu\in[\delta(\nu+\frac{1}{2}),\nu+\frac{1}{2}]. Here

14.15.26
a=\frac{\left(\left(\nu+\mu+\frac{1}{2}\right)\left|\nu-\mu+\frac{1}{2}\right|%
\right)^{{1/2}}}{\nu+\frac{1}{2}},
\alpha=\left(\frac{2\left|\nu-\mu+\frac{1}{2}\right|}{\nu+\frac{1}{2}}\right)^%
{{1/2}},

and the variable \zeta is defined implicitly by

and

when a>0, and

14.15.29\zeta^{2}=-\mathop{\ln\/}\nolimits\!\left(1-x^{2}\right),-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\leq x<1 is mapped one-to-one to the interval -\alpha\leq\zeta<\infty, with the points x=-a, x=a, and x=1 corresponding to \zeta=-\alpha, \zeta=\alpha, and \zeta=\infty, respectively. When a=0 the interval -1<x<1 is mapped one-to-one to the interval -\infty<\zeta<\infty, with the points x=-1, 0, and 1 corresponding to \zeta=-\infty, 0, and \infty, respectively.

Next, as \nu\to\infty,

uniformly with respect to x\in(-1,1) and \mu\in[\nu+\frac{1}{2},(1/\delta)(\nu+\frac{1}{2})]. Here \zeta is defined implicitly by

when a>0, which maps the interval -1<x<1 one-to-one to the interval -\infty<\zeta<\infty: the points x=-1 and x=1 correspond to \zeta=-\infty and \zeta=\infty, 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 \mathop{\mathsf{Q}^{{\mp\mu}}_{{\nu}}\/}\nolimits\!\left(x\right), 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).