OFFSET
0,3
COMMENTS
This g.f. arose in a discussion with R. K. Guy about divisibility sequences.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..200
Index entries for linear recurrences with constant coefficients, signature (13,-36,13,-1).
FORMULA
a(0)=0, a(1)=1, a(2)=13, a(3)=133, a(n)=13*a(n-1)-36*a(n-2)+ 13*a(n-3)- a(n-4). - Harvey P. Dale, Mar 10 2015
MATHEMATICA
CoefficientList[Series[x / (x^4 - 13 x^3 + 36 x^2 - 13 x + 1), {x, 0, 40}], x] (* Vincenzo Librandi, Jun 19 2013 *)
LinearRecurrence[{13, -36, 13, -1}, {0, 1, 13, 133}, 30] (* Harvey P. Dale, Mar 10 2015 *)
PROG
(Maxima) makelist(coeff(taylor(x/(x^4-13*x^3+36*x^2-13*x+1), x, 0, n), x, n), n, 0, 22); /* Bruno Berselli, Jun 05 2011 */
(PARI) a=[0, 1, 13, 133]; for(i=1, 99, a=concat(a, 13*a[#a]-36*a[#a-1]+13*a[#a-2]-a[#a-3])); a \\ Charles R Greathouse IV, Jun 05 2011
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Mar 13 2011
STATUS
approved