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

A254795
Numerators of the convergents of the generalized continued fraction 2 + 1^2/(4 + 3^2/(4 + 5^2/(4 + ... ))).
3
2, 9, 54, 441, 4410, 53361, 747054, 12006225, 216112050, 4334247225, 95353438950, 2292816782025, 59613236332650, 1671463434096225, 50143903022886750, 1606276360166472225, 54613396245660055650, 1967688541203928475625, 74772164565749282073750
OFFSET
0,1
COMMENTS
Brouncker gave the generalized continued fraction expansion 4/Pi = 1 + 1^2/(2 + 3^2/(2 + 5^2/(2 + ... ))). The sequence of convergents begins [1/1, 3/2, 15/13, 105/76, ... ]. The numerators of the convergents are in A001147, the denominators in A024199.
In extending Brouckner's result, Osler showed that 2 + 1^2/(4 + 3^2/(4 + 5^2/(4 + ... ))) = L^2/Pi, where L is the lemniscate constant A062539. The sequence of convergents to Osler's continued fraction begins [2/1, 9/4, 54/25, 441/200, 4410/2025, ...]. Here we list the (unreduced) numerators of these convergents. See A254796 for the sequence of denominators. See A254794 for the decimal expansion of L^2/Pi.
FORMULA
a(2*n-1) = ( A008545(n) )^2 = ( Product {k = 0..n-1} 4*k + 3 )^2.
a(2*n) = (4*n + 2)*( A008545(n) )^2 = (4*n + 2)*( Product {k = 0..n-1} 4*k + 3 )^2.
a(n) = 4*a(n-1) + (2*n - 1)^2*a(n-2) with a(0) = 2, a(1) = 9.
a(2*n) = (4*n + 2)*a(2*n-1); a(2*n+1) = (4*n + 4)*a(2*n) + a(2*n-1).
MAPLE
a[0] := 2: a[1] := 9:
for n from 2 to 18 do a[n] := 4*a[n-1] + (2*n-1)^2*a[n-2] end do:
seq(a[n], n = 0 .. 18);
CROSSREFS
KEYWORD
nonn,frac,easy
AUTHOR
Peter Bala, Feb 23 2015
STATUS
approved