OFFSET
0,3
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..200
Index entries for linear recurrences with constant coefficients, signature (0, 0, 0, 0, 0, 0, 0, 3386, 0, 0, 0, 0, 0, 0, 0, -1).
FORMULA
G.f.: -(x^2-x-1)*(x^4+3*x^2+1)*(x^8+22*x^4+1) / (x^16-3386*x^8+1). - Colin Barker, Dec 06 2013
a(n) = 3386*a(n-8) - a(n-16) for n>15. - Vincenzo Librandi, Jan 19 2014
MAPLE
convert(sqrt(658), confrac, 30, cvgts): denom(cvgts); # Wesley Ivan Hurt, Dec 07 2013
MATHEMATICA
Denominator[Convergents[Sqrt[658], 30]] (* Vincenzo Librandi, Jan 19 2014 *)
PROG
(Magma) I:=[1, 1, 2, 3, 20, 23, 43, 66, 3343, 3409, 6752, 10161, 67718, 77879, 145597, 223476]; [n le 16 select I[n] else 3386*Self(n-8)-Self(n-16): n in [1..40]]; // Vincenzo Librandi, Jan 19 2014
CROSSREFS
KEYWORD
nonn,frac,easy
AUTHOR
EXTENSIONS
More terms from Colin Barker, Dec 06 2013
STATUS
approved