OFFSET
0,2
COMMENTS
Value of a(5) is incorrectly given as 66693 in Henrici paper. - Sean A. Irvine, Jun 20 2013
REFERENCES
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
Peter Henrici, Automatic computations with power series, J. Assoc. Comput. Mach. 3 (1956), 10-15.
FORMULA
Let f(x) = [Sum_{k>=1}(3/(2*k+1)) * x^(2*k+1)]^(1/3) = x + (1/5)*x^3 + (18/175) * x^5 + ...; let g(x) be the Lagrange inversion of f(x), g(x) = REVERT(f(x)) = 1 - (1/5) * x^3 + (3/175) * x^5 + .... Then a(n) = denominator((2 * n + 1) * coeff(g(x), 2*n+1)). - Sean A. Irvine, Jun 20 2013
MATHEMATICA
nmax = 14;
S = Sum[(3/(2k+1)) x^(2k+1), {k, 1, Infinity}]^(1/3) + O[x]^(3nmax) // Normal // Simplify[#, x > 0]& // InverseSeries[# + O[x]^(3nmax), x]&;
a[n_] := Denominator[(2n+1) SeriesCoefficient[S, {x, 0, 2n+1}]];
a /@ Range[0, nmax] (* Jean-François Alcover, Oct 01 2020 *)
CROSSREFS
KEYWORD
nonn,frac
AUTHOR
EXTENSIONS
More terms from Sean A. Irvine, Jun 20 2013
STATUS
approved