login
A200949
E.g.f. A(x) satisfies A''(x)+A'(x)=exp(A(x))+1.
0
1, 1, 0, 2, 2, 10, 28, 130, 564, 3102, 17464, 114282, 784008, 5987482, 48279164, 422742854, 3909897740, 38675519334, 403121938616, 4449327971330, 51570742407808, 628575203791026, 8017017724008788, 106974723312803902, 1488628872502049348
OFFSET
1,4
PROG
(Maxima)
Composita(n, k, F):=block([], if k=1 then F(n) else if n=k then (F(1))^n else sum(F(i+1)*Composita(n-i-1, k-1, F), i, 0, n-k));
SolveDiff21(G, key):=block([b], array(b, 100), b[1]:1, b[2]:1/2, bfun(n):=b[n], makelist(b[n]:-b[n-1]/n+1/(n*(n-1))*sum(Composita(n-2, k, bfun)*G(k) , k, 1, n-2), n, 3, 20),
if key=1 then makelist(n!*b[n], n, 1, 20) else makelist(b[n], n, 1, 20));
f(n):=1/n!;
SolveDiff21(f, 1);
CROSSREFS
Sequence in context: A309753 A102345 A151364 * A001885 A300641 A078433
KEYWORD
nonn
AUTHOR
Vladimir Kruchinin, Nov 25 2011
STATUS
approved