cp's OEIS Frontend

A289739 Expansion of solution to dy/dx = y + exp(y).

%I A289739 #18 Aug 21 2017 13:21:36
%S A289739 0,1,2,5,17,79,474,3468,29799,293528,3258373,40234231,546921835,
%T A289739 8115147998,130503876054,2260929219675,41979302557200,831593152814251,
%U A289739 17506400133530765,390278100156698627,9185223726173708408,227578002295869672508,5921091852493279814589
%N A289739 Expansion of solution to dy/dx = y + exp(y).
%C A289739 See A266539 for more details.
%H A289739 Robert Israel, <a href="/A289739/b289739.txt">Table of n, a(n) for n = 0..435</a>
%F A289739 E.g.f. y(x) = log(A(x)) and y'(x) = B(x) where A(x), B(x) are as in A266539.
%F A289739 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
%e A289739 E.g.f. = x + 2*x^2/2! + 5*x^3/3! + 17*x^4/4! + ...
%p A289739 S:= dsolve({diff(y(x),x) = y(x) + exp(y(x)), y(0)=0},y(x),series,order=31):
%p A289739 seq(coeff(rhs(S),x,j)*j!,j=0..30); # _Robert Israel_, Aug 09 2017
%t A289739 a[ n_] := If[ n < 0, 0, n! SeriesCoefficient[ InverseSeries[ Series[Integrate[ 1 / (x + Exp[x]), x], {x, 0, n}]], {x, 0, n}]];
%o A289739 (PARI) {a(n) = if( n<0, 0, my(A = O(x)); for(k=1, n, A = intformal(A + exp(A))); n! * polcoeff(A, n))};
%o A289739 (PARI) {a(n) = if( n<0, 0, n! * polcoeff( serreverse( intformal( 1 / (exp(x + x * O(x^n)) + x))), n))};
%Y A289739 Cf. A266329.
%K A289739 nonn
%O A289739 0,3
%A A289739 _Michael Somos_, Aug 09 2017