gfun
poltorec
determine the recurrence satisfied by a polynomial in holonomic sequences
Calling Sequence
Parameters
Description
Examples
poltorec(P, listrec, list_unknowns, u(n))
P
-
polynomial in n and u1(n), u2(n), ... and possibly their shifts (u1(n+1), u2(n+1), ...) and repeated shifts
listrec
list containing, for each of u1(n), u2(n), ..., either a linear recurrence equation it satisfies or a set containing the equation together with initial conditions
list_unknowns
list of sequences [u1n,u2n,...]
u
name; holonomic sequence name
n
name; variable of the holonomic sequence u
The poltorec(P, listrec, list_unknowns, u(n)) command returns the recurrence satisfied by the polynomial P.
If u1n, u2n, ... are holonomic sequence solutions of listrec[1], listrec[2], ..., the poltorec function returns a linear recurrence equation satisfied by Pn,u1n,....
withgfun:
rec1≔u10=1,u1n+1=n+1u1n:
rec2≔u20=1,u21=1,u2n+2=2u2n+1−3nu2n:
poltorecu1n2+2u1nu2n,rec1,rec2,u1n,u2n,un
−3n7−39n6−192n5−462n4−579n3−363n2−90nun+5n5+54n4+209n3+354n2+254n+60un+1+−n4−12n3−46n2−62n−15un+2+n2+4nun+3,u0=3,u1=3,u2=12,u3=48
Cassini's identity:
fib≔F0=1,F1=1,Fn+2=Fn+1+Fn:
poltorecFn+2Fn−Fn+12,fib,Fn,fn
fn+1+fn,f0=1
See Also
gfun['rec+rec']
gfun[`rec*rec`]
gfun[parameters]
gfun[poltodiffeq]
Download Help Document