The Heun Confluent function
The derivative of the Heun Confluent function
HeunC(α, β, γ, δ, η, z)
HeunCPrime(α, β, γ, δ, η, z)
The HeunC function is the solution of the Heun Confluent equation. Following the first reference (at the end), the equation and the conditions at the origin satisfied by HeunC are
This Heun (singly) Confluent equation is obtained from the Heun General equation through a confluence process, that is, a process where two singularities coalesce, performed by redefining parameters and taking limits, resulting in a single (typically irregular) singularity. The Heun Confluent equation thus has two regular singularities and one irregular singularity, and includes as particular cases both the 2F1 and 1F1 hypergeometric equations. The solution to the 2F1 equation,
DEtools[hyperode]( hypergeom([a,b],[c],z), y(z) ) = 0;
can then be expressed in terms of HeunC functions
and the same for the 1F1 hypergeometric confluent equation
DEtools[hyperode]( hypergeom([a],[c],z), y(z) ) = 0;
HeunC, thus, contains as particular cases all the hypergeometric functions of the 2F1 and 1F1 classes - some of these specializations are listed at the end of the Examples section.
Two other important non-hypergeometric case of Heun's Confluent equation, are the "spheroidal wave function" equation
diff(y(z),z,z) + 2*(gamma+1)*z*diff(y(z),z)/(z^2-1) + (4*delta*z^2-c)/(z^2-1)*y(z) = 0;
obtained from Heun's Confluent equation taking α=0,β=−12,η=1−γ−c4 and changing z -> z2;
and the rational form of Mathieu's equation,
diff(y(z),z,z) + z/(z^2-1)*diff(y(z),z) + (2*delta*(2*z^2-1)-a)/(z^2-1)*y(z) = 0;
obtained from the spheroidal wave function equation above by taking c=a+2⁢δ and further specializing γ=−12.
The HeunC(α,β,γ,δ,η, z) function is a local (Frobenius) solution to Heun's Confluent equation, computed as a power series expansion around the origin, a regular singular point. The series converges for z<1, where the second regular singularity is located. An analytic continuation of HeunC is obtained through identities, relating the values of the function in different regions of the z plane, for given values of the other parameters, or by expanding the solution around 1, the other regular singularity, and overlapping the series. General formulas relating these series expansions at different singularities and for arbitrary values of the other parameters, however, are not known at present.
A special case happens when the parameters entering HeunC are such that the function is, simultaneously, a Frobenius solution around the two regular singularities and hence analytic in a domain containing both of them. In such a case the series expansion for HeunC truncates and the function becomes a polynomial. A necessary (not sufficient) condition for this case is that δ=−n+γ+β+22⁢α, with n a positive integer, and η has one of a finite number of characteristic values, so that the function is a polynomial of degree n.
Heun's Confluent equation,
CHE ≔ ⅆ2ⅆz2⁢y⁡z=−z2⁢α+−2−β−gamma+α⁢z+1+β⁢ⅆⅆz⁢y⁡zz⁢z−1+1⁢−β−gamma−2⁢α−2⁢δ⁢z+β+1⁢α+−gamma−1⁢β−gamma−2⁢η⁢y⁡z2⁢z⁢z−1
can be transformed into another version of itself, that is, an equation with two regular singularities and one irregular singularity respectively located at 0,1,∞, through transformations of the form
where λ2=1, μ2=1 and ν2=1. Under this transformation, the HeunC parameters transform according to α -> α⁢λ, β -> β⁢μ and γ -> γ⁢ν. These transformations form a group of eight elements and imply on a number of identities, among which you have
Changing z -> 1−t also results in a HeunC equation with the singularities located at 0,1,∞; this permits rewriting the solution to the CHE in different manners. For example, the general solution returned by default by dsolve is
When β is an integer, however, these two "independent" solutions are not independent, yet a second pair of independent solutions can be constructed exploring this invariance in form under z -> 1−t
For z different from 1, the 2F1 and the confluent 1F1 hypergeometric functions are related to HeunC by
When δ=−n+γ+β+22⁢α, with n a positive integer, the nth + 1 coefficient in the series expansion is a polynomial in η of order n+1. If δ is a root of that polynomial, that coefficient is zero and with it all the following ones; the series then truncates and HeunC is a polynomial. For example, the necessary condition for a polynomial form is
Considering the first non-trivial case, for n=1, the function is
HC ≔ subs⁡n=1,
So the coefficient of z2 in the series expansion is
Q ≔ simplify⁡series⁡HC,z,3,size
c2 ≔ coeff⁡Q,z,2
solving for η, requesting from solve to return using RootOf, you have
_EnvExplicit ≔ false
substituting in HC we have
HC_polynomial ≔ subs⁡,HC
When the function admits a polynomial form, as is the case of HC_polynomial by construction, to obtain the actual polynomial of degree n (in this case n=1) use
Decarreau, A.; Dumont-Lepage, M.C.; Maroni, P.; Robert, A.; and Ronveaux, A. "Formes Canoniques de Equations confluentes de l'equation de Heun." Annales de la Societe Scientifique de Bruxelles. Vol. 92 I-II, (1978): 53-78.
Ronveaux, A. ed. Heun's Differential Equations. Oxford University Press, 1995.
Slavyanov, S.Y., and Lay, W. Special Functions, A Unified Theory Based on Singularities. Oxford Mathematical Monographs, 2000.
Download Help Document