The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A289739 Expansion of solution to dy/dx = y + exp(y). 2
 0, 1, 2, 5, 17, 79, 474, 3468, 29799, 293528, 3258373, 40234231, 546921835, 8115147998, 130503876054, 2260929219675, 41979302557200, 831593152814251, 17506400133530765, 390278100156698627, 9185223726173708408, 227578002295869672508, 5921091852493279814589 (list; graph; refs; listen; history; text; internal format)
 OFFSET 0,3 COMMENTS See A266539 for more details. LINKS Robert Israel, Table of n, a(n) for n = 0..435 FORMULA E.g.f. y(x) = log(A(x)) and y'(x) = B(x) where A(x), B(x) are as in A266539. a(n) ~ c^n * (n-1)!, where c = 1/Integral_{x=0..infinity} 1/(x + exp(x)) dx = 1.2400861064984976662394901721056528110217273471501174317019052800276... - Vaclav Kotesovec, Aug 21 2017 EXAMPLE E.g.f. = x + 2*x^2/2! + 5*x^3/3! + 17*x^4/4! + ... MAPLE S:= dsolve({diff(y(x), x) = y(x) + exp(y(x)), y(0)=0}, y(x), series, order=31): seq(coeff(rhs(S), x, j)*j!, j=0..30); # Robert Israel, Aug 09 2017 MATHEMATICA a[ n_] := If[ n < 0, 0, n! SeriesCoefficient[ InverseSeries[ Series[Integrate[ 1 / (x + Exp[x]), x], {x, 0, n}]], {x, 0, n}]]; PROG (PARI) {a(n) = if( n<0, 0, my(A = O(x)); for(k=1, n, A = intformal(A + exp(A))); n! * polcoeff(A, n))}; (PARI) {a(n) = if( n<0, 0, n! * polcoeff( serreverse( intformal( 1 / (exp(x + x * O(x^n)) + x))), n))}; CROSSREFS Cf. A266329. Sequence in context: A187245 A361243 A302194 * A243337 A259622 A054499 Adjacent sequences: A289736 A289737 A289738 * A289740 A289741 A289742 KEYWORD nonn AUTHOR Michael Somos, Aug 09 2017 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified April 18 15:35 EDT 2024. Contains 371780 sequences. (Running on oeis4.)