The mhypergeom_formal_sol command returns formal solutions with m-hypergeometric series coefficients for the given homogeneous linear ordinary differential equation with polynomial coefficients.

If ode is an expression, then it is equated to zero.

The command returns an error message if the differential equation ode does not satisfy the following conditions.

ode must be homogeneous and linear in var

The coefficients of ode must be polynomial in the independent variable of var, for example, $x$, over the rational number field which can be extended by one or more parameters.

A homogeneous linear ordinary differential equation with coefficients that are polynomials in $x$ has a basis of formal solutions (see DEtools[formal_sol]). A formal solution contains a finite number of power series $\sum _{n=0}^{\mathrm{\infty }}v\left(n\right)T^n$ where $T$ is a parameter and the sequence $v\left(n\right)$ satisfies a linear recurrence (homogeneous or inhomogeneous).

This command selects so-called "nice" solutions, that is, solutions that contain series  where $v\left(n+m\right)=p\left(n\right)v\left(n\right)$ for all sufficiently large $n$, where $m\ge 1$ is an integer and $p\left(n\right)$ is a rational function.

The command determines an integer $N\ge 0$ such that $v\left(nm\right)$, $v\left(nm+1\right)$, ..., $v\left(nm+m-1\right)$ can be represented in the form of hypergeometric terms (see SumTools[Hypergeometric], LREtools):

for  $0\le i<m$ and all $n\ge N$.

'parameter'=T

Specifies the name T that is used to denote $\mathrm{\lambda }{x}^{\frac{1}{r}}$ where $\mathrm{\lambda }$ is a constant and $r$ is called the ramification index. If this option is given, then the command expresses the formal solutions in terms of T and returns a list of lists each of which is of the form [formal solution, relation between T and x]. Otherwise, it returns the formal solutions in terms of $x^{\frac{1}{r}}$.

x=a or 'point'=a

Specifies the expansion point a. It can be an algebraic number, depending rationally on some parameters, or $\mathrm{\infty }$.

The default is $a=0$.

'free'=C

Specifies a base name C to use for free variables C[0], C[1], etc. The default is the global name  _C. Note that the number of free variables may be less than the order of the given equation.

'indices'=[n,k]

Specifies names for dummy variables. The default values are the global names _n and _k. The name n is used as the summation index in the power series. The name k is used as the product index in ( * ).

'outputHGT'=name

Specifies the form of representation of hypergeometric terms.  The default value is 'active'.

'inert' - the hypergeometric term ( * ) is represented by an inert product, except for $\prod _{k=N}^{n-1}1$, which is simplified to $1$.

'rcf1' or 'rcf2' - the hypergeometric term is represented in the first or second minimal representation, respectively (see ConjugateRTerm).

'active' - the hypergeometric term is represented by non-inert products which, if possible, are computed (see product).

$\mathrm{with}\left(\mathrm{Slode}\right)\&colon;$

$\mathrm{ode}≔\left(x^2+1\right)x\mathrm{diff}\left(y\left(x\right)\&comma;x\&comma;x\&comma;x\right)+3\left(2x^2+1\right)\mathrm{diff}\left(y\left(x\right)\&comma;x\&comma;x\right)-12y\left(x\right)$

$\mathrm{ode}≔\left(x^2+1\right)x\left(\frac{{\&DifferentialD;}^3}{\&DifferentialD;x^3}y\left(x\right)\right)+3\left(2x^2+1\right)\left(\frac{{\&DifferentialD;}^2}{\&DifferentialD;x^2}y\left(x\right)\right)-12y\left(x\right)$

$\mathrm{mhypergeom\_formal\_sol}\left(\mathrm{ode}\&comma;y\left(x\right)\right)$

$\frac{\left(2x^3+x\right){\mathrm{\_C}}_0+\frac{{\mathrm{\_C}}_1\left(\sum _{\mathrm{\_n}=1}^{\mathrm{\infty }}\frac{\mathrm{\Gamma }\left(\mathrm{\_n}-\frac{3}{2}\right){\left(\mathrm{-1}\right)}^{\mathrm{\_n}}{x}^{2\mathrm{\_n}}}{\mathrm{\Gamma }\left(\mathrm{\_n}\right)}\right)}{2\sqrt{\mathrm{\pi }}}}{x}$