OFFSET
1,2
FORMULA
E.g.f. A(x) satisfies: A(x) = asin(x*exp(A(x))).
a(n) ~ 2^((n-1)/2) * exp(n*Pi/4) * n^(n-1) / exp(n). - Vaclav Kotesovec, Jan 10 2014
EXAMPLE
E.g.f.: A(x) = x + 2*x^2/2! + 10*x^3/3! + 80*x^4/4! + 884*x^5/5! + 12480*x^6/6! +...
where A( sin(x)/exp(x) ) = x.
MATHEMATICA
Rest[CoefficientList[InverseSeries[Series[Sin[x]/E^x, {x, 0, 20}], x], x] * Range[0, 20]!] (* Vaclav Kotesovec, Jan 10 2014 *)
PROG
(PARI) {a(n)=local(X=x+x*O(x^n)); n!*polcoeff(serreverse(sin(X)/exp(X)), n)}
for(n=1, 25, print1(a(n), ", "))
(PARI) {a(n)=local(A=x); for(i=1, n, A=asin(x*exp(A+x*O(x^n)))); n!*polcoeff(A, n)}
for(n=1, 25, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jul 14 2013
STATUS
approved