OFFSET
0,3
COMMENTS
a(2) = 268 explains the comment in A021085 that "The decimal expansion of Sum_{n>=1} floor(n * tanh(Pi))/10^n is the same as that of 1/81 for the first 268 decimal places [Borwein et al.]".
REFERENCES
J. M. Borwein, D. H. Bailey and R. Girgensohn, Experimentation in Mathematics, A K Peters, Ltd., Natick, MA, 2004. x+357 pp. See p. 13.
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1920
FORMULA
a(n) = A060402(n)*a(n-1) + a(n-2) for n >= 2. - Nathaniel Johnston, May 10 2011
MAPLE
lim:=50: with(numtheory): cfr := cfrac(tanh(Pi), lim+10, 'quotients'): q[0]:=1:q[1]:=cfr[2]: printf("%d, %d, ", q[0], q[1]): for n from 2 to lim do q[n]:=cfr[n+1]*q[n-1]+q[n-2]: printf("%d, ", q[n]): od: # Nathaniel Johnston, May 10 2011
MATHEMATICA
a[0] := 1; a[1] := 1; A060402:= ContinuedFraction[Tanh[Pi], 100];
a[n_]:= a[n] = A060402[[n + 1]]*a[n - 1] + a[n - 2]; Join[{1, 1}, Table[a[n], {n, 2, 75}]] (* G. C. Greubel, Apr 05 2018 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
N. J. A. Sloane, May 09 2011
EXTENSIONS
a(4)-a(24) from Nathaniel Johnston, May 10 2011
STATUS
approved