OFFSET
1,1
COMMENTS
This equation is used for worked examples in the Robertson paper.
(1/8)*a(n)^2 -5 is a term of A152741. [Bruno Berselli, Aug 17 2013]
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
John P. Robertson, Solving the generalized Pell equation x^2 - Dy^2 = N
Index entries for linear recurrences with constant coefficients, signature (0,0,0,1298,0,0,0,-1).
FORMULA
G.f.: -4*x*(x-1)*(3*x^6+13*x^5+68*x^4+260*x^3+68*x^2+13*x+3) / ((x^4-36*x^2-1)*(x^4+36*x^2-1)).
a(n) = 1298*a(n-4)-a(n-8).
MATHEMATICA
CoefficientList[Series[-4 (x - 1) (3 x^6 + 13 x^5 + 68 x^4 + 260 x^3 + 68 x^2 + 13 x+3) / ((x^4 - 36 x^2 - 1) (x^4 + 36 x^2 - 1)), {x, 0, 30}], x] (* Vincenzo Librandi, Aug 17 2013 *)
LinearRecurrence[{0, 0, 0, 1298, 0, 0, 0, -1}, {12, 40, 220, 768, 14808, 51700, 285520, 996852}, 30] (* Harvey P. Dale, Apr 22 2024 *)
PROG
(PARI) Vec(-4*x*(x-1)*(3*x^6+13*x^5+68*x^4+260*x^3+68*x^2+13*x+3)/((x^4-36*x^2-1)*(x^4+36*x^2-1)) + O(x^100))
(Magma) I:=[12, 40, 220, 768, 14808, 51700, 285520, 996852]; [n le 8 select I[n] else 1298*Self(n-4)-Self(n-8): n in [1..30]]; // Vincenzo Librandi, Aug 17 2013
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Colin Barker, Aug 16 2013
STATUS
approved