OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (4,-5,5,-3,1).
FORMULA
G.f.: (1 + x^2)^2/(1 - 4*x + 5*x^2 - 5*x^3 + 3*x^4 - x^5).
MATHEMATICA
<<DiscreteMath`; Do[Print[SeriesTerm[(1 + q^2)^2/(1-4q+5q^2-5q^3+3q^4-q^5), {q, 0, n}]], {n, 0, 50}] (* Ryan Propper, Sep 27 2005 *)
LinearRecurrence[{4, -5, 5, -3, 1}, {1, 4, 13, 37, 101}, 51] (* G. C. Greubel, Apr 19 2021 *)
PROG
(Magma)
R<x>:=PowerSeriesRing(Integers(), 50);
Coefficients(R!( (1+x^2)^2/(1-4*x+5*x^2-5*x^3+3*x^4-x^5) )); // G. C. Greubel, Apr 19 2021
(Sage)
def A097551_list(prec):
P.<x> = PowerSeriesRing(ZZ, prec)
return P( (1+x^2)^2/(1-4*x+5*x^2-5*x^3+3*x^4-x^5) ).list()
A097551_list(50) # G. C. Greubel, Apr 19 2021
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
D n Verma, Aug 16 2004
EXTENSIONS
More terms from Ryan Propper, Sep 27 2005
STATUS
approved