HeunT - Maple Programming Help

HeunT

The Heun Triconfluent function

HeunTPrime

The derivative of the Heun Triconfluent function

 Calling Sequence HeunT($\mathrm{\alpha }$, $\mathrm{\beta }$, $\mathrm{\gamma }$, z) HeunTPrime($\mathrm{\alpha }$, $\mathrm{\beta }$, $\mathrm{\gamma }$, z)

Parameters

 $\mathrm{\alpha }$ - algebraic expression $\mathrm{\beta }$ - algebraic expression $\mathrm{\gamma }$ - algebraic expression z - algebraic expression

Description

 • The HeunT function is the solution of the Heun Triconfluent equation. Following the first reference (at the end), the equation and the conditions at the origin satisfied by HeunT are
 ${\mathrm{HeunT}}{}\left({\mathrm{α}}{,}{\mathrm{β}}{,}{\mathrm{γ}}{,}{z}\right){=}{\mathrm{DESol}}{}\left(\left\{\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{z}}^{{2}}}{}{\mathrm{_Y}}{}\left({z}\right){-}\left({3}{}{{z}}^{{2}}{+}{\mathrm{γ}}\right){}\left(\frac{{ⅆ}}{{ⅆ}{z}}{}{\mathrm{_Y}}{}\left({z}\right)\right){-}\left(\left({-}{\mathrm{β}}{+}{3}\right){}{z}{-}{\mathrm{α}}\right){}{\mathrm{_Y}}{}\left({z}\right)\right\}{,}\left\{{\mathrm{_Y}}{}\left({z}\right)\right\}{,}\left\{{\mathrm{_Y}}{}\left({0}\right){=}{1}{,}{\mathrm{D}}{}\left({\mathrm{_Y}}\right){}\left({0}\right){=}{0}\right\}\right)$ (1)
 • The HeunT($\mathrm{\alpha }$,$\mathrm{\beta }$,$\mathrm{\gamma }$,z) function is a local solution to Heun's Triconfluent equation, computed as a standard power series expansion around the origin, a regular point. Because the single singularity is located at $\mathrm{\infty }$, this series converges in the whole complex plane.
 • The Triconfluent Heun Equation (THE) above is obtained from the Doubleconfluent Heun Equation (DHE) through a confluence process, that is, a process where two singularities coalesce, performed by redefining parameters and taking limits. In this case the two irregular singularities of the DHE are coalesced into one irregular singularity at $\mathrm{\infty }$. The resulting Heun Triconfluent equation, thus, has the structure of singularities f of the 0F1 hypergeometric equation and so can be related to the Airy functions.
 • A special case happens when in HeunT($\mathrm{\alpha }$,$\mathrm{\beta }$,$\mathrm{\gamma }$,z), the second parameter satisfies $\mathrm{\beta }=3\left(1+n\right)$, where $n$ is a positive integer. In this case the $n$th+1, $n$th+2 and $n$th+3 coefficients form a polynomial system for the remaining parameters $\mathrm{\alpha }$ and $\mathrm{\gamma }$; when this system is identically satisfied all the subsequent coefficients cancel too and the series truncates, resulting in a polynomial form of degree $n$ for HeunT. Remark: for $n=0$ this situation leads to a constant, for $n=1$ HeunT will also be a constant since its series expansion satisfies $\mathrm{HeunT}\text{'}$ at 0 = 0 and for $n=2$ the polynomial system for $\mathrm{\alpha }$ and $\mathrm{\gamma }$ is inconsistent. So the non-trivial polynomial forms of HeunT are of degree $3\le n$.

Examples

Heun's Triconfluent equation,

 > $\mathrm{THE}≔\frac{{ⅆ}^{2}}{ⅆ{z}^{2}}y\left(z\right)=\left(3{z}^{2}+\mathrm{γ}\right)\left(\frac{ⅆ}{ⅆz}y\left(z\right)\right)+\left(\left(3-\mathrm{β}\right)z-\mathrm{α}\right)y\left(z\right)$
 ${\mathrm{THE}}{≔}\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{z}}^{{2}}}{}{y}{}\left({z}\right){=}\left({3}{}{{z}}^{{2}}{+}{\mathrm{γ}}\right){}\left(\frac{{ⅆ}}{{ⅆ}{z}}{}{y}{}\left({z}\right)\right){+}\left(\left({3}{-}{\mathrm{β}}\right){}{z}{-}{\mathrm{α}}\right){}{y}{}\left({z}\right)$ (2)

can be transformed into another version of itself, that is, an equation with one regular and one irregular singularities respectively located at 0 and $\mathrm{\infty }$ through transformations of the form

 > $z=\mathrm{κ}t,y\left(z\right)={ⅇ}^{\frac{1t\left(\mathrm{κ}-1\right)\left({\mathrm{κ}}^{2}+\mathrm{κ}+1\right)\left(\mathrm{γ}+{t}^{2}\right)}{2{\mathrm{κ}}^{3}}}u\left(t\right)$
 ${z}{=}{\mathrm{κ}}{}{t}{,}{y}{}\left({z}\right){=}{{ⅇ}}^{\frac{{1}}{{2}}{}\frac{{t}{}\left({\mathrm{κ}}{-}{1}\right){}\left({{\mathrm{κ}}}^{{2}}{+}{\mathrm{κ}}{+}{1}\right){}\left({{t}}^{{2}}{+}{\mathrm{γ}}\right)}{{{\mathrm{κ}}}^{{3}}}}{}{u}{}\left({t}\right)$ (3)

where $\left\{t,u\left(t\right)\right\}$ are new variables and ${\mathrm{\kappa }}^{6}=1$. Under this transformation, the HeunT parameters transform according to $\mathrm{\alpha }$ -> $\frac{\mathrm{\alpha }}{{\mathrm{\kappa }}^{2}}$, $\mathrm{\beta }$ -> $\frac{\mathrm{\beta }}{{\mathrm{\kappa }}^{3}}$, $\mathrm{\gamma }$ -> $\frac{\mathrm{\gamma }}{{\mathrm{\kappa }}^{4}}$. These transformations form a group of six elements and imply on identities, among which you have

 > $\mathrm{FunctionAdvisor}\left(\mathrm{identities},\mathrm{HeunT}\right)$
 $\left[\left[{\mathrm{HeunT}}{}\left({\mathrm{α}}{,}{\mathrm{β}}{,}{\mathrm{γ}}{,}{z}\right){=}{\mathrm{HeunT}}{}\left({j}{}{\mathrm{α}}{,}{\mathrm{β}}{,}{{j}}^{{2}}{}{\mathrm{γ}}{,}{j}{}{z}\right){,}{\mathrm{And}}{}\left({{j}}^{{3}}{=}{1}\right)\right]{,}\left[{\mathrm{HeunT}}{}\left({\mathrm{α}}{,}{\mathrm{β}}{,}{\mathrm{γ}}{,}{z}\right){=}{\mathrm{HeunT}}{}\left({\mathrm{α}}{,}{-}{\mathrm{β}}{,}{\mathrm{γ}}{,}{-}{z}\right){}{{ⅇ}}^{{{z}}^{{3}}}{,}{\mathrm{And}}{}\left({\mathrm{γ}}{=}{0}\right)\right]\right]$ (4)

When, in HeunT($\mathrm{\alpha }$,$\mathrm{\beta }$,$\mathrm{\gamma }$,z), $\mathrm{\beta }=3\left(1+n\right)$, where $n$ is a positive integer, the $n$th+1, $n$th+2 and $n$th+3 coefficients form a polynomial system for the remaining parameters $\mathrm{\alpha }$ and $\mathrm{\gamma }$. When this system is identically satisfied all the subsequent coefficients cancel too and the series truncates, resulting in a polynomial form of degree $n$ for HeunT. For example, this is the necessary condition for a polynomial form

 > $\mathrm{HeunT}\left(\mathrm{α},3n+3,\mathrm{γ},z\right)$
 ${\mathrm{HeunT}}{}\left({\mathrm{α}}{,}{3}{}{n}{+}{3}{,}{\mathrm{γ}}{,}{z}\right)$ (5)

Considering the first non-trivial case, for $n=3$, the function is

 > $\mathrm{HT}≔\mathrm{subs}\left(n=3,\right)$
 ${\mathrm{HT}}{≔}{\mathrm{HeunT}}{}\left({\mathrm{α}}{,}{12}{,}{\mathrm{γ}}{,}{z}\right)$ (6)

So the coefficients of ${z}^{m}$ for $m$ equal to 4, 5, and 6 in the series expansion are

 > $Q≔\mathrm{simplify}\left(\mathrm{series}\left(\mathrm{HT},z,7\right),\mathrm{size}\right)$
 ${Q}{≔}{1}{-}\frac{{1}}{{2}}{}{\mathrm{α}}{}{{z}}^{{2}}{+}\left({-}\frac{{1}}{{6}}{}{\mathrm{γ}}{}{\mathrm{α}}{-}\frac{{3}}{{2}}\right){}{{z}}^{{3}}{+}\left(\frac{{1}}{{24}}{}{{\mathrm{α}}}^{{2}}{-}\frac{{1}}{{24}}{}{{\mathrm{γ}}}^{{2}}{}{\mathrm{α}}{-}\frac{{3}}{{8}}{}{\mathrm{γ}}\right){}{{z}}^{{4}}{-}\frac{{1}}{{120}}{}\left({{\mathrm{γ}}}^{{2}}{-}{2}{}{\mathrm{α}}\right){}\left({\mathrm{α}}{}{\mathrm{γ}}{+}{9}\right){}{{z}}^{{5}}{+}\left({-}\frac{{1}}{{720}}{}{{\mathrm{α}}}^{{3}}{+}\frac{{1}}{{240}}{}{{\mathrm{γ}}}^{{2}}{}{{\mathrm{α}}}^{{2}}{+}\frac{{1}}{{720}}{}\left({-}{{\mathrm{γ}}}^{{4}}{+}{27}{}{\mathrm{γ}}\right){}{\mathrm{α}}{-}\frac{{1}}{{80}}{}{{\mathrm{γ}}}^{{3}}\right){}{{z}}^{{6}}{+}{\mathrm{O}}\left({{z}}^{{7}}\right)$ (7)
 > $\mathrm{c4},\mathrm{c5},\mathrm{c6}≔\mathrm{coeff}\left(Q,z,4\right),\mathrm{coeff}\left(Q,z,5\right),\mathrm{coeff}\left(Q,z,6\right)$
 ${\mathrm{c4}}{,}{\mathrm{c5}}{,}{\mathrm{c6}}{≔}\frac{{1}}{{24}}{}{{\mathrm{α}}}^{{2}}{-}\frac{{1}}{{24}}{}{{\mathrm{γ}}}^{{2}}{}{\mathrm{α}}{-}\frac{{3}}{{8}}{}{\mathrm{γ}}{,}{-}\frac{{1}}{{120}}{}\left({{\mathrm{γ}}}^{{2}}{-}{2}{}{\mathrm{α}}\right){}\left({\mathrm{α}}{}{\mathrm{γ}}{+}{9}\right){,}{-}\frac{{1}}{{720}}{}{{\mathrm{α}}}^{{3}}{+}\frac{{1}}{{240}}{}{{\mathrm{γ}}}^{{2}}{}{{\mathrm{α}}}^{{2}}{+}\frac{{1}}{{720}}{}\left({-}{{\mathrm{γ}}}^{{4}}{+}{27}{}{\mathrm{γ}}\right){}{\mathrm{α}}{-}\frac{{1}}{{80}}{}{{\mathrm{γ}}}^{{3}}$ (8)

solving for $\mathrm{\alpha }$ and $\mathrm{\gamma }$, requesting from solve to return using RootOf, you have

 > $\mathrm{_EnvExplicit}≔\mathrm{false}$
 ${\mathrm{_EnvExplicit}}{≔}{\mathrm{false}}$ (9)
 > $\mathrm{subs}\left(\mathrm{ga}=\mathrm{γ},\left[\mathrm{solve}\left(\mathrm{subs}\left(\mathrm{γ}=\mathrm{ga},\left\{\mathrm{c4},\mathrm{c5},\mathrm{c6}\right\}\right),\left\{\mathrm{α},\mathrm{ga}\right\}\right)\right]\right)$
 $\left[\left\{{\mathrm{α}}{=}{0}{,}{\mathrm{γ}}{=}{0}\right\}{,}\left\{{\mathrm{α}}{=}\frac{{1}}{{2}}{}{{\mathrm{RootOf}}{}\left({{\mathrm{_Z}}}^{{3}}{+}{36}\right)}^{{2}}{,}{\mathrm{γ}}{=}{\mathrm{RootOf}}{}\left({{\mathrm{_Z}}}^{{3}}{+}{36}\right)\right\}\right]$ (10)

substituting for instance the first of these two solutions in HT we have

 > $\mathrm{HT_polynomial}≔\mathrm{subs}\left({1}_{},\mathrm{HT}\right)$
 ${\mathrm{HT_polynomial}}{≔}{\mathrm{HeunT}}{}\left({0}{,}{12}{,}{0}{,}{z}\right)$ (11)

When the function admits a polynomial form, as is the case of HT_polynomial by construction, to obtain the actual polynomial of degree n (in this case n=3) use

 > $\genfrac{}{}{0}{}{\phantom{\mathrm{HeunT}=\mathrm{HeunT}:-\mathrm{SpecialValues}:-\mathrm{Polynomial}}}{}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}|\phantom{\rule[-0.0ex]{0.1em}{0.0ex}}\genfrac{}{}{0}{}{\phantom{}}{\mathrm{HeunT}=\mathrm{HeunT}:-\mathrm{SpecialValues}:-\mathrm{Polynomial}}$
 ${1}{-}\frac{{3}}{{2}}{}{{z}}^{{3}}$ (12)

References

 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.