login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A041920
Numerators of continued fraction convergents to sqrt(482).
2
21, 22, 461, 483, 20747, 21230, 445347, 466577, 20041581, 20508158, 430204741, 450712899, 19360146499, 19810859398, 415577334459, 435388193857, 18701881476453, 19137269670310, 401447274882653, 420584544552963, 18065998146107099, 18486582690660062
OFFSET
0,1
FORMULA
G.f.: (21 +22*x +461*x^2 +483*x^3 +461*x^4 -22*x^5 +21*x^6 -x^7)/(1 -966*x^4 +x^8). - Vincenzo Librandi, Nov 12 2013
a(n) = 966*a(n-4) - a(n-8). - Vincenzo Librandi, Nov 12 2013
MATHEMATICA
Numerator[Convergents[Sqrt[482], 30]] (* or *) CoefficientList[Series[(21 + 22 x + 461 x^2 + 483 x^3 + 461 x^4 - 22 x^5 + 21 x^6 - x^7)/(1 - 966 x^4 + x^8), {x, 0, 30}], x] (* Vincenzo Librandi, Nov 12 2013 *)
LinearRecurrence[{0, 0, 0, 966, 0, 0, 0, -1}, {21, 22, 461, 483, 20747, 21230, 445347, 466577}, 30] (* Harvey P. Dale, May 28 2023 *)
PROG
(Magma) I:=[21, 22, 461, 483, 20747, 21230, 445347, 466577]; [n le 8 select I[n] else 966*Self(n-4)-Self(n-8): n in [1..30]]; // Vincenzo Librandi, Nov 12 2013
CROSSREFS
Sequence in context: A041918 A068325 A041219 * A041921 A041922 A104669
KEYWORD
nonn,frac,easy,less
AUTHOR
EXTENSIONS
More terms and second link corrected by Colin Barker, Dec 28 2013
STATUS
approved