HeunG - Maple Programming Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Conversions : Function : HeunG

HeunG

The Heun general function

HeunGPrime

The derivative of the Heun general function

 

Calling Sequence

Parameters

Description

Examples

References

Calling Sequence

HeunG(a, q, α, β, γ, δ, z)

HeunGPrime(a, q, α, β, γ, δ, z)

Parameters

a

-

algebraic expression

q

-

algebraic expression

α

-

algebraic expression

β

-

algebraic expression

γ

-

algebraic expression

δ

-

algebraic expression

z

-

algebraic expression

Description

• 

The HeunG function is the solution of the Heun General equation. Following the first reference (at the end), the equation and the conditions at the origin satisfied by HeunG are

FunctionAdvisor(definition, HeunG);

HeunGa,q,α,β,γ,δ,z=DESolⅆ2ⅆz2_Yzα+β+1z2+δγaα+δβ1z+γaⅆⅆz_Yzzz1z+aαβzq_Yzzz1z+a,_Yz,_Y0=1,D_Y0=qγa

(1)
  

Heun's equation is an extension of the 2F1 hypergeometric equation in that it is a second-order Fuchsian equation with four regular singular points. The 2F1 equation has three regular singularities. The HeunG function, thus, contains as particular cases all the functions of the hypergeometric 2F1 class.

• 

Another important case of Heun's equation, for q=ah4,α=ν2,β=12+ν2,δ=12,γ=12, is Lame's equation in algebraic form,

diff(y(z),z,z) + 1/2*(1/z + 1/(z-1) + 1/(z-a))*diff(y(z),z) + (a*h-nu*(nu+1)*z)/(4*z*(z-1)*(z-a))*y(z) = 0;

ⅆ2ⅆz2yz+121z+1z1+1zaⅆⅆzyz+14ahν1+νzyzzz1za=0

(2)
  

where the parameter ν is called the order of the equation and many special features arise when ν is an integer. Lame's equation arises in the separation of variables in Laplace's equation.

• 

The HeunG(a,q,α,β,γ,δ,z) function is a local (Frobenius) solution to Heun's equation, computed as a power series expansion around the origin, a regular singular point. The radius of convergence of this series is zs where zs is the location of the singularity closer to the origin, either a or 1. An analytic continuation of HeunG is obtained through identities, relating the values of the function in different regions, or by expanding the solution around the other singularities (a, 1 or ). For example, the value of the function outside a circle of radius 2 is computable from the value of another HeunG function inside that circle using

HeunG(a,q,alpha,beta,gamma,delta,z) = (1-z)^(-alpha)*HeunG(a/(a-1),(-q+gamma*alpha*a)/(a-1),alpha,alpha-delta+1,gamma,alpha-beta+1,z/(z-1));

HeunGa,q,α,β,γ,δ,z=1zαHeunGaa1,aαγqa1,α,αδ+1,γ,αβ+1,zz1

(3)
• 

For certain values of the parameters, it can happen that the expansion around the origin is also a Frobenius solution around the next adjacent singularity, so the function is analytic in some domain including both singularities. In the literature, the term Heun function is sometimes reserved for these cases, important in physical applications.

• 

A more special situation happens when the parameters entering HeunG are such that the function is, simultaneously, a Frobenius solution around three adjacent singularities and hence analytic in a domain containing all of them. In such a case the solution will also be a Frobenius solution around the fourth singularity and HeunG will be a polynomial. A necessary (not sufficient) condition for this case is that α=n, with n a positive integer, and q has one of a finite number of characteristic values, in which case the function is a polynomial of degree n.

Examples

Heun's equation can be transformed into another version of itself, that is, an equation with four regular singularities three of which are located at 0,1,, by any of 24 Mobius transformations z -> fz of the independent variable z; these forms of fz are

Matrixz,1z,1z,11z,zz1,z1z,za,aza,az,aaz,zza,zaz,za1a,z1a1,1aza,a1z1,zaz1,z1za,zaaz1,a1zza,az1za,az1a1z,za1az,1azza

z1z1z11zzz1z1zzaazaazaazzzazazza1az1a11azaa1z1zaz1z1zazaaz1a1zzaaz1zaaz1a1zza1az1azza

(4)

Note that the location of the fourth singularity resulting from these transformations, say λ, is in general different from a. The six possible values of λ are a,1a,11a,aa1,1a and a1a.

Taking into account that Heun's equation has 4 regular singularities, that at around each one it is possible to construct 2 Frobenius solutions, and that there exist these 24 transformations mapping the equation into one of the same type, the solution to Heun's equation can be written in 192 different manners. This situation is equivalent to the one of the 2F1 hypergeometric equation, where instead of 4 there are 3 regular singularities, instead of 24 transformations there are only 6, and so the solution can be written in 24 different manners.

This group of transformations of order 24 admitted by Heun's equation also leads to a rather large number of identities for the function solution HeunG, among which you have

FunctionAdvisoridentities,HeunG

HeunGa,q,α,β,γ,δ,z=1z1δHeunGa,qδ1γa,βδ+1,αδ+1,γ,2δ,z,HeunGa,q,α,β,γ,δ,z=1zaαβ+γ+δHeunGa,qγα+βγδ,β+γ+δ,α+γ+δ,γ,δ,z,Anda0,HeunGa,q,α,β,γ,δ,z=1z1δ1zaαβ+γ+δHeunGa,qγδ1a+α+βγδ,β+γ+1,α+γ+1,γ,2δ,z,Anda0,HeunGa,q,α,β,γ,δ,z=HeunG1a,qa,α,β,γ,α+βγδ+1,za,Anda0,HeunGa,q,α,β,γ,δ,z=1zαHeunGaa1,aαγqa1,α,αδ+1,γ,αβ+1,zz1,Anda1,z1,HeunGa,q,α,β,γ,δ,z=1zaαHeunG1a,αγq,α,β+γ+δ,γ,δ,1azza,Anda0,za,HeunGa,q,α,β,γ,δ,z=1zaαHeunG11a,αγ+qa1,α,β+γ+δ,γ,αβ+1,zza,Anda0,a1,za

(5)

Some hypergeometric special cases of HeunG are

FunctionAdvisorspecialize,HeunG,hypergeom

HeunGa&comma;q&comma;&alpha;&comma;&beta;&comma;&gamma;&comma;&delta;&comma;z&equals;hypergeom&alpha;&comma;&beta;&comma;&alpha;&delta;&plus;1&plus;&beta;&comma;z&comma;Anda&equals;0&comma;q&equals;0&comma;HeunGa&comma;q&comma;&alpha;&comma;&beta;&comma;&gamma;&comma;&delta;&comma;z&equals;hypergeom&alpha;&comma;&beta;&comma;&gamma;&comma;z&comma;Anda&equals;1&comma;q&equals;&alpha;&beta;orAndq&equals;a&alpha;&beta;&comma;&delta;&equals;&alpha;&plus;&beta;&gamma;&plus;1&comma;HeunGa&comma;q&comma;&alpha;&comma;&beta;&comma;&gamma;&comma;&delta;&comma;z&equals;hypergeom12&alpha;&beta;22&alpha;&plus;&beta;&gamma;&plus;&gamma;2&plus;4q&plus;12&gamma;&plus;12&alpha;12&beta;&comma;12&alpha;&beta;22&alpha;&plus;&beta;&gamma;&plus;&gamma;2&plus;4q&plus;12&gamma;12&alpha;&plus;12&beta;&comma;&gamma;&comma;zz112&alpha;&beta;22&alpha;&plus;&beta;&gamma;&plus;&gamma;2&plus;4q12&gamma;&plus;12&alpha;&plus;12&beta;&ExponentialE;12I&pi;&alpha;&beta;22&alpha;&plus;&beta;&gamma;&plus;&gamma;2&plus;4q&gamma;&plus;&alpha;&plus;&beta;&comma;Anda&equals;1&comma;0&Im;z&comma;HeunGa&comma;q&comma;&alpha;&comma;&beta;&comma;&gamma;&comma;&delta;&comma;z&equals;hypergeom12&beta;&comma;12&alpha;&comma;&gamma;&comma;zz2&comma;Anda&equals;2&comma;q&equals;&alpha;&beta;&comma;&delta;&equals;&alpha;&plus;&beta;2&gamma;&plus;1&comma;1<&real;z&comma;&real;z<1&comma;HeunGa&comma;q&comma;&alpha;&comma;&beta;&comma;&gamma;&comma;&delta;&comma;z&equals;hypergeom13&beta;&comma;13&alpha;&comma;12&comma;14zz32&comma;Anda&equals;4&comma;q&equals;&alpha;&beta;&comma;&gamma;&equals;12&comma;&delta;&equals;23&alpha;&plus;23&beta;&comma;1<&real;z&comma;&real;z<1&comma;&pi;<&Im;z&comma;&Im;z<&pi;

(6)

When α=n, with n a positive integer, the nth + 1 coefficient in the series expansion is a polynomial in q of order n&plus;1. If q is a root of that polynomial, that coefficient is zero and with it all the following ones, so the series truncates and HeunG is a polynomial. For example, for α=−1

HGHeunGa&comma;q&comma;1&comma;&beta;&comma;g&comma;&delta;&comma;z

HG:=HeunGa&comma;q&comma;1&comma;&beta;&comma;g&comma;&delta;&comma;z

(7)

QsimplifyseriesHG&comma;z&comma;3&comma;size

Q:=1&plus;qgaz&plus;12q2&plus;g&plus;&delta;a&delta;&plus;&beta;q&plus;&beta;gaga2g&plus;1z2&plus;Oz3

(8)

So the coefficient of z2 is

c2coeffQ&comma;z&comma;2

c2:=12q2&plus;g&plus;&delta;a&delta;&plus;&beta;q&plus;&beta;gaga2g&plus;1

(9)

solving for q, requesting from solve to return using RootOf we have

_EnvExplicitfalse

_EnvExplicit:=false

(10)

q&equals;solvec2&comma;q

q&equals;RootOf_Z2&plus;a&delta;&plus;ag&plus;&beta;&delta;_Z&plus;&beta;ga

(11)

substituting in HG we have

HG_polynomialsubs&comma;HG

HG_polynomial:=HeunGa&comma;RootOf_Z2&plus;a&delta;&plus;ag&plus;&beta;&delta;_Z&plus;&beta;ga&comma;1&comma;&beta;&comma;g&comma;&delta;&comma;z

(12)

When the function admits a polynomial form, as is the case of HG_polynomial by construction, to obtain the actual polynomial of degree n (in this case n&equals;1) use

HeunG&equals;HeunG:-SpecialValues:-Polynomial|HeunG&equals;HeunG:-SpecialValues:-Polynomial

1&plus;RootOf_Z2&plus;a&delta;&plus;ag&plus;&beta;&delta;_Z&plus;&beta;gazga

(13)

References

  

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.

See Also

FunctionAdvisor

Heun

HeunB

HeunC

HeunD

HeunT

hypergeom

 


Download Help Document

Was this information helpful?



Please add your Comment (Optional)
E-mail Address (Optional)
What is ? This question helps us to combat spam