return the nth order ODE having a given hypergeometric pFq or MeijerG function as solution
hyperode(H, y(x), n, [S1, S2, ..., Sn])
hypergeometric function; can be hypergeom⁡...,...,f⁡x or MeijerG⁡...,...,f⁡x, or a product of the form g⁡x⁢H
dependent variable of the ODE; can be any unknown function of one variable
(optional) integer; differential order of the ODE to be returned
[S1, ..., Sn]
(optional) list with expressions depending on x; solutions of the ODE to be returned
Hypergeometric functions pFq⁡a1,...,ap,b1,...,bq,f⁡x, or their generalization known as MeijerG functions, can be defined in different ways. One possible way is to define them as the solutions to some linear ODEs. In connection with that approach, the hyperode routine receives as input a hypergeometric function H as a parameter, where H is constructed using one of hypergeom or MeijerG, as either
H([...], [...], f(x)), where f⁡x is an arbitrary expression depending on x, or
the product of an arbitrary expression g⁡x times H([...], [...], f(x)).
Using this input, hyperode returns the general nth order ODE having as a solution that hypergeometric function H or that product of H times g⁡x.
By identifying the general ODE underlying the definition of a hypergeometric function, this command can be of help when studying properties of these relevant functions and facilitate the understanding of related algorithms. In other cases, this ODE representation of pFq or MeijerG opens the way to computations (for example, differential elimination) that can only be performed with polynomial differential objects. See dpolyform, casesplit, DifferentialAlgebra, and rifsimp.
The ODE returned by hyperode, having a solution y⁡x=g⁡x⁢pFq⁡a1,...,ap,bi,...,bq,f⁡x (where pFq is expressed using hypergeom) is built by using the formula
where D is the differential operator ddx, followed by a change of variables
and finally by the renaming of the variables according to
The output in this case is a linear ODE where the differential order is the largest of p and q+1.
The ODE returned by hyperode, when the input is a MeijerG generalized hypergeometric function, is built using the same process, but the formula used in the first step is:
where in above we assume, without loss of generality, that p is less than or equal to q, and ai and bi represent, respectively, the parameters entering the first and second lists of parameters in MeijerG. The output in this case is a qth order linear ODE.
When the differential order of the ODE to be returned is not indicated, the one implied by the given pFq or MeijerG function is used. When the differential order is indicated as an extra argument, and provided it is greater than the one implied by the hypergeometric function, the ODE is constructed in the same way as just outlined except that before proceeding, the lists [a,...], [b,...] are augmented by introducing into both of them the necessary number of additional parameters, all equal to zero.
Since the transformation shown above with arbitrary f⁡x and g⁡x represents the structure invariance group for the linear ODEs (that is, it is enough to map any linear ODE into any other one of the same differential order), then by keeping both f⁡x and g⁡x arbitrary, the ODE returned is "the most general linear ODE of a given order, written in such a way that its solution is expressed in terms of the product of a hypergeometric function (pFq or MeijerG) times an arbitrary function g(x)." Therefore, any linear ODE and its solution in terms of a hypergeometric function can be obtained from this general form by appropriately choosing f and g (see the last example).
When an optional list of additional ODE solutions [S1x,S2x,...,Snx] is given, the returned ODE is built by first using the formula explained above, and then applying to it operators of the form
for each Sn⁡x, where the functions Fn⁡x are adjusted so that in addition to
for all Sn⁡x,
is also a solution of the returned ODE.
Note: This building process raises the differential order of the returned ODE by n.
y⁡x⁢will now be displayed as⁢y
derivatives with respect to⁢x⁢of functions of one variable will now be displayed with '
The general linear ODE family having as solution "2F1" (two indices in the first list and one index in the other list)
is of second order since the number of operands of the first list is two, and the number of operands of the second list plus one is not greater than two. This ODE is known as the hypergeometric ODE (or Gauss ODE) and is given by
The meaning of "1F0" can be determined by looking at the solution of the corresponding linear ODE (which is of order max( 1, 0+1 ) = 1).
Hence, up to a constant c__1 (which may depend on a), 1F0 is equal to the right-hand side of the ODE solution above. An alternative approach to this result is to first convert 1F0 to an infinite sum and then perform the summation:
Equivalently, but perhaps a more black-box approach, would be a one-step conversion.
Consider the following statements in order to determine the meaning of the "0F0" hypergeometric function.
In particular, consider the following "0F1",
or use the one-step conversion routine.
The general third order ODE equivalent to the second order hypergeometric ODE is obtained by departing from the general hypergeometric function: "3F2", in turn obtained from 2F1 by adding one arbitrary parameter to each of the two lists.
Test these results by using odetest.
The next third order ODE has as solutions both
hyperode can be used in the same way with (generalized hypergeometric) MeijerG functions as input
Comparing this output with the one for a pFq function
it is apparent that, for some values of a, b, c, and d, there is a relation between the pFq and the more general MeijerG functions. This relation is actually given by:
Finally, the most general second order linear ODE, can be written in terms of two functions
f⁡x⁢will now be displayed as⁢f
g⁡x⁢will now be displayed as⁢g
such that the ODE solution is written in terms of f and g as:
sol ≔ y⁡x=g⁡x⁢hypergeom⁡a,b,c,f⁡x
Such a form of the general 2nd order linear ODE is obtained as follows:
ode ≔ collect⁡isolate⁡,ⅆ2ⅆx2⁢y⁡x,diff,normal
Consider for instance some particular values of f and g.
The following is the solution to this ODE.
This is an example of an ODE family with radicals and its solution in terms of hypergeometric functions.
Marsden, J.E.; Sirovich, L.; and Antman, S.S. eds. Texts in Applied Mathematics. 56 vols. New York: Springer-Verlag, 1991. Vol. 8: Hypergeometric Functions and Their Applications.
Mathai, A.M. A Handbook of Generalized Special Functions for Statistical and Physical Sciences. Oxford: Clarendon Press, 1993.
Download Help Document