SumTools[Hypergeometric]
BottomSequence
bottom sequence of a hypergeometric term
Calling Sequence
Parameters
Description
Examples
References
Compatibility
BottomSequence(T, x, opt)
T
-
hypergeometric term in x
x
name
opt
(optional) equation of the form primitive=true or primitive=false
Consider T as an analytic function in x satisfying a linear difference equation p⁡x⁢T⁡x+1+q⁡x⁢T⁡x=0, where p⁡x and q⁡x are polynomials in x. For h∈ℤ and any integer k, let ck,h be the h-th coefficient of the Laurent series expansion for T at x=k. An integer m is called depth of T if ck,h=0 for all h<m and all integers k, and ck,m≠0 for some k∈ℤ.
The bottom sequence of T is the doubly infinite sequence bx defined as bx=cx,m for all integers x, where m is the depth of T. The command BottomSequence(T, x) returns the bottom sequence of T in form of an expression representing a function of (integer values of) x. Typically, this is a piecewise expression.
The bottom sequence bx is defined at all integers x and satisfies the same difference equation p⁡x⁢bx+1+q⁡x⁢bx=0 as T.
If T is Gosper-summable and S=v⁢T is its indefinite sum found by Gosper's algorithm, then the depth of S is also m. If the optional argument primitive=true (or just primitive) is specified, the command returns a pair v,u, where v is the bottom sequence of T and u is the bottom sequence of S or FAIL if T is not Gosper-summable.
Note that this command rewrites expressions of the form nk in terms of GAMMA functions Γ⁡n+1Γ⁡k+1⁢Γ⁡n−k+1.
If assumptions of the form x0<x and/or x<x1 are made, the depth and the bottom of T are computed with respect to the given interval instead of −∞..∞.
with⁡SumToolsHypergeometric:
T≔n⁢n!
b,s≔BottomSequence⁡T,n,primitive
b,s≔−−1n⁢nΓ⁡−nn≤−100≤n,−−1nΓ⁡−nn≤−100≤n
Note that b is not equivalent to T:
eval⁡b,n=1
0
eval⁡T,n=1
1
eval⁡b,n=−1
−1
eval⁡T,n=−1
Error, numeric exception: division by zero
However, b satisfies the same difference equation as T:
expand⁡n⁢eval⁡T,n=n+1−n+12⁢T
z≔n⁢eval⁡b,n=n+1−n+12⁢b
z≔n⁢−−1n+1⁢n+1Γ⁡−n−1n≤−200≤n+1−n+12⁢−−1n⁢nΓ⁡−nn≤−100≤n
simplify⁡zassumingn≤−2
simplify⁡zassuming0≤n
eval⁡z,n=−1
s is an indefinite sum of b:
z≔eval⁡s,n=n+1−s−b
z≔−−1n+1Γ⁡−n−1n≤−200≤n+1−−−1nΓ⁡−nn≤−100≤n−−−1n⁢nΓ⁡−nn≤−100≤n
Now assume that 0≤n:
b,s≔BottomSequence⁡T,n,primitiveassuming0≤n
b,s≔n⁢Γ⁡n+1,Γ⁡n+1
With that assumption, b and T are equivalent, and s is an indefinite sum of both:
simplify⁡b−T
simplify⁡eval⁡s,n=n+1−s−b
Example of a hypergeometric term with parameters:
T≔Γ⁡−nn−k
BottomSequence⁡T,n
0n≤−1−1n−n+k⁢Γ⁡n+10≤n
Note that k is considered non-integer.
BottomSequence⁡T,nassumingk::nonnegint
BottomSequence⁡eval⁡T,k=2,n
0n≤1−12n=203≤n
T≔binomial⁡2⁢n−3,n4n
T≔2⁢n−3n4n
b,s≔0n≤−112n=0−18n=14−n⁢n−2⁢Γ⁡2⁢n−12⁢Γ⁡n2⁢n2≤n,0n≤012n=14−n⁢n+1⁢Γ⁡2⁢n−1Γ⁡n22≤n
S.A. Abramov, M. Petkovsek. "Analytic solutions of linear difference equations, formal series, and bottom summation." Proc. of CASC'07, (2007): 1-10.
S.A. Abramov, M. Petkovsek. "Gosper's Algorithm, Accurate Summation, and the Discrete Newton-Leibniz Formula." Proceedings of ISSAC'05, (2005): 5-12.
The SumTools[Hypergeometric][BottomSequence] command was introduced in Maple 15.
For more information on Maple 15 changes, see Updates in Maple 15.
See Also
assuming
binomial
SumTools[DefiniteSum][SummableSpace]
SumTools[Hypergeometric][Gosper]
Download Help Document