login
A042528
Numerators of continued fraction convergents to sqrt(793).
2
28, 169, 704, 4393, 246712, 1484665, 6185372, 38596897, 2167611604, 13044266521, 54344677688, 339112332649, 19044635306032, 114606924168841, 477472331981396, 2979440916057217, 167326163631185548, 1006936422703170505
OFFSET
0,1
FORMULA
G.f.: (28 +169*x +704*x^2 +4393*x^3 +704*x^4 -169*x^5 +28*x^6 -x^7)/(1 -8786*x^4 +x^8). - Vincenzo Librandi, Nov 26 2013
a(n) = 8786*a(n-4) - a(n-8). - Vincenzo Librandi, Nov 26 2013
MATHEMATICA
Numerator[Convergents[Sqrt[793], 30]] (* or *) CoefficientList[Series[(28 + 169 x + 704 x^2 + 4393 x^3 + 704 x^4 - 169 x^5 + 28 x^6 - x^7)/(1 - 8786 x^4 + x^8), {x, 0, 30}], x] (* Vincenzo Librandi, Nov 25 2013 *)
PROG
(Magma) I:=[28, 169, 704, 4393, 246712, 1484665, 6185372, 38596897]; [n le 8 select I[n] else 8786*Self(n-4)-Self(n-8): n in [1..30]]; // Vincenzo Librandi, Nov 26 2013
CROSSREFS
Cf. A042529.
Sequence in context: A305270 A366954 A129136 * A345215 A219821 A341160
KEYWORD
nonn,cofr,frac,easy
AUTHOR
STATUS
approved