Digital Library of Mathematical Functions
About the Project
NIST
10 Bessel FunctionsModified Bessel Functions

§10.41 Asymptotic Expansions for Large Order

Contents

§10.41(ii) Uniform Expansions for Real Variable

As \nu\to\infty through positive real values,

uniformly for 0<z<\infty. Here

10.41.7\eta=(1+z^{2})^{{\frac{1}{2}}}+\mathop{\ln\/}\nolimits\frac{z}{1+(1+z^{2})^{{%
\frac{1}{2}}}},
10.41.8p=(1+z^{2})^{{-\frac{1}{2}}},

where the branches assume their principal values. Also, U_{k}(p) and V_{k}(p) are polynomials in p of degree 3k, given by U_{0}(p)=V_{0}(p)=1, and

10.41.9
U_{{k+1}}(p)=\tfrac{1}{2}p^{2}(1-p^{2})U_{k}^{{\prime}}(p)+\frac{1}{8}\int_{0}%
^{p}(1-5t^{2})U_{k}(t)dt,
V_{{k+1}}(p)=U_{{k+1}}(p)-\tfrac{1}{2}p(1-p^{2})U_{k}(p)-p^{2}(1-p^{2})U_{k}^{%
{\prime}}(p),k=0,1,2,\ldots.

For k=1,2,3,

10.41.10
U_{1}(p)=\tfrac{1}{24}(3p-5p^{3}),
U_{2}(p)=\tfrac{1}{1152}(81p^{2}-462p^{4}+385p^{6}),
U_{3}(p)=\tfrac{1}{4\;14720}\*(30375p^{3}-3\;69603p^{5}+7\;65765p^{7}-4\;25425%
p^{9}),
10.41.11
V_{1}(p)=\tfrac{1}{24}(-9p+7p^{3}),
V_{2}(p)=\tfrac{1}{1152}(-135p^{2}+594p^{4}-455p^{6}),
V_{3}(p)=\tfrac{1}{4\;14720}\*(-42525p^{3}+4\;51737p^{5}-8\;83575p^{7}+4\;7547%
5p^{9}).

For U_{4}(p), U_{5}(p), U_{6}(p), see Bickley et al. (1952, p. xxxv).

For numerical tables of \eta=\eta(z) and the coefficients U_{k}(p), V_{k}(p), see Olver (1962, pp. 43–51).

§10.41(iii) Uniform Expansions for Complex Variable

The expansions (10.41.3)–(10.41.6) also hold uniformly in the sector |\mathop{\mathrm{ph}\/}\nolimits z|\leq\tfrac{1}{2}\pi-\delta (<\tfrac{1}{2}\pi), with the branches of the fractional powers in (10.41.3)–(10.41.8) extended by continuity from the positive real z-axis.

Figures 10.41.1 and 10.41.2 show corresponding points of the mapping of the z-plane and the \eta-plane. The curve E_{1}BE_{2} in the z-plane is the upper boundary of the domain \mathbf{K} depicted in Figure 10.20.3 and rotated through an angle -\tfrac{1}{2}\pi. Thus B is the point z=c, where c is given by (10.20.18).

For derivations of the results in this subsection, and also error bounds, see Olver (1997b, pp. 374–378). For extensions of the regions of validity in the z-plane and extensions to complex values of \nu see Olver (1997b, pp. 378–382).

See accompanying text
Figure 10.41.1: z-plane. Magnify
See accompanying text
Figure 10.41.2: \eta-plane. Magnify

For expansions in inverse factorial series see Dunster et al. (1993).

§10.41(iv) Double Asymptotic Properties

In the case of (10.41.13) with positive real values of z the result is a consequence of the error bounds given in Olver (1997b, pp. 377–378). Then by expanding the quantities \eta, (1+z^{2})^{{-\frac{1}{4}}}, and U_{k}(p), k=0,1,\ldots,\ell-1, and rearranging, we arrive at an expansion of the right-hand side of (10.41.13) in powers of z^{{-1}}. Moreover, because of the uniqueness property of asymptotic expansions (§2.1(iii)) this expansion must agree with (10.40.2), with z replaced by \nu z, up to and including the term in z^{{-(\ell-1)}}. It also enjoys the same sector of validity.

To establish (10.41.12) we substitute into (10.34.3), with m=0 and z replaced by \nu z, by means of (10.41.13) observing that when |z| is large the effect of replacing z by ze^{{\pm\pi i}} is to replace \eta, (1+z^{2})^{{\frac{1}{4}}}, and p by -\eta, \pm i(1+z^{2})^{{\frac{1}{4}}}, and -p, respectively.

§10.41(v) Double Asymptotic Properties (Continued)

Similar analysis can be developed for the uniform asymptotic expansions in terms of Airy functions given in §10.20. We first prove that for the expansions (10.20.6) for the Hankel functions \mathop{{H^{{(1)}}_{{\nu}}}\/}\nolimits\!\left(\nu z\right) and \mathop{{H^{{(2)}}_{{\nu}}}\/}\nolimits\!\left(\nu z\right) the z-asymptotic property applies when z\to\pm i\infty, respectively. This is a consequence of the error bounds associated with these expansions. We then extend the validity of this property from z\to\pm i\infty to z\to\infty in the sector -\pi+\delta\leq\mathop{\mathrm{ph}\/}\nolimits z\leq 2\pi-\delta in the case of \mathop{{H^{{(1)}}_{{\nu}}}\/}\nolimits\!\left(\nu z\right), and to z\to\infty in the sector -2\pi+\delta\leq\mathop{\mathrm{ph}\/}\nolimits z\leq\pi-\delta in the case of \mathop{{H^{{(2)}}_{{\nu}}}\/}\nolimits\!\left(\nu z\right). This is done by re-expansion with the aid of (10.20.10), (10.20.11), and §10.41(ii), followed by comparison with (10.17.5) and (10.17.6), with z replaced by \nu z. Lastly, we substitute into (10.4.4), again with z replaced by \nu z. The final results are:

as z\to\infty in |\mathop{\mathrm{ph}\/}\nolimits z|\leq\pi-\delta, or equivalently as \zeta\to\infty in |\mathop{\mathrm{ph}\/}\nolimits(-\zeta)|\leq\tfrac{2}{3}\pi-\delta, for fixed \ell (\geq 0) and fixed \nu (>0).

It needs to be noted that the results (10.41.14) and (10.41.15) do not apply when z\to 0+ or equivalently \zeta\to+\infty. This is because A_{k}(\zeta) and \zeta^{{-\frac{1}{2}}}B_{k}(\zeta),k=0,1,\ldots, do not form an asymptotic scale (§2.1(v)) as \zeta\to+\infty; see Olver (1997b, pp. 422–425).