Skip to content

Conventions and parameter maps for Heun functions

The fastest way to make a wrong statement about a Heun function is to write only its name and not its equation.

For Heun functions, the symbol alone is rarely enough. In practice you must also know

  1. the differential equation,
  2. the parameter order,
  3. the normalization or base point, and
  4. the analytic-continuation convention (for example a branch cut, or a local power series around an ordinary point).

This page is the site-wide reference for those choices on heun.xyz. It is not claiming that one universal convention exists. It does the opposite: it makes our conventions explicit and gives the translations you need when moving among

  • the attached reference volume Heun’s Differential Equations (Ronveaux, ed., 1995),
  • the NIST Digital Library of Mathematical Functions (DLMF),
  • the Wolfram Language documentation, and
  • Maple’s built-in Heun functions.

FunctionSingularities in a standard formLocal object used as the functionBase point on heun.xyzConvergence / global issue
HeunGfour regular singularities 0,1,a,0,1,a,\inftyFrobenius solution analytic at 00z=0z=0power series radius $\min(1,
HeunCregular at 0,10,1; irregular at \infty of rank 11Frobenius solution analytic at 00z=0z=0power series radius 11; principal branch cut convention matters
HeunDtwo irregular singularities (standardly at 0,0,\infty or, in Maple, at ±1\pm1)local analytic solution fixed at an ordinary pointz=1z=1 in the site conventionno Frobenius normalization at an irregular singular point
HeunBregular at 00; irregular at \infty of rank 22Frobenius solution analytic at 00z=0z=0entire in zz for the usual local solution
HeunTonly an irregular singularity at \infty of rank 33power-series solution fixed by two initial conditionsz=0z=0entire in zz; Stokes behavior appears only at infinity

We use the standard general Heun equation

y(z)+(γz+δz1+ϵza)y(z)+αβzqz(z1)(za)y(z)=0,y''(z)+\left(\frac{\gamma}{z}+\frac{\delta}{z-1}+\frac{\epsilon}{z-a}\right)y'(z) +\frac{\alpha\beta z-q}{z(z-1)(z-a)}\,y(z)=0,

with

ϵ=α+βγδ+1,a0,1.\epsilon=\alpha+\beta-\gamma-\delta+1, \qquad a\neq 0,1.

Our symbol

HeunG(a,q,α,β,γ,δ;z)\mathrm{HeunG}(a,q,\alpha,\beta,\gamma,\delta;z)

means the solution analytic at z=0z=0 and normalized by

HeunG(a,q,α,β,γ,δ;0)=1.\mathrm{HeunG}(a,q,\alpha,\beta,\gamma,\delta;0)=1.

This is the easiest family conventionally because the major CAS agree on it: Maple and the Wolfram Language use the same parameter order and the same local normalization.

Practical reading rule. If a paper writes Hl(a,q;α,β,γ,δ;z)Hl(a,q;\alpha,\beta,\gamma,\delta;z), that is usually the same local object as our HeunG, up to notation.


We use the DLMF / Wolfram confluent Heun equation

y(z)+(γz+δz1+ϵ)y(z)+αzqz(z1)y(z)=0.y''(z)+\left(\frac{\gamma}{z}+\frac{\delta}{z-1}+\epsilon\right)y'(z) +\frac{\alpha z-q}{z(z-1)}\,y(z)=0.

Our symbol

HeunC(q,α,γ,δ,ϵ;z)\mathrm{HeunC}(q,\alpha,\gamma,\delta,\epsilon;z)

means the Frobenius solution at z=0z=0 normalized by

HeunC(q,α,γ,δ,ϵ;0)=1.\mathrm{HeunC}(q,\alpha,\gamma,\delta,\epsilon;0)=1.

The local exponents are 00 and 1γ1-\gamma at z=0z=0, and 00 and 1δ1-\delta at z=1z=1. The Taylor series at z=0z=0 has radius of convergence 11. In the principal Wolfram convention the branch cut runs from z=1z=1 to \infty.

Maple uses

HeunC(αM,βM,γM,δM,ηM,z)\mathrm{HeunC}(\alpha_M,\beta_M,\gamma_M,\delta_M,\eta_M,z)

for the local solution of

U(z)+(αM+βM+1z+γM+1z1)U(z)+(μz+νz1)U(z)=0,U''(z)+\left(\alpha_M+\frac{\beta_M+1}{z}+\frac{\gamma_M+1}{z-1}\right)U'(z) +\left(\frac{\mu}{z}+\frac{\nu}{z-1}\right)U(z)=0,

where

μ=12(αMβMγM+αMβMβMγM)ηM,\mu=\frac12\bigl(\alpha_M-\beta_M-\gamma_M+\alpha_M\beta_M-\beta_M\gamma_M\bigr)-\eta_M, ν=12(αM+βM+γM+αMγM+βMγM)+δM+ηM.\nu=\frac12\bigl(\alpha_M+\beta_M+\gamma_M+\alpha_M\gamma_M+\beta_M\gamma_M\bigr)+\delta_M+\eta_M.

Then the equivalent site parameters are

ϵ=αM,γ=βM+1,δ=γM+1,\epsilon=\alpha_M, \qquad \gamma=\beta_M+1, \qquad \delta=\gamma_M+1, q=μ,α=μ+ν.q=\mu, \qquad \alpha=\mu+\nu.

Conversely, given the site parameters (q,α,γ,δ,ϵ)(q,\alpha,\gamma,\delta,\epsilon),

αM=ϵ,βM=γ1,γM=δ1,\alpha_M=\epsilon, \qquad \beta_M=\gamma-1, \qquad \gamma_M=\delta-1, δM=αϵ(γ+δ)2,ηM=12+γ(ϵδ)2q.\delta_M=\alpha-\frac{\epsilon(\gamma+\delta)}{2}, \qquad \eta_M=\frac12+\frac{\gamma(\epsilon-\delta)}{2}-q.

We use the Wolfram/DLMF-style double-confluent equation

z2y(z)+(γ+δz+ϵz2)y(z)+(αzq)y(z)=0,z^2y''(z)+\bigl(\gamma+\delta z+\epsilon z^2\bigr)y'(z)+(\alpha z-q)y(z)=0,

or equivalently

y(z)+(γz2+δz+ϵ)y(z)+αzqz2y(z)=0.y''(z)+\left(\frac{\gamma}{z^2}+\frac{\delta}{z}+\epsilon\right)y'(z) +\frac{\alpha z-q}{z^2}\,y(z)=0.

Our symbol

HeunD(q,α,γ,δ,ϵ;z)\mathrm{HeunD}(q,\alpha,\gamma,\delta,\epsilon;z)

means the local analytic solution fixed at the ordinary point z=1z=1 by

y(1)=1,y(1)=0.y(1)=1, \qquad y'(1)=0.

This choice is not cosmetic. In this convention z=0z=0 is an irregular singular point, so there is no Frobenius normalization there.

DLMF writes the doubly confluent equation in the reduced form

y(z)+(δDz2+γDz+1)y(z)+αDzqDz2y(z)=0.y''(z)+\left(\frac{\delta_D}{z^2}+\frac{\gamma_D}{z}+1\right)y'(z) +\frac{\alpha_D z-q_D}{z^2}\,y(z)=0.

If ϵ0\epsilon\neq 0, a rescaling of the independent variable reduces the site convention to the DLMF normalization with the constant term in yy' fixed to 11. After that rescaling, the two conventions differ only by a relabeling of the coefficients of 1/z21/z^2 and 1/z1/z.

Maple also has a built-in HeunD, but it uses a different normal form: the irregular singularities are placed at x=±1x=\pm1, the origin is an ordinary point, and the built-in function is expanded around x=0x=0. So Maple HeunD is not a simple parameter relabeling of the site convention.

Translation rule. For HeunD, never compare software output by symbol name alone. Compare by ODE and normalization.


We use the 4-parameter Ronveaux / Maple biconfluent form

zy(z)+(1+αβz2z2)y(z)+((γα2)z12[δ+(1+α)β])y(z)=0.z\,y''(z)+(1+\alpha-\beta z-2z^2)y'(z) +\Bigl((\gamma-\alpha-2)z-\tfrac12\bigl[\delta+(1+\alpha)\beta\bigr]\Bigr)y(z)=0.

Our symbol

HeunB(α,β,γ,δ;z)\mathrm{HeunB}(\alpha,\beta,\gamma,\delta;z)

means the solution analytic at z=0z=0 and normalized by

HeunB(α,β,γ,δ;0)=1.\mathrm{HeunB}(\alpha,\beta,\gamma,\delta;0)=1.

This convention is extremely convenient because

  • the exponents at z=0z=0 are transparent,
  • the local power series is straightforward,
  • polynomial truncation is easy to state,
  • and Maple uses the same 4-parameter structure.

Exact map from the site convention to Wolfram HeunB

Section titled “Exact map from the site convention to Wolfram HeunB”

Wolfram uses the 5-parameter equation

zy(z)+(γW+δWz+ϵWz2)y(z)+(αWzqW)y(z)=0,z\,y''(z)+\bigl(\gamma_W+\delta_W z+\epsilon_W z^2\bigr)y'(z)+(\alpha_W z-q_W)y(z)=0,

with the regular solution normalized by y(0)=1y(0)=1.

The exact parameter map from the site convention to Wolfram’s parameters is

ϵW=2,γW=1+α,δW=β,\epsilon_W=-2, \qquad \gamma_W=1+\alpha, \qquad \delta_W=-\beta, αW=γα2,qW=12(δ+(1+α)β).\alpha_W=\gamma-\alpha-2, \qquad q_W=\frac12\bigl(\delta+(1+\alpha)\beta\bigr).

DLMF uses a different reduced biconfluent form,

y(z)(γDz+δD+z)y(z)+αDzqDzy(z)=0,y''(z)-\left(\frac{\gamma_D}{z}+\delta_D+z\right)y'(z) +\frac{\alpha_D z-q_D}{z}\,y(z)=0,

which is equivalent to the site convention only after a rescaling of the independent variable and a linear redefinition of parameters.

So for HeunB the correct question is not “which symbol is right?”, but “which ODE am I actually solving?”


We use a reduced 3-parameter triconfluent form,

y(z)(γ+3z2)y(z)+[α+(β3)z]y(z)=0.y''(z)-\bigl(\gamma+3z^2\bigr)y'(z)+\bigl[\alpha+(\beta-3)z\bigr]y(z)=0.

Our symbol

HeunT(α,β,γ;z)\mathrm{HeunT}(\alpha,\beta,\gamma;z)

means the entire solution fixed by

y(0)=1,y(0)=0.y(0)=1, \qquad y'(0)=0.

This convention is well adapted to entire-function methods, Taylor recurrences, and polynomial solutions.

Wolfram uses the 5-parameter triconfluent equation

y(z)+(γW+δWz+ϵWz2)y(z)+(αWzqW)y(z)=0,y''(z)+\bigl(\gamma_W+\delta_W z+\epsilon_W z^2\bigr)y'(z) +(\alpha_W z-q_W)y(z)=0,

with the power-series normalization

y(0)=1,y(0)=0.y(0)=1, \qquad y'(0)=0.

If you want the site convention with the same independent variable and no further gauge transformation, then it is the specialization

γW=γ,δW=0,ϵW=3,\gamma_W=-\gamma, \qquad \delta_W=0, \qquad \epsilon_W=-3, αW=β3,qW=α.\alpha_W=\beta-3, \qquad q_W=-\alpha.

DLMF uses another 3-parameter triconfluent standard form,

y(z)+(γD+z)zy(z)+(αDzqD)y(z)=0,y''(z)+(\gamma_D+z)z\,y'(z)+(\alpha_D z-q_D)y(z)=0,

which belongs to the same triconfluent family but is not the same literal ODE as the site convention. In the triconfluent literature, different normal forms are common, and one often passes among them by normal-form reductions, rescalings, and gauge transformations.

So for HeunT, more than for any other member of the family, you should compare by the ODE itself and not by the function name alone.


A compact crosswalk of the five site conventions

Section titled “A compact crosswalk of the five site conventions”
Site pageHouse symbol on heun.xyzClosest direct CAS matchWhat to watch for
GeneralHeunG(a,q,\alpha,\beta,\gamma,\delta;z)Maple and Wolfram both match directlyremember the Fuchs relation for ϵ\epsilon
ConfluentHeunC(q,\alpha,\gamma,\delta,\epsilon;z)Wolfram matches directlyMaple uses (αM,βM,γM,δM,ηM)(\alpha_M,\beta_M,\gamma_M,\delta_M,\eta_M) and a different ODE presentation
Double confluentHeunD(q,\alpha,\gamma,\delta,\epsilon;z)Wolfram matches directlyMaple HeunD uses a different canonical variable placement and a different base point
BiconfluentHeunB(\alpha,\beta,\gamma,\delta;z)Maple is the closest direct matchWolfram uses a 5-parameter form; DLMF uses a reduced form
TriconfluentHeunT(\alpha,\beta,\gamma;z)reduced 3-parameter literature / Maple-style workflowsWolfram uses a 5-parameter form; DLMF uses another reduced normal form

Which parameter is the accessory parameter?

Section titled “Which parameter is the accessory parameter?”

This matters in spectral problems.

  • For HeunG, the accessory parameter is qq.
  • For HeunC, the accessory parameter is again qq in the DLMF / Wolfram convention; Maple packages the same freedom into ηM\eta_M together with the other coefficients.
  • For HeunD, the accessory parameter is naturally qq in the Wolfram / DLMF-like convention.
  • For HeunB, in the 4-parameter Ronveaux / Maple form the role most often played by an accessory parameter is carried by δ\delta.
  • For HeunT, different normal forms distribute the free coefficients differently; in the reduced 3-parameter form used on this site there is no universally preferred “one true accessory parameter” in the same way as for HeunG or HeunC.

Practical rule. The accessory parameter is the one not fixed by the local exponent data and usually determined by a global boundary condition.


Polynomial conditions: the site-wide rule of thumb

Section titled “Polynomial conditions: the site-wide rule of thumb”

The detailed truncation formulas depend strongly on the chosen convention, but the overall pattern is stable.

  • HeunG: one exponent parameter must be a nonpositive integer, and the accessory parameter must satisfy a finite algebraic condition.
  • HeunC: a linear truncation condition plus a determinant or accessory-root condition.
  • HeunB: a linear truncation condition plus a root condition in the accessory-like parameter.
  • HeunD: the usual double-confluent normalizations do not lead to the same simple Frobenius-polynomial story, because the canonical equation has no regular singular point at which a polynomial truncation is naturally defined.
  • HeunT: polynomial solutions occur when one coefficient parameter takes a positive-integer value and the remaining parameters satisfy a finite determinant condition.

So when you quote a “Heun polynomial condition,” always quote it in the same convention as the ODE you are using.


The common traps this site is designed to eliminate

Section titled “The common traps this site is designed to eliminate”

1. Confusing a local solution with a global eigenfunction

Section titled “1. Confusing a local solution with a global eigenfunction”

HeunG, HeunC, HeunB, and the ordinary-point normalizations of HeunD and HeunT are first and foremost local normalized solutions. In applications you often want something else: a solution satisfying conditions at multiple singular points or in multiple Stokes sectors. That usually turns the accessory parameter into an eigenvalue.

2. Comparing Maple and Wolfram values without converting parameters

Section titled “2. Comparing Maple and Wolfram values without converting parameters”

This is harmless for HeunG, but dangerous for HeunC, very dangerous for HeunD, and often fatal for HeunB and HeunT.

3. Forgetting that irregular singularities change what “normalization” means

Section titled “3. Forgetting that irregular singularities change what “normalization” means”

At a regular singular point you can define a Frobenius solution. At an irregular singular point you usually cannot. Then the function is defined instead by initial data at an ordinary point or by an asymptotic condition in a Stokes sector.

4. Treating branch issues and Stokes issues as if they were the same thing

Section titled “4. Treating branch issues and Stokes issues as if they were the same thing”
  • HeunG and HeunC have branch-cut conventions for principal values.
  • HeunD, HeunB, and HeunT are better thought of in terms of local analytic continuation and Stokes phenomena near irregular singular points.

5. Switching conventions in the middle of a derivation

Section titled “5. Switching conventions in the middle of a derivation”

Pick one convention, do the algebra in that convention, and only translate at the beginning or at the end.


What every function page on heun.xyz should state up front

Section titled “What every function page on heun.xyz should state up front”

Every page on this site should begin with a compact definition box containing the following items.

  1. Exact ODE. No shorthand, no hidden parameter constraints.
  2. Normalization. For example y(0)=1y(0)=1 or y(1)=1y(1)=1, y(1)=0y'(1)=0.
  3. Singularity structure. Regular vs irregular, and rank when irregular.
  4. Expansion point. The point around which the local series is defined.
  5. Closest CAS match. Maple, Wolfram, both, or neither directly.
  6. Parameter map note. A one-line warning if the same symbol means something else elsewhere.

That single box already removes most avoidable confusion.


  • A. Ronveaux (ed.), Heun’s Differential Equations, Oxford University Press, 1995. In particular Part B (confluent), Part C (double confluent), Part D (biconfluent), and Part E (triconfluent).
  • NIST Digital Library of Mathematical Functions, Chapter 31, especially §31.2 and §31.12: https://dlmf.nist.gov/31 and https://dlmf.nist.gov/31.12

Use this page first, then go to the function-specific guide:

  • HeunG for the general equation,
  • HeunC for the singly confluent case,
  • HeunD for the double-confluent case,
  • HeunB for the biconfluent case,
  • HeunT for the triconfluent case.

When those pages mention a “site convention,” they mean exactly the conventions fixed here.