[go: up one dir, main page]

login
A042208
Numerators of continued fraction convergents to sqrt(630).
2
25, 251, 12575, 126001, 6312625, 63252251, 3168925175, 31752504001, 1590794125225, 15939693756251, 798575481937775, 8001694513134001, 400883301138637825, 4016834705899512251, 201242618596114250375
OFFSET
0,1
FORMULA
G.f.: (25 +251*x +25*x^2 -x^3)/(x^4 -502*x^2 +1). - Vincenzo Librandi, Nov 19 2013
a(n) = 502*a(n-2) - a(n-4). - Vincenzo Librandi, Nov 19 2013
MATHEMATICA
Numerator[Convergents[Sqrt[630], 30]] (* or *) CoefficientList[Series[(25 + 251 x + 25 x^2 - x^3)/(x^4 - 502 x^2 + 1), {x, 0, 20}], x] (* Vincenzo Librandi, Nov 19 2013 *)
LinearRecurrence[{0, 502, 0, -1}, {25, 251, 12575, 126001}, 20] (* Harvey P. Dale, Apr 01 2018 *)
PROG
(Magma) I:=[25, 251, 12575, 126001]; [n le 4 select I[n] else 502*Self(n-2)-Self(n-4): n in [1..30]]; // Vincenzo Librandi, Nov 19 2013
CROSSREFS
Cf. A042209.
Sequence in context: A023073 A022685 A308492 * A143009 A090022 A298068
KEYWORD
nonn,cofr,frac,easy
AUTHOR
STATUS
approved