|
|
A041018
|
|
Numerators of continued fraction convergents to sqrt(13).
|
|
14
|
|
|
3, 4, 7, 11, 18, 119, 137, 256, 393, 649, 4287, 4936, 9223, 14159, 23382, 154451, 177833, 332284, 510117, 842401, 5564523, 6406924, 11971447, 18378371, 30349818, 200477279, 230827097, 431304376, 662131473
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,1
|
|
LINKS
|
Index entries for linear recurrences with constant coefficients, signature (0,0,0,0,36,0,0,0,0,1).
|
|
FORMULA
|
(End)
G.f.: (3 + 4*x + 7*x^2 + 11*x^3 + 18*x^4 + 11*x^5 - 7*x^6 + 4*x^7 - 3*x^8 + x^9)/(1 - 36*x^5 - x^10). - Peter J. C. Moses, Jul 29 2013
|
|
MAPLE
|
a[0]:=3: a[-1]:=1: b(0):=6: b(1):=1; b(2):=1: b(3):=1: b(4):=1:
for n from 1 to 100 do k:=n mod 5:
a[n]:=b(k)*a[n-1]+a[n-2]:
printf("%12d", a[n]):
|
|
MATHEMATICA
|
CoefficientList[Series[(3 + 4*x + 7*x^2 + 11*x^3 + 18*x^4 + 11*x^5 - 7*x^6 + 4*x^7 - 3*x^8 + x^9)/(1 - 36*x^5 - x^10), {x, 0, 50}], x] (* Stefano Spezia, Aug 31 2018 *)
|
|
CROSSREFS
|
Cf. A010122 (continued fraction for sqrt(13)).
|
|
KEYWORD
|
nonn,frac,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|