About the Project
28 Mathieu Functions and Hill’s EquationHill’s Equation

§28.31 Equations of Whittaker–Hill and Ince

Contents
  1. §28.31(i) Whittaker–Hill Equation
  2. §28.31(ii) Equation of Ince; Ince Polynomials
  3. §28.31(iii) Paraboloidal Wave Functions

§28.31(i) Whittaker–Hill Equation

Hill’s equation with three terms

28.31.1 W^{\prime\prime}+\left(A+B\cos\left(2z\right)-\tfrac{1}{2}(kc)^{2}\cos\left(4z%
\right)\right)W=0

and constant values of A,B,k, and c, is called the Equation of Whittaker–Hill. It has been discussed in detail by Arscott (1967) for k^{2}<0, and by Urwin and Arscott (1970) for k^{2}>0.

§28.31(ii) Equation of Ince; Ince Polynomials

When k^{2}<0, we substitute

in (28.31.1). The result is the Equation of Ince:

Formal 2\pi-periodic solutions can be constructed as Fourier series; compare §28.4:

28.31.4 w_{\mathit{e},s}(z)=\sum_{\ell=0}^{\infty}A_{2\ell+s}\cos(2\ell+s)z,s=0,1,
28.31.5 w_{\mathit{o},s}(z)=\sum_{\ell=0}^{\infty}B_{2\ell+s}\sin(2\ell+s)z,s=1,2,

where the coefficients satisfy

28.31.6
-2\eta A_{0}+(2+p)\xi A_{2}=0,
p\xi A_{0}+(4-\eta)A_{2}+\left(\tfrac{1}{2}p+2\right)\xi A_{4}=0,
(\tfrac{1}{2}p-\ell+1)\xi A_{2\ell-2}+\left(4\ell^{2}-\eta\right)A_{2\ell}+(%
\tfrac{1}{2}p+\ell+1)\xi A_{2\ell+2}=0,\ell\geq 2,
28.31.7
\left(1-\eta+\left(\tfrac{1}{2}p+\tfrac{1}{2}\right)\xi\right)A_{1}+\left(%
\tfrac{1}{2}p+\tfrac{3}{2}\right)\xi A_{3}=0,
(\tfrac{1}{2}p-\ell+\tfrac{1}{2})\xi A_{2\ell-1}+\left((2\ell+1)^{2}-\eta%
\right)A_{2\ell+1}+(\tfrac{1}{2}p+\ell+\tfrac{3}{2})\xi A_{2\ell+3}=0,\ell\geq 1,
28.31.8
\left(1-\eta-\left(\tfrac{1}{2}p+\tfrac{1}{2}\right)\xi\right)B_{1}+\left(%
\tfrac{1}{2}p+\tfrac{3}{2}\right)\xi B_{3}=0,
(\tfrac{1}{2}p-\ell+\tfrac{1}{2})\xi B_{2\ell-1}+\left((2\ell+1)^{2}-\eta%
\right)B_{2\ell+1}+(\tfrac{1}{2}p+\ell+\tfrac{3}{2})\xi B_{2\ell+3}=0,\ell\geq 1,
28.31.9
(4-\eta)B_{2}+\left(\tfrac{1}{2}p+2\right)\xi B_{4}=0,
(\tfrac{1}{2}p-\ell+1)\xi B_{2\ell-2}+(4\ell^{2}-\eta)B_{2\ell}+(\tfrac{1}{2}p%
+\ell+1)\xi B_{2\ell+2}=0,\ell\geq 2.

When p is a nonnegative integer, the parameter \eta can be chosen so that solutions of (28.31.3) are trigonometric polynomials, called Ince polynomials. They are denoted by

28.31.10 \begin{array}[]{cl}C_{2n}^{2m}(z,\xi)&\mbox{with $p=2n$},\\
C_{2n+1}^{2m+1}(z,\xi)&\mbox{with $p=2n+1$},\end{array}
28.31.11 \begin{array}[]{cl}S_{2n+1}^{2m+1}(z,\xi)&\mbox{with $p=2n+1$},\\
S_{2n+2}^{2m+2}(z,\xi)&\mbox{with $p=2n+2$},\end{array}

and m=0,1,\dots,n in all cases.

The values of \eta corresponding to C_{p}^{m}(z,\xi), S_{p}^{m}(z,\xi) are denoted by a_{p}^{m}(\xi), b_{p}^{m}(\xi), respectively. They are real and distinct, and can be ordered so that C_{p}^{m}(z,\xi) and S_{p}^{m}(z,\xi) have precisely m zeros, all simple, in 0\leq z<\pi. The normalization is given by

ambiguities in sign being resolved by requiring C_{p}^{m}(x,\xi) and {S_{p}^{m}}^{\prime}(x,\xi) to be continuous functions of x and positive when x=0.

For \xi\to 0, with x fixed,

If p\to\infty and \xi\to 0 in such a way that p\xi\to 2q, then in the notation of §§28.2(v) and 28.2(vi)

28.31.15
a_{p}^{m}(\xi)\to a_{m}(q),
b_{p}^{m}(\xi)\to b_{m}(q).

For proofs and further information, including convergence of the series (28.31.4), (28.31.5), see Arscott (1967).

§28.31(iii) Paraboloidal Wave Functions

With (28.31.10) and (28.31.11),

28.31.16 \mathit{hc}_{p}^{m}(z,\xi)=e^{-\frac{1}{4}\xi\cos\left(2z\right)}C_{p}^{m}(z,%
\xi),
28.31.17 \mathit{hs}_{p}^{m}(z,\xi)=e^{-\frac{1}{4}\xi\cos\left(2z\right)}S_{p}^{m}(z,%
\xi),

are called paraboloidal wave functions. They satisfy the differential equation

with \eta=a_{p}^{m}(\xi), \eta=b_{p}^{m}(\xi), respectively.

For proofs and further integral equations see Urwin (1964, 1965).

Asymptotic Behavior

For \xi>0, the functions \mathit{hc}_{p}^{m}(z,\xi), \mathit{hs}_{p}^{m}(z,\xi) behave asymptotically as multiples of \exp\left(-\tfrac{1}{4}\xi\cos\left(2z\right)\right)\left(\cos z\right)^{p} as z\to\pm\mathrm{i}\infty. All other periodic solutions behave as multiples of \exp\left(\tfrac{1}{4}\xi\cos\left(2z\right)\right)(\cos z)^{-p-2}.

For \xi>0, the functions \mathit{hc}_{p}^{m}(z,-\xi), \mathit{hs}_{p}^{m}(z,-\xi) behave asymptotically as multiples of \exp\left(\tfrac{1}{4}\xi\cos\left(2z\right)\right)(\cos z)^{-p-2} as z\to\tfrac{1}{2}\pi\pm\mathrm{i}\infty. All other periodic solutions behave as multiples of \exp\left(-\tfrac{1}{4}\xi\cos\left(2z\right)\right)\left(\cos z\right)^{p}.