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”).

A042744
Numerators of continued fraction convergents to sqrt(903).
2
30, 601, 36090, 722401, 43380150, 868325401, 52142904210, 1043726409601, 62675727480270, 1254558276015001, 75336172288380330, 1507978004043621601, 90554016414905676390, 1812588306302157149401, 108845852394544334640450
OFFSET
0,1
FORMULA
G.f.: (30 +601*x +30*x^2 -x^3)/(1 -1202*x^2 +x^4). - Vincenzo Librandi, Dec 03 2013
a(n) = 1202*a(n-2) - a(n-4). - Vincenzo Librandi, Dec 03 2013
MATHEMATICA
Numerator[Convergents[Sqrt[903], 30]] (* or *) CoefficientList[Series[(30 + 601 x + 30 x^2 - x^3)/(1 - 1202 x^2 + x^4), {x, 0, 30}], x] (* Vincenzo Librandi, Dec 03 2013 *)
LinearRecurrence[{0, 1202, 0, -1}, {30, 601, 36090, 722401}, 20] (* Harvey P. Dale, May 25 2022 *)
PROG
(Magma) I:=[30, 601, 36090, 722401]; [n le 4 select I[n] else 1202*Self(n-2)-Self(n-4): n in [1..30]]; // Vincenzo Librandi, Dec 03 2013
CROSSREFS
Cf. A042745.
Sequence in context: A131206 A020982 A024436 * A020980 A051303 A020975
KEYWORD
nonn,cofr,frac,easy
AUTHOR
STATUS
approved