LREtools[HypergeometricTerm]
HypergeometricSolution
return the hypergeometric solution of linear difference equation depending on a hypergeometric term
Calling Sequence
Parameters
Description
Examples
HypergeometricSolution(eq, var)
HypergeometricSolution(eq, var, term, cert)
eq
-
linear difference equation depending on a hypergeometric term
var
function variable for which to solve, for example, z(n)
term
hypergeometric term
cert
certificate of the right-hand side
The HypergeometricSolution(eq, var, term, cert) command returns the hypergeometric 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.
In the current state, the procedure handles only equations with hypergeometric right-hand sides. Here the term "hypergeometric" means that for the right-hand side R⁡x,t⁡x the ratio R⁡x+1,t⁡x+1R⁡x is a rational function in x and t (such a ratio is referred to as "certificate"). For the two argument case the procedure computes and checks the certificate. For the four argument case the certificate should be given directly in cert parameter and equation eq should be homogeneous.
The search for a hypergeometric solution is reduced to searching a rational solution by change of the function variable y⁡x=z⁡x⁢R⁡x,t⁡x, where R⁡x,t⁡x is the right-hand side and z⁡x is a new function variable that should be rational.
The solution is the function, corresponding to var. The solution may involve arbitrary constants of the form, for example, _c1 and _c2.
with⁡LREtoolsHypergeometricTerm:
eq≔y⁡x+1−y⁡x−2x+1x+1⁢x2−5⁢2x+1x+1⁢x−27⁢2x+1x+1+x2⁢2xx+3⁢2xx⁢x+23⁢2xx
sol≔HypergeometricSolution⁡eq,y⁡x
sol≔x2+3⁢x+23⁢2x+1x+1⁢x2+5⁢2x+1x+1⁢x+27⁢2x+1x+1−x2⁢2xx−3⁢2xx⁢x−23⁢2xx2⁢t⁢x2+10⁢x⁢t−x2+54⁢t−3⁢x−23,t,4
IsSolution⁡sol,eq,y⁡x
true
simplify⁡SubstituteTerm⁡sol,x
2xx⁢x2+3⁢x+23
eq≔2x⁢y⁡x+1−y⁡x2x−2x⁢x+12⁢2x+1x+1+x2⁢2xx2x
sol≔x2⁢2x+1x+1⁢2x2⁢x2+2⁢2x+1x+1⁢2x2⁢x+2x2⁢2x+1x+1−x2⁢2xx2⁢t2⁢x2+4⁢t2⁢x+2⁢t2−x2,t,4
x2⁢2xx
eq0≔t⁢y⁡x+1−y⁡xt
cert≔t2⁢32⁢t4⁢x2+128⁢t4⁢x+128⁢t4−x2−2⁢x−1−x2+2⁢t4⁢x2+4⁢t4⁢x+2⁢t4
cert≔t2⁢32⁢t4⁢x2+128⁢t4⁢x+128⁢t4−x2−2⁢x−12⁢t4⁢x2+4⁢t4⁢x+2⁢t4−x2
sol≔HypergeometricSolution⁡eq0,y⁡x,t,2,cert
sol≔x22⁢t4⁢x2+4⁢t4⁢x+2⁢t4−x2,2⁢t2⁢32⁢t4⁢x2+128⁢t4⁢x+128⁢t4−x2−2⁢x−12⁢t4⁢x2+4⁢t4⁢x+2⁢t4−x2,t,2
eq≔eq0−Product⁡subs⁡x=n,t=2n,cert,n=0..x−1
eq≔t⁢y⁡x+1−y⁡xt−∏n=0x−1⁡2n2⁢32⁢n2⁢2n4+128⁢n⁢2n4+128⁢2n4−n2−2⁢n−12⁢n2⁢2n4+4⁢n⁢2n4+2⁢2n4−n2
s≔sol11⁢Product⁡subs⁡x=n,t=2n,sol12,n=0..x−1
s≔x2⁢∏n=0x−1⁡2⁢2n2⁢32⁢n2⁢2n4+128⁢n⁢2n4+128⁢2n4−n2−2⁢n−12⁢n2⁢2n4+4⁢n⁢2n4+2⁢2n4−n22⁢t4⁢x2+4⁢t4⁢x+2⁢t4−x2
IsSolution⁡s,t,2,eq,y⁡x
See Also
LREtools[HypergeometricTerm][RationalSolution]
Download Help Document