OFFSET
0,2
FORMULA
b(0)=1. b(n) = -sum{k=1 to n} binomial(n,k) H(k+1) b(n-k)/(k+1).
EXAMPLE
1/(1 + x * 3/(2 * 2) + x^2 * 11/(6 * 6) + x^3 * 25/(12 * 24) +...) = 1 -x * 3/4 + x^2 * 37/72 -x^3 * 29/96 ...
MATHEMATICA
b[0] = 1; b[n_] := b[n] = -Sum[Binomial[n, k] *HarmonicNumber[k + 1]*b[n - k]/(k + 1), {k, n}]; Denominator[Array[b, 20, 0]] (* Ray Chandler, Feb 19 2007 *)
CROSSREFS
KEYWORD
frac,nonn
AUTHOR
Leroy Quet, Feb 13 2007
EXTENSIONS
Extended by Ray Chandler, Feb 19 2007
STATUS
approved