|
MATHEMATICA
|
a[n_] := a[n] = 1 + 4 Sum[Binomial[n, k] a[k] a[n - k - 1], {k, 0, n - 1}]; Table[a[n], {n, 0, 17}]
nmax = 17; CoefficientList[Series[Exp[x]/Sqrt[9 - 8 Exp[x]], {x, 0, nmax}], x] Range[0, nmax]!
Table[Sum[Sum[Binomial[n, k] StirlingS2[k, j] 4^j (2 j - 1)!!, {j, 0, k}], {k, 0, n}], {n, 0, 17}]
|