BottomSequence - Maple Help
For the best experience, we recommend viewing online help using Google Chrome or Mozilla Firefox.

Online Help

All Products    Maple    MapleSim


SumTools[Hypergeometric]

  

BottomSequence

  

bottom sequence of a hypergeometric term

 

Calling Sequence

Parameters

Description

Examples

References

Compatibility

Calling Sequence

BottomSequence(T, x, opt)

Parameters

T

-

hypergeometric term in x

x

-

name

opt

-

(optional) equation of the form primitive=true or primitive=false

Description

• 

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 −∞..∞.

Examples

> 

with⁡SumToolsHypergeometric&colon;

> 

T≔n⁢n!

T≔n⁢n!

(1)
> 

b,s≔BottomSequence⁡T&comma;n&comma;primitive

b,s≔−−1n⁢nΓ⁡−nn≤−100≤n,−−1nΓ⁡−nn≤−100≤n

(2)

Note that b is not equivalent to T:

> 

eval⁡b&comma;n=1

0

(3)
> 

eval⁡T&comma;n=1

1

(4)
> 

eval⁡b&comma;n=−1

−1

(5)
> 

eval⁡T&comma;n=−1

Error, numeric exception: division by zero

However, b satisfies the same difference equation as T:

> 

expand⁡n⁢eval⁡T&comma;n=n+1−n+12⁢T

0

(6)
> 

z≔n⁢eval⁡b&comma;n=n+1−n+12⁢b

z≔n⁢−−1n+1⁢n+1Γ⁡−1−nn≤−200≤n+1−n+12⁢−−1n⁢nΓ⁡−nn≤−100≤n

(7)
> 

simplify⁡zassumingn≤−2

0

(8)
> 

simplify⁡zassuming0≤n

0

(9)
> 

eval⁡z&comma;n=−1

0

(10)

s is an indefinite sum of b:

> 

z≔eval⁡s&comma;n=n+1−s−b

z≔−−1n+1Γ⁡−1−nn≤−200≤n+1−−−1nΓ⁡−nn≤−100≤n−−−1n⁢nΓ⁡−nn≤−100≤n

(11)
> 

simplify⁡zassumingn≤−2

0

(12)
> 

simplify⁡zassuming0≤n

0

(13)
> 

eval⁡z&comma;n=−1

0

(14)

Now assume that 0≤n:

> 

b,s≔BottomSequence⁡T&comma;n&comma;primitiveassuming0≤n

b,s≔n⁢Γ⁡n+1,Γ⁡n+1

(15)

With that assumption, b and T are equivalent, and s is an indefinite sum of both:

> 

simplify⁡b−T

0

(16)
> 

simplify⁡eval⁡s&comma;n=n+1−s−b

0

(17)

Example of a hypergeometric term with parameters:

> 

T≔Γ⁡−nn−k

T≔Γ⁡−nn−k

(18)
> 

BottomSequence⁡T&comma;n

0n≤−1−1n−n+k⁢Γ⁡n+10≤n

(19)

Note that k is considered non-integer.

> 

BottomSequence⁡T&comma;nassumingk::nonnegint

Warning, the assumptions about variable(s) k are ignored

0n≤−1−1n−n+k⁢Γ⁡n+10≤n

(20)
> 

BottomSequence⁡eval⁡T&comma;k=2&comma;n

0n≤1−12n=203≤n

(21)
> 

T≔binomial⁡2⁢n−3&comma;n4n

T≔2⁢n−3n4n

(22)
> 

b,s≔BottomSequence⁡T&comma;n&comma;primitive

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

(23)

References

  

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.

Compatibility

• 

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]

SumTools[Hypergeometric][Gosper]