Digital Library of Mathematical Functions
About the Project
NIST
31 Heun FunctionsProperties

§31.8 Solutions via Quadratures

For half-odd-integer values of the exponent parameters:

the Hermite–Darboux method (see Whittaker and Watson (1927, pp. 570–572)) can be applied to construct solutions of (31.2.1) expressed in quadratures, as follows.

Denote \mathbf{m}=(m_{0},m_{1},m_{2},m_{3}) and \lambda=-4q. Then

are two independent solutions of (31.2.1). Here \Psi_{{g,N}}(\lambda,z) is a polynomial of degree g in \lambda and of degree N=m_{0}+m_{1}+m_{2}+m_{3} in z, that is a solution of the third-order differential equation satisfied by a product of any two solutions of Heun’s equation. The degree g is given by

The variables \lambda and \nu are two coordinates of the associated hyperelliptic (spectral) curve \Gamma:\nu^{2}=\prod_{{j=1}}^{{2g+1}}(\lambda-\lambda_{j}). (This \nu is unrelated to the \nu in §31.6.) Lastly, \lambda_{j}, j=1,2,\ldots,2g+1, are the zeros of the Wronskian of w_{+}(\mathbf{m};\lambda;z) and w_{-}(\mathbf{m};\lambda;z).

By automorphisms from §31.2(v), similar solutions also exist for m_{0},m_{1},m_{2},m_{3}\in\Integer, and \Psi_{{g,N}}(\lambda,z) may become a rational function in z. For instance,

and

For \mathbf{m}=(m_{0},0,0,0), these solutions reduce to Hermite’s solutions (Whittaker and Watson (1927, §23.7)) of the Lamé equation in its algebraic form. The curve \Gamma reflects the finite-gap property of Equation (31.2.1) when the exponent parameters satisfy (31.8.1) for m_{j}\in\Integer. When \lambda=-4q approaches the ends of the gaps, the solution (31.8.2) becomes the corresponding Heun polynomial. For more details see Smirnov (2002).

The solutions in this section are finite-term Liouvillean solutions which can be constructed via Kovacic’s algorithm; see §31.14(ii).