9 Airy and Related FunctionsAiry Functions

§9.8 Modulus and Phase

Contents
  1. §9.8(i) Definitions
  2. §9.8(ii) Identities
  3. §9.8(iii) Monotonicity
  4. §9.8(iv) Asymptotic Expansions

§9.8(i) Definitions

Throughout this section x is real and nonpositive.

9.8.1 Ai⁡(x) =M⁡(x)⁢sin⁡θ⁡(x),
9.8.2 Bi⁡(x) =M⁡(x)⁢cos⁡θ⁡(x),
9.8.3 M⁡(x) =Ai2⁡(x)+Bi2⁡(x),
9.8.4 θ⁡(x) =arctan⁡(Ai⁡(x)/Bi⁡(x)).
9.8.5 Ai′⁡(x) =N⁡(x)⁢sin⁡ϕ⁡(x),
9.8.6 Bi′⁡(x) =N⁡(x)⁢cos⁡ϕ⁡(x),
9.8.7 N⁡(x) =Ai′2⁡(x)+Bi′2⁡(x),
9.8.8 ϕ⁡(x) =arctan⁡(Ai′⁡(x)/Bi′⁡(x)).

Graphs of M⁡(x) and N⁡(x) are included in §9.3(i). The branches of θ⁡(x) and ϕ⁡(x) are continuous and fixed by θ⁡(0)=−ϕ⁡(0)=16⁢π. (These definitions of θ⁡(x) and ϕ⁡(x) differ from Abramowitz and Stegun (1964, Chapter 10), and agree more closely with those used in Miller (1946) and Olver (1997b, Chapter 11).)

In terms of Bessel functions, and with ξ=23⁢|x|3/2,

9.8.9 |x|1/2⁢M2⁡(x) =12⁢ξ⁢(J1/32⁡(ξ)+Y1/32⁡(ξ)),
9.8.10 |x|−1/2⁢N2⁡(x) =12⁢ξ⁢(J2/32⁡(ξ)+Y2/32⁡(ξ)),

§9.8(ii) Identities

Primes denote differentiations with respect to x, which is continued to be assumed real and nonpositive.

9.8.13 M⁡(x)⁢N⁡(x)⁢sin⁡(θ⁡(x)−ϕ⁡(x))=π−1,
9.8.14 M2⁡(x)⁢θ′⁡(x) =−π−1,
N2⁡(x)⁢ϕ′⁡(x) =π−1⁢x,
N⁡(x)⁢N′⁡(x) =x⁢M⁡(x)⁢M′⁡(x),
9.8.15 N2⁡(x) =M′2⁡(x)+M2⁡(x)⁢θ′2⁡(x)=M′2⁢(x)+π−2⁢M−2⁡(x),
9.8.16 x2⁢M2⁡(x) =N′2⁡(x)+N2⁡(x)⁢ϕ′2⁡(x)=N′2⁡(x)+π−2⁢x2⁢N−2⁡(x),
9.8.17 tan⁡(θ⁡(x)−ϕ⁡(x))=1/(π⁢M⁡(x)⁢M′⁡(x))=−M⁡(x)⁢θ′⁡(x)/M′⁡(x),
9.8.18 M′′⁡(x)=x⁢M⁡(x)+π−2⁢M−3⁡(x),
M2′′′⁡(x)−4⁢x⁢M2′⁡(x)−2⁢M2⁡(x)=0,
9.8.19 θ′2⁡(x)+12⁢(θ′′′⁡(x)/θ′⁡(x))−34⁢(θ′′⁡(x)/θ′⁡(x))2=−x.

§9.8(iii) Monotonicity

As x increases from −∞ to 0 each of the functions M⁡(x), M′⁡(x), |x|−1/4⁢N⁡(x), M⁡(x)⁢N⁡(x), θ′⁡(x), ϕ′⁡(x) is increasing, and each of the functions |x|1/4⁢M⁡(x), θ⁡(x), ϕ⁡(x) is decreasing.

§9.8(iv) Asymptotic Expansions

As x→−∞

9.8.20 M2⁡(x) ∼1π⁢(−x)1/2⁢∑k=0∞1⋅3⋅5⁢⋯⁢(6⁢k−1)k!⁢(96)k⁢1x3⁢k,
9.8.21 N2⁡(x) ∼(−x)1/2π⁢∑k=0∞1⋅3⋅5⁢⋯⁢(6⁢k−1)k!⁢(96)k⁢1+6⁢k1−6⁢k⁢1x3⁢k,
9.8.22 θ⁡(x) ∼π4+23⁢(−x)3/2⁢(1+532⁢1x3+11056144⁢1x6+8282565536⁢1x9+12820 31525587 20256⁢1x12+⋯),
9.8.23 ϕ⁡(x) ∼−π4+23⁢(−x)3/2⁢(1−732⁢1x3−14636144⁢1x6−4 952713 27680⁢1x9−2065 3042983 88608⁢1x12−⋯).

The remainder after n terms does not exceed the (n+1)th term in absolute value and is of the same sign, provided that n≥0 for (9.8.20), (9.8.22) and (9.8.23), and n≥1 for (9.8.21).

For higher terms in (9.8.22) and (9.8.23) see Fabijonas et al. (2004). Also, approximate values (25S) of the coefficients of the powers x−15, x−18, …, x−56 are available in Sherry (1959).