LREtools[HypergeometricTerm]
RationalSolution
return the rational solution of linear difference equation depending on a hypergeometric term
Calling Sequence
Parameters
Description
Examples
References
RationalSolution(eq, var, term)
eq
-
linear difference equation depending on a hypergeometric term
var
function variable for which to solve, for example, z(n)
term
hypergeometric term
The RationalSolution(eq, var, term) command returns the rational solution of the linear difference equation eq. If such a solution does not exist, the function returns NULL.
The hypergeometric term in the linear difference equation is specified by a name, for example, t. The meaning of the term is defined by the parameter term. It can be specified directly in the form of an equation, for example, t=n!, or specified as a list consisting of the name of term variable and the consecutive term ratio, for example, t,n+1.
If the third parameter is omitted, then the input equation can contain a hypergeometric term directly (not a name). In this case, the procedure extracts the term from the equation, transforms the equation to the form with a name representing a hypergeometric term, and then solves the transformed equation.
The term "rational solution" means a solution y⁡x in Q⁡x⁡t. (See PolynomialSolution for the meaning of "polynomial solution".) Here we use the term "denominator" which is q in Qxt to mean that q⁢y is in Qxt,t−1 .
The search for a rational solution is based on finding a universal denominator which is u in Qxt such that u⁢y is in Qxt,t−1 for any rational solution y. By replacing y with Yu in the given equation, we reduce the problem to searching for a polynomial solution.
The solution is the function, corresponding to var. The solution involves arbitrary constants of the form, for example, _c1 and _c2.
with⁡LREtoolsHypergeometricTerm:
eq≔y⁡n+1⁢1+n+1⁢t−y⁡n⁢1+t
RationalSolution⁡eq,y⁡n,t=n!
_C11+t,t,n+1
eq≔numer⁡y⁡n+2⁢n+2+n+2⁢n+1⁢t1+n+2⁢n+1⁢t−y⁡n⁢n+t1+t
eq≔y⁡n+2⁢n2⁢t2−y⁡n⁢n3⁢t−y⁡n⁢n2⁢t2+y⁡n+2⁢n2⁢t+3⁢y⁡n+2⁢n⁢t2−3⁢y⁡n⁢n2⁢t−3⁢y⁡n⁢n⁢t2+4⁢y⁡n+2⁢n⁢t+2⁢y⁡n+2⁢t2−2⁢y⁡n⁢n⁢t−2⁢y⁡n⁢t2+y⁡n+2⁢n+4⁢y⁡n+2⁢t−y⁡n⁢n−y⁡n⁢t+2⁢y⁡n+2
_C1⁢1+tn+t,t,n+1
eq≔y⁡n+1⁢2⁢t+n+1−y⁡n⁢t+n
eq≔y⁡n+1⁢2⁢t+n+1−y⁡n⁢n+t
RationalSolution⁡eq,y⁡n,t=2n
_C1n+t,t,2
eq≔y⁡n+1⁢2n+1+n+1−y⁡n⁢2n+n
RationalSolution⁡eq,y⁡n
_C1t+n,t,2
Abramov, S.A., and Bronstein, M. "Hypergeometric dispersion and the orbit problem." Proc. ISSAC 2000.
Bronstein, M. "On solutions of Linear Ordinary Difference Equations in their Coefficients Field." INRIA Research Report. No. 3797. November 1999.
See Also
LREtools[HypergeometricTerm][HGDispersion]
LREtools[HypergeometricTerm][PolynomialSolution]
LREtools[HypergeometricTerm][RationalSolution]
LREtools[HypergeometricTerm][SubstituteTerm]
LREtools[HypergeometricTerm][UniversalDenominator]
Download Help Document