This presentation is a continuation of the reconsideration of solutions to second-order linear differential equations with polynomial coefficients. The general form of the equation is

${A}_{2}\left(z\right)\frac{{d}^{2}F}{d{z}^{2}}+{A}_{1}\left(z\right)\frac{dF}{dz}+{A}_{0}\left(z\right)F=0\phantom{\rule{5em}{0ex}}{A}_{i}=\sum _{k}{a}_{ik}{z}^{k}$

Solutions to cases of this equation are controlled by the singularity structure of the equation, namely whether singular points exist and what the s-rank of each singular point is. Regular singular points are all assigned an s-rank of unity. The s-rank of an irregular finite singular point is given by

$R\left({z}_{\mathrm{s}}\right)=max[{s}_{2}-{s}_{1}\phantom{\rule{.5em}{0ex}},\phantom{\rule{.5em}{0ex}}\frac{{s}_{2}-{s}_{0}}{2}]$

where *s*_{2}, *s*_{1} and *s*_{0} designate the lowest powers of the independent variable with respect to the singularity in each coefficient polynomial. The s-rank of an irregular singular point at infinity is given by

$R\left(\infty \right)=max[{r}_{1}-{r}_{2}+2\phantom{\rule{.5em}{0ex}},\phantom{\rule{.5em}{0ex}}\frac{{r}_{0}-{r}_{2}+4}{2}]$

where *r*_{2}, *r*_{1} and *r*_{0} designate the highest powers of the independent variable in each coefficient polynomial.

The general process of categorizing cases of the equation can be extended to coefficient polynomials of arbitrary degree, but in this part of the presentation will be limited to constant, linear or quadratic polynomials. As seen at the end of the first part of the presentation, a quadratic coefficient on the first derivative of the equation is already enough complexity to go beyond functions in common use in physics.

It should be noted that the two coefficient polynomials

${a}_{0}+{a}_{1}z+{a}_{2}{z}^{2}\phantom{\rule{5em}{0ex}}{a}_{0}+{a}_{2}{z}^{2}$

make the same contribution to the values of s-ranks, since the formulae for the values depend on the lowest and highest powers at the singularities. That means that solutions to equations with coefficients of the second form can always be found from those with coefficients of the first form by setting *a*_{1} equal to zero.

A tabulation of the singularity structure of each case of constant, linear and quadratic coefficients is available here. The term “multisymbol” denotes the collection of s-ranks for all singular points. They are first given by increasing complexity of the coefficients, from constant on the second derivative through quadratics on all derivatives. The multisymbols are then grouped by identical symbols: these equations will all have a similar solution.

The tabulation includes sets of transformed multisymbols, for a change to a squared variable, a cubed variable, a square root variable, a cube root variable and a cubed square root variable. Three of these are transformations considered in another presentation. These three transformations will apply to only one finite singular point at a time, leaving any others unaffected.

The remaining transformations are simple enough to describe that no tabulation is needed. Splitting off a power from the solutions will at most add a regular singular point at either a finite location or infinity, but not change the singularity structure otherwise. Splitting off an exponential of the form
${e}^{{z}^{k}}$
will not change finite singular points, but will add an irregular singular point at infinity of s-rank *k*+1 unless there is already something equivalent or worse there.

Keep in mind when using this tabulation that the exponents *s*_{1} and *s*_{0} represent factors with respect to the finite singularity. If the polynomial on the second derivative has a leading constant term, then there will be one finite singular point not at the origin for a linear polynomial, and one or two singular points not at the origin for a quadratic polynomial. The exponents for the other coefficients are in that case meant with respect to locations of these finite singular points not at the origin.

Also keep in mind that discriminants are assumed to be nonzero, so that any factored quadratic coefficients will have distinct roots. Cases of zero discriminant need to be addressed separately, and are not included in the tabulation.

The tabulation is useful in several ways. It first shows the number of different types of special functions that can occur in solutions to these differential equations, where each unique multisymbol corresponds to a different function. After identifying the multisymbols for known special functions, it is then simple to know which equations require functions more complicated than usually encountered for their general solutions. The sets of transformed multisymbols indicate how a change of variable can cast a given equation in terms of a known equation, and explain the connections previously found between solutions.

First consider identification of multisymbols in terms of known functions. It was pointed out in the first part of the presentation that the second-order equation with constant coefficients

${a}_{20}\frac{{d}^{2}F}{d{z}^{2}}+{a}_{10}\frac{dF}{dz}+{a}_{00}F=0$

has a multisymbol

${a}_{22}{z}^{2}\frac{{d}^{2}F}{d{z}^{2}}+{a}_{11}z\frac{dF}{dz}+{a}_{00}F=0$

has the multisymbol

Now consider three kinds of hypergeometric functions: Gauss, confluent and reduced confluent. The differential equations for these three functions are

$z(1-z)\frac{{d}^{2}F}{d{z}^{2}}+[c-(a+b+1\left)z\right]\frac{dF}{dz}-abF=0$

$z\phantom{\rule{.2em}{0ex}}\frac{{d}^{2}F}{d{z}^{2}}+(c-z)\frac{dF}{dz}-aF=0$

$z\phantom{\rule{.2em}{0ex}}\frac{{d}^{2}F}{d{z}^{2}}+c\frac{dF}{dz}-F=0$

and the solutions are denoted by
${}_{2}F_{1}(a,b;c;z)$ ,
${}_{1}F_{1}(a,c;z)$
and
${}_{0}F_{1}(c;z)$ ,
respectively. A second linearly independent solution to each equation is found by setting
*k* to remove a resulting coefficient on each third term inverse in the independent variable, and identifying the new parameters for each function. The result is that these equations have similar pairs of linearly independent functions:

${}_{2}F_{1}(a,b;c;z)\phantom{\rule{4em}{0ex}}{z}^{1-c}{}_{2}F_{1}(a-c+1,b-c+1;2-c;z)$

${}_{1}F_{1}(a,c;z)\phantom{\rule{4em}{0ex}}{z}^{1-c}{}_{1}F_{1}(a-c+1,2-c;z)$

${}_{0}F_{1}(c;z)\phantom{\rule{4em}{0ex}}{z}^{1-c}{}_{0}F_{1}(2-c;z)$

The equation for the Gauss hypergeometric function has regular singular points at zero and one, as well as a regular singular point at infinity, and thus has the multisymbol

The equation for the confluent hypergeometric function has a regular singular point at the origin and an irregular singular point at infinity, with the multisymbol

{ L , L , C } , { L , C , L } and { L , L , L }

of similar solutions identified in the previous presentation. The general solution includes in each case one part with an explicit power as a factor, but since the origin is already a regular singular point the process of splitting off a power does not change the singularity structure. The various solutions also have different exponential factors, but since the singular point at infinity is already irregular with the same s-rank as an exponential function the process of splitting off an exponential does not change the singularity structure either.

This last statement has a direct implication: the confluent hypergeometric function has a simple identity that does not occur for the other hypergeometric functions. Setting

$z\phantom{\rule{.2em}{0ex}}\frac{{d}^{2}f}{d{z}^{2}}+(c+z)\frac{df}{dz}+(c-a)f=0$

The solution of this last equation is the basis for an identity that allows one to change the sign of the independent variable in this function by multiplying by an exponential:

${e}^{z}{}_{1}F_{1}(b,c;-z)={}_{1}F_{1}(c-b,c;z)$

This identity occurs because the coefficients of the first and second terms of the equation are both linear, and splitting off the exponential does not alter the form of the equation. Another way of saying this is that splitting off the exponential does not alter the singularity structure, since the singular point at infinity already has the same s-rank as the exponential.

There is another special relationship between the confluent hypergeometric function and the exponential: when the two parameters of the former are equal, its differential equation becomes

$z\phantom{\rule{.2em}{0ex}}\frac{{d}^{2}F}{d{z}^{2}}+(c-z)\frac{dF}{dz}-cF={z}^{1-c}\frac{d}{dz}\left[{z}^{c}\right(\frac{dF}{dz}-F\left)\right]=0$

and one of the solutions is an exponential itself. This does not mean, however, that the singularity structure of the general solution has changed. The second linearly independent solution is a power and an exponential multiplying an incomplete gamma function, so that the regular singular point at the origin is still there.

While the exponential function shares one equal s-rank with the confluent hypergeometric function, the former does not have a singular point at the origin, and this is why the solution to an equation

The equation for the reduced confluent hypergeometric function has a regular singular point at the origin and an irregular singular point at infinity, but with the multisymbol

There are two ways to attempt to find the relationship. The process of splitting off an exponential from the reduced confluent hypergeometric function would change its multisymbol to *a*_{11} and *a*_{01} equal to zero, a subcase that invalidates the solution found before. This explains why certain constants could not simply be set equal to zero in many of the previous solutions: that would alter the singularity structure of the equation, necessitating a different function for its solution.

The apparent relationship just described is a caution that while multisymbols indicate differences between differential equations, there may be circumstances where they do fail to do so. This happens whenever there is a transformation that simplifies the given equation, reducing the s-rank of one or more singular points. In this particular example, one says that an equation with the multisymbol

The other way to find the relationship is to change to a square root variable, which also alters the multisymbol to

$u\phantom{\rule{.2em}{0ex}}\frac{{d}^{2}F}{d{u}^{2}}+(2c-u)\frac{dF}{du}-4uF=0$

As for the previous general solution, this equation can be brought into the form of the confluent equation by including an exponential factor. The final result is the identity

${}_{0}F_{1}(c;z)={e}^{-2\sqrt{z}}{}_{1}F_{1}(c-\frac{1}{2}\phantom{\rule{.3em}{0ex}},\phantom{\rule{.3em}{0ex}}2c-1\phantom{\rule{.3em}{0ex}};\phantom{\rule{.3em}{0ex}}4\sqrt{z})$

An equation of the form

The next group of solutions to understand is

{ C , L , C } , { C , L , L } , { C , C , Q } and { C , L , Q }

The equations all have the multisymbol *Special Functions: A Unified Theory Based on Singularities*, but this does not yet appear to be standard usage. Since this volume is an extremely well-organized enumeration of generalized special functions, following the usage of what is sure to become a standard reference seems prudent.

The biconfluent hypergeometric equation, as given on page 61 of the book,

$\frac{{d}^{2}F}{d{z}^{2}}-z\phantom{\rule{.2em}{0ex}}\frac{dF}{dz}-aF=0$

also goes by the name of the Hermite differential equation. Its solution, along with solutions of more complicated equations to follow, will be denoted as

The biconfluent hypergeometric function bears a relationship to the confluent function similar to that between the reduced confluent and confluent, but this time with change to a squared variable that alters the multisymbol to

$u\phantom{\rule{.2em}{0ex}}\frac{{d}^{2}F}{d{u}^{2}}+(\frac{1}{2}-\frac{u}{2})\frac{dF}{du}-\frac{a}{4}F=0$

which is a confluent hypergeometric equation as expected, so that one can make the immediate identification

${F}^{\{;3\}}(a;z)={}_{1}F_{1}(\frac{a}{2}\phantom{\rule{.2em}{0ex}},\phantom{\rule{.2em}{0ex}}\frac{1}{2}\phantom{\rule{.2em}{0ex}};\phantom{\rule{.2em}{0ex}}\frac{{z}^{2}}{2})$

This relationship explains why the previous general solutions to the equations in this group all have a squared argument in the confluent hypergeometric function, and why all have a parameter of one half. It is now clear why the two constants *a*_{11} and *a*_{02} cannot both be set equal to zero for the equation

Before considering the last case

${a}_{22}(z+\frac{{a}_{21}}{2{a}_{22}}{)}^{2}\frac{{d}^{2}F}{d{z}^{2}}+({a}_{10}+{a}_{11}z)\frac{dF}{dz}+{a}_{00}F=0$

This equation has a double singularity at
$z=-\frac{{a}_{21}}{2{a}_{22}}$
where the two previous singularities have converged. This is what is meant by a confluent equation: one where singularities flow together. The confluence of the singularities has turned them from two regular singularities into one irregular singularity. There is additionally a regular singularity at infinity, so that the multisymbol is

Now on to the equation

$\frac{{d}^{2}F}{d{z}^{2}}-z\phantom{\rule{.2em}{0ex}}F=0$

which is more commonly known as the Airy differential equation. The equation has one irregular singular point at infinity and a multisymbol of

While the Airy functions are well-enough known in physics to count among standard functions, they too have a relationship with confluent hypergeometric functions. A change to a cubed square root variable alters the multisymbol to

$u\phantom{\rule{.2em}{0ex}}\frac{{d}^{2}F}{d{u}^{2}}+\frac{1}{3}\frac{dF}{du}-\frac{4u}{9}F=0$

As before this equation can be brought into the form of the confluent equation by including an exponential factor, with the resulting identities

${F}^{\{;2.5\}}\left(z\right)=exp(-\frac{2}{3}{z}^{3/2}){}_{1}F_{1}(\frac{1}{6}\phantom{\rule{.3em}{0ex}},\phantom{\rule{.3em}{0ex}}\frac{1}{3}\phantom{\rule{.3em}{0ex}};\phantom{\rule{.3em}{0ex}}\frac{4}{3}{z}^{3/2})={}_{0}F_{1}(\frac{2}{3}\phantom{\rule{.3em}{0ex}};\phantom{\rule{.3em}{0ex}}\frac{{z}^{3}}{9})$

The outer equality can also be determined directly with a change to an cubed variable by setting

These relationships explain the unusual argument in the general solution of the equation

Since the Airy functions do not have finite singularities with multisymbol *z* appears in the Airy functions, instead of the cube-root branch point . Similarly, the confluent hypergeometric function
*z*. The identical circumstance occurs for the biconfluent hypergeometric function, so that while
*z*. Nifty!

Identification of the remaining multisymbols will require the introduction of functions beyond those commonly known. Begin by restating the equation

${a}_{20}\frac{{d}^{2}F}{d{z}^{2}}+({a}_{10}+{a}_{11}z+{a}_{12}{z}^{2})\frac{dF}{dz}+{a}_{00}F=0$

This is a case of what Slavyanov and Lay call the triconfluent Heun equation. There is one irregular singular point at infinity and the multisymbol is

${a}_{20}\frac{{d}^{2}F}{d{z}^{2}}+({a}_{10}+{a}_{11}z+{a}_{12}{z}^{2})\frac{dF}{dz}+({a}_{00}+{a}_{01}z+{a}_{02}{z}^{2}+{a}_{03}{z}^{3}+{a}_{04}{z}^{4})F=0$

because both of these equations have the same singularity structure. This structure will not be altered by splitting off an exponential with terms up to a cubic, so that the transformation

$\frac{{d}^{2}F}{d{z}^{2}}+({b}_{10}+{b}_{11}z+{b}_{12}{z}^{2})\frac{dF}{dz}+({b}_{00}+{b}_{01}z)F=0$

First change the sign of the independent variable with $u=-z$ ,

$\frac{{d}^{2}F}{d{u}^{2}}-({b}_{10}-{b}_{11}u+{b}_{12}{u}^{2})\frac{dF}{du}+({b}_{00}-{b}_{01}u)F=0$

then factor the coefficient of the second term:

${b}_{10}-{b}_{11}u+{b}_{12}{u}^{2}={b}_{12}(u-{u}_{+})(u-{u}_{-})\phantom{\rule{3em}{0ex}}{u}_{\pm}=\frac{{b}_{11}\pm \sqrt{{b}_{11}^{2}-4{b}_{12}{b}_{10}}}{2{b}_{12}}$

Now construct a Möbius transformation to move the roots to two as yet unspecified points equidistant from the origin of the independent variable:

$\begin{array}{c}(u-{u}_{+})(w+{w}_{0})=(u-{u}_{-})(w-{w}_{0})\\ u=\frac{{u}_{+}-{u}_{-}}{2{w}_{0}}w+\frac{{u}_{+}+{u}_{-}}{2}\end{array}$

The differential equation in terms of the new variable is

$\begin{array}{l}\frac{{d}^{2}F}{d{w}^{2}}-{b}_{12}\frac{({u}_{+}-{u}_{-}{)}^{3}}{8{w}_{0}^{3}}({w}^{2}-{w}_{0}^{2})\frac{dF}{dw}\\ \phantom{\rule{5em}{0ex}}+\frac{({u}_{+}-{u}_{-}{)}^{2}}{4{w}_{0}^{2}}[{b}_{00}-{b}_{01}(\frac{{u}_{+}-{u}_{-}}{2{w}_{0}}w+\frac{{u}_{+}+{u}_{-}}{2}\left)\right]F=0\end{array}$

Choosing ${w}_{0}={b}_{12}^{1/3}\frac{{u}_{+}-{u}_{-}}{2}$ to make the overall factor on the second term unity brings this to the form of the triconfluent Heun equation as given on page 104 of Slavyanov and Lay,

$\frac{{d}^{2}F}{d{z}^{2}}-(t+{z}^{2})\frac{dF}{dz}+(\lambda -az)F=0$

and a solution to the equation in question is ${F}^{\{;4\}}[a;t;\lambda ;w(z\left)\right]$ , where in this case

$\begin{array}{c}a=\frac{{b}_{01}}{{b}_{12}}\phantom{\rule{4em}{0ex}}t=-\frac{{b}_{11}^{2}-4{b}_{12}{b}_{10}}{2{b}_{12}^{4/3}}\\ \lambda ={b}_{12}^{-2/3}({b}_{00}-\frac{{b}_{11}{b}_{01}}{2{b}_{12}})\\ w\left(z\right)=-{b}_{12}^{1/3}(z+\frac{{b}_{11}}{2{b}_{12}})\end{array}$

The triconfluent Heun function is not yet commonly known in physics. It occurs, for example, in quantum mechanics for the so-called quartic oscillator, when a fourth-order correction term is added to the harmonic oscillator potential. Mathematica does not currently handle Heun functions other than as Mathieu functions, which Slavyanov and Lay call reduced singly confluent Heun functions. The mathematical package Maple does include most Heun functions, although with deficiencies in numerical evaluation.

The triconfluent Heun function provides a solution to the group of equations

{ C , Q , C }, { C , Q , L } and { C , Q , Q }

since they all have the multisymbol

{ L , Q , C }, { L , Q , L } and { L , Q , Q }

all have the multisymbol

$z\phantom{\rule{.2em}{0ex}}\frac{{d}^{2}F}{d{z}^{2}}+(c-tz-{z}^{2})\frac{dF}{dz}+(\lambda -az)F=0$

The two equations

{ L , C , Q } and { L , L , Q}

both have the multisymbol

$z\phantom{\rule{.2em}{0ex}}\frac{{d}^{2}F}{d{z}^{2}}+c\phantom{\rule{.2em}{0ex}}\frac{dF}{dz}+(\lambda -tz-{z}^{2})F=0$

The five equations

{ Q , C , Q }, { Q , L , Q }, { Q , Q , C }, { Q , Q , L } and { Q , Q , Q }

all have the multisymbol

$z(z-1)\frac{{d}^{2}F}{d{z}^{2}}+\left[c\right(z-1)-tz(z-1)+dz]\frac{dF}{dz}+(\lambda -taz)F=0$

The two equations

{ Q , C , L } and { Q , L , L }

both have the multisymbol

$z(z-1)\frac{{d}^{2}F}{d{z}^{2}}+\left[c\right(z-1)+dz]\frac{dF}{dz}+(\lambda -tz)F=0$

In enumerating all of the groups, it has been assumed that quadratic coefficient has zero discriminant and thus two distinct roots. The cases of two equal roots are covered by the multisymbols

${z}^{2}\frac{{d}^{2}F}{d{z}^{2}}+(t+cz-{z}^{2})\frac{dF}{dz}+(\lambda -az)F=0$

The multisymbol

${a}_{22}{z}^{2}\frac{{d}^{2}F}{d{z}^{2}}+{a}_{10}\frac{dF}{dz}+{a}_{01}zF=0$

Changing to an inverse variable with $z=\frac{1}{w}$ , the equation becomes

${a}_{22}{w}^{2}\frac{{d}^{2}F}{d{w}^{2}}+(2{a}_{22}w-{a}_{10}{w}^{2})\frac{dF}{dw}+\frac{{a}_{01}}{w}F=0$

The linear part of the coefficient on the second term can be removed by splitting off a power factor with an exponent of minus one-half, which will add constant and linear parts to the coefficient on the third term. After rescaling the independent variable, one will have a case of the reduced doubly confluent Heun equation as given on page 107 of Slavyanov and Lay,

${z}^{2}\frac{{d}^{2}F}{d{z}^{2}}-{z}^{2}\frac{dF}{dz}+(\lambda -az-\frac{t}{z})F=0$

and the multisymbol

The tabulations of transformed multisymbols can be used to indicate the possibility of additional identities, not only among the Heun functions but also between Heun and hypergeometric functions. As with previous identities, a function with fewer arbitrary parameters corresponds to one with fixed values for one or more of the parameters. The relationships among Heun functions are only simple for fixed values of certain parameters, and so are more restricted in this sense than previous identities.

For example, consider a change to a square root variable in the confluent hypergeometric equation, which alters its multisymbol from

$u\phantom{\rule{.2em}{0ex}}\frac{{d}^{2}F}{d{u}^{2}}+(2c-1-2{u}^{2})\frac{dF}{du}-4auF=0$

This is a form of the Heun biconfluent equation with both of the parameters *t* and λ equal to zero. If the function with all parameters nonzero is notated by

${}_{1}F_{1}(a,c;z)={F}^{\{1;2\}}(a,c;z)={F}^{\{1;3\}}[2a,2c-1;0;0;(2z{)}^{1/2}]$

which means that the confluent hypergeometric function is a case of the biconfluent Heun function with a square root argument. Since the other types of confluent hypergeometric can be written in terms of the confluent function, they can also be written as biconfluent Heun functions with appropriate arguments. In particular, multiplying the biconfluent hypergeometric equation by a single power of the independent variable gives a case of the biconfluent Heun equation, so one can immediately identify

${F}^{\{;3\}}(a;z)={F}^{\{1;3\}}(a,0;0;0;z)$

which is consistent with the previous identification between the two hypergeometric functions.

Since the confluent hypergeometric function has a simple relation to the biconfluent Heun function, it is expected that the reduced confluent hypergeometric function will have a simple relation to the reduced biconfluent Heun function. Inspecting available transformations, a change to a cube root variable alters the multisymbols appropriately. With

$u\phantom{\rule{.2em}{0ex}}\frac{{d}^{2}F}{d{u}^{2}}+(3c-2)\frac{dF}{du}-9{u}^{2}F=0$

which is a form of the reduced biconfluent Heun equation with both of the parameters *t* and λ equal to zero. With notation as before one can write the identity

${}_{0}F_{1}(c;z)={F}^{\{1;1.5\}}(c;z)={F}^{\{1;2.5\}}[3c-2;0;0;(9z{)}^{1/3}]$

The reduced confluent hypergeometric function is related to the confluent hypergeometric function by an exponential and a square root argument. Comparing the identities just found, the reduced biconfluent Heun function will be equal to the product of an exponential with argument of exponent three halves and a biconfluent Heun function with argument of exponent three quarters. With
*t* = λ = 0

$u\phantom{\rule{.2em}{0ex}}\frac{{d}^{2}F}{d{u}^{2}}+\frac{(4c-1)}{3}\frac{dF}{du}-\frac{16}{9}{u}^{3}F=0$

For a solution of the form $F\left(u\right)={e}^{-k{u}^{2}}f\left(u\right)$ the equation is then

$u\phantom{\rule{.2em}{0ex}}\frac{{d}^{2}F}{d{u}^{2}}+[\frac{(4c-1)}{3}-4k{u}^{2}]\frac{dF}{du}+[4{k}^{2}{u}^{3}-\frac{16}{9}{u}^{3}-\frac{4(2c+1)}{3}ku]F=0$

and for an appropriate choice of the free constant one has the identity

${F}^{\{1;2.5\}}[c;0;0;z]=exp(-\frac{2}{3}{z}^{3/2}){F}^{\{1;3\}}[\frac{2c+1}{3}\phantom{\rule{.2em}{0ex}},\phantom{\rule{.2em}{0ex}}\frac{4c-1}{3}\phantom{\rule{.2em}{0ex}};0;0;\sqrt{\frac{8}{3}}\phantom{\rule{.3em}{0ex}}{z}^{3/4}]$

which is completely consistent with previous identities. If the parameters *t* and λ were not set equal to zero, then the coefficient of the third term in the reduced biconfluent Heun equation would include

Since the triconfluent Heun equation appears readily enough in physics, it is interesting to develop its identities to the hypergeometric functions and other Heun functions. This will be done for the restricted case *t* = λ = 0

$\frac{{d}^{2}F}{d{z}^{2}}-{z}^{2}\frac{dF}{dz}-azF=0$

A change to a cubed variable alters the multisymbol from

$u\frac{{d}^{2}F}{d{u}^{2}}+(\frac{2}{3}-\frac{u}{3})\frac{dF}{du}-\frac{a}{9}F=0$

and one has the immediate identity

${F}^{\{;4\}}(a;0;0;z)={}_{1}F_{1}(\frac{a}{3}\phantom{\rule{.3em}{0ex}},\phantom{\rule{.3em}{0ex}}\frac{2}{3}\phantom{\rule{.3em}{0ex}};\phantom{\rule{.3em}{0ex}}\frac{{z}^{3}}{3})$

Knowing the corresponding linearly independent hypergeometric function indicates that the second linearly independent solution to this case of the triconfluent Heun equation is

The confluent hypergeometric function in this identity can represent a reduced confluent hypergeometric function when the first parameter is half of the second. With the restriction *a* = 1

${F}^{\{;4\}}(1;0;0;z)=exp\left(\frac{{z}^{3}}{6}\right){}_{0}F_{1}(\frac{5}{6}\phantom{\rule{.3em}{0ex}};\phantom{\rule{.3em}{0ex}}\frac{{z}^{6}}{144})$

The reduced confluent hypergeometric function can represent an Airy function when its parameter is equal to two thirds or four thirds, which is not the case here. The triconfluent Heun function apparently cannot be directly related to Airy functions.

Finally, consider relationships of the triconfluent Heun function to the biconfluent and reduced biconfluent Heun functions, again all for *t* = λ = 0

$u\frac{{d}^{2}F}{d{u}^{2}}+(\frac{1}{3}-\frac{2}{3}{u}^{2})\frac{dF}{du}-\frac{4}{9}auF=0$

and the identity between the two functions is

${F}^{\{;4\}}[a;0;0;z]={F}^{\{1;3\}}[\frac{2}{3}a\phantom{\rule{.2em}{0ex}},\phantom{\rule{.2em}{0ex}}\frac{1}{3}\phantom{\rule{.2em}{0ex}};0;0;\sqrt{\frac{2}{3}}\phantom{\rule{.3em}{0ex}}{z}^{3/2}]$

The biconfluent Heun function here can only represent a reduced biconfluent Heun function when the parameters have the relative values
*a* = 1

A change to a squared variable alters the triconfluent multisymbol from

$\frac{{d}^{2}f}{d{z}^{2}}+[(1-a)z-\frac{1}{4}{z}^{4}]f=0$

Now with

$u\frac{{d}^{2}f}{d{u}^{2}}+\frac{1}{2}\frac{df}{du}+\frac{1}{4}[(1-a)\sqrt{z}-\frac{1}{4}{u}^{2}]f=0$

and it is clear that *a* = 1

${F}^{\{;4\}}[1;0;0;z]=exp\left(\frac{{z}^{3}}{6}\right){F}^{\{1;2.5\}}[\frac{1}{2}\phantom{\rule{.2em}{0ex}};0;0;\frac{{z}^{2}}{{16}^{1/3}}]$

which is completely consistent with the identity between the biconfluent and reduced biconfluent Heun functions.

Summary table of distinct multisymbols for second-order differential equations with polynomial coefficients up to quadratic and their corresponding general solution functions:

Multisymbol | Solution Function |

{ ; 2 } | Exponential |

{ ; 2.5 } | Reduced biconfluent hypergeometric
or Airy aka confluent hypergeometric with cubed square root argument aka reduced confluent hypergeometric with cubed argument |

{ ; 3 } | Biconfluent hypergeometric or Hermite aka confluent hypergeometric with squared argument aka biconfluent Heun |

{ ; 4 } | Triconfluent Heun |

{ 1 ; 1 } | Single power |

{ 1 ; 1.5 } | Reduced confluent hypergeometric aka confluent hypergeometric with square root argument aka biconfluent Heun with cube root argument |

{ 1 ; 2 } | Confluent hypergeometric aka biconfluent Heun with square root argument |

{ 1 ; 2.5 } | Reduced biconfluent Heun |

{ 1 ; 3 } | Biconfluent Heun |

{ 1 , 1 ; 1 } | Gauss hypergeometric |

{ 1 , 1 ; 1.5 } | Reduced singly confluent Heun |

{ 1 , 1 ; 2 } | Singly confluent Heun |

{ 2 ; 1 } | Confluent hypergeometric with inverse argument |

{ 2 ; 1.5 } | Reduced doubly confluent Heun with inverse argument |

{ 2 ; 2 } | Doubly confluent Heun |

The solutions can be accompanied by powers and exponentials that do not alter the singularity structure.

*Uploaded 2013.07.31 — Updated 2013.08.07*
analyticphysics.com