cp's OEIS Frontend

A228173 E.g.f. satisfies: A(x) = 1+x - exp(-A(x)^2).

%I A228173 #13 Jan 11 2014 10:43:32
%S A228173 1,2,12,108,1320,20280,374640,8072400,198465120,5475284640,
%T A228173 167285321280,5600184004800,203602252613760,7978382871338880,
%U A228173 334767145102790400,14952953514231532800,707221717016278464000,35242469168705967168000,1841491290250262851200000
%N A228173 E.g.f. satisfies: A(x) = 1+x - exp(-A(x)^2).
%C A228173 Note that a(30) is negative. - _Vaclav Kotesovec_, Sep 16 2013
%H A228173 Vaclav Kotesovec, <a href="/A228173/b228173.txt">Table of n, a(n) for n = 1..370</a>
%H A228173 Vaclav Kotesovec, <a href="/A228173/a228173.jpg">Graph of convergence to limit for 1000 terms</a>
%F A228173 E.g.f. A(x) satisfies:
%F A228173 (1) A(x - 1 + exp(-x^2)) = x.
%F A228173 (2) A(x) = x + Sum_{n>=1} d^(n-1)/dx^(n-1) (1 - exp(-x^2))^n / n!.
%F A228173 (3) A(x) = x*exp( Sum_{n>=1} d^(n-1)/dx^(n-1) (1 - exp(-x^2))^n/x / n! ).
%F A228173 Lim sup n->infinity (|a(n)|/n!)^(1/n) = 1/abs(-1-(LambertW(-1/2)-1) / sqrt(-2*LambertW(-1/2))) = 3.19002880735268... - _Vaclav Kotesovec_, Jan 11 2014
%e A228173 E.g.f.: A(x) = x + 2*x^2/2! + 12*x^3/3! + 108*x^4/4! + 1320*x^5/5! +...
%e A228173 where
%e A228173 exp(-A(x)^2) = 1 - 2*x^2/2! - 12*x^3/3! - 108*x^4/4! - 1320*x^5/5! -...
%e A228173 The e.g.f. equals the series:
%e A228173 A(x) = x + (1 - exp(-x^2)) + d/dx (1 - exp(-x^2))^2/2! + d^2/dx^2 (1 - exp(-x^2))^3/3! + d^3/dx^3 (1 - exp(-x^2))^4/4! + d^4/dx^4 (1 - exp(-x^2))^5/5! +...
%e A228173 Also,
%e A228173 log(A(x)/x) = (1 - exp(-x^2))/x + d/dx (1 - exp(-x^2))^2/(2!*x) + d^2/dx^2 (1 - exp(-x^2))^3/(3!*x) + d^3/dx^3 (1 - exp(-x^2))^4/(4!*x) +...
%t A228173 Rest[CoefficientList[InverseSeries[Series[x-1+E^(-x^2),{x,0,20}],x],x]*Range[0,20]!] (* _Vaclav Kotesovec_, Sep 16 2013 *)
%o A228173 (PARI) {a(n)=n!*polcoeff(serreverse(x-1+exp(-x^2+x*O(x^n))), n)}
%o A228173 for(n=1, 25, print1(a(n), ", "))
%o A228173 (PARI) {Dx(n, F)=local(D=F); for(i=1, n, D=deriv(D)); D}
%o A228173 {a(n)=local(A=x); A=x+sum(m=1, n, Dx(m-1, (1 - exp(-x^2+x*O(x^n)))^m)/m!); n!*polcoeff(A, n)}
%o A228173 for(n=1, 25, print1(a(n), ", "))
%o A228173 (PARI) {Dx(n, F)=local(D=F); for(i=1, n, D=deriv(D)); D}
%o A228173 {a(n)=local(A=x+x^2+x*O(x^n)); A=x*exp(sum(m=1, n, Dx(m-1, (1 - exp(-x^2+x*O(x^n)))^m/x)/m!)+x*O(x^n)); n!*polcoeff(A, n)}
%o A228173 for(n=1, 25, print1(a(n), ", "))
%Y A228173 Cf. A218652.
%K A228173 sign
%O A228173 1,2
%A A228173 _Paul D. Hanna_, Aug 14 2013