cp's OEIS Frontend

A224788 E.g.f. satisfies: A(x) = exp( Integral A(x)/(1 - x*A(x)^2) dx ).

%I A224788 #23 Feb 16 2025 08:33:19
%S A224788 1,1,3,18,168,2142,34704,682740,15810372,421339176,12702393792,
%T A224788 427435993512,15881634963216,645804320863680,28527455317884336,
%U A224788 1360332028008819360,69645942884911181184,3810436222004101378656,221867131720533800409216,13698420738298341356760768
%N A224788 E.g.f. satisfies: A(x) = exp( Integral A(x)/(1 - x*A(x)^2) dx ).
%C A224788 Compare to: C(x) = exp( Integral C(x)^2/(1 - x*C(x)^2) dx ), which is satisfied by: C(x) = (1-sqrt(1-4*x))/(2*x) (Catalan numbers, A000108).
%C A224788 Compare to: W(x) = exp( Integral W(x)/(1 - x*W(x)) dx ), which is satisfied by: W(x) = LambertW(-x)/(-x)  =  Sum_{n>=0} (n+1)^(n-1)*x^n/n!.
%H A224788 Vaclav Kotesovec, <a href="/A224788/b224788.txt">Table of n, a(n) for n = 0..300</a>
%H A224788 Eric Weisstein, <a href="https://mathworld.wolfram.com/ExponentialIntegral.html">MathWorld: Exponential Integral</a>
%F A224788 E.g.f. derivative: A'(x) = A(x)^2 / (1-x*A(x)^2). - _Vaclav Kotesovec_, Feb 19 2014
%F A224788 a(n) ~ n^(n-1) / (sqrt(2) * exp(n) * r^(n+1/4)), where r = 0.28396034297... is the root of the equation Ei(1/sqrt(r)) - Ei(1) + exp(1) = (r+sqrt(r)) * exp(1/sqrt(r)), where Ei is the Exponential Integral. - _Vaclav Kotesovec_, Feb 19 2014
%e A224788 E.g.f.: A(x) = 1 + x + 3*x^2/2! + 18*x^3/3! + 168*x^4/4! + 2142*x^5/5! +...
%e A224788 where
%e A224788 log(A(x)) = x + 2*x^2/2! + 11*x^3/3! + 99*x^4/4! + 1236*x^5/5! + 19752*x^6/6! +...
%e A224788 A(x)/(1-x*A(x)^2) = 1 + 2*x + 11*x^2/2! + 99*x^3/3! + 1236*x^4/4! + 19752*x^5/5! +...
%t A224788 a = ConstantArray[0,21]; a[[1]]=1; a[[2]]=1; Do[a[[n+2]] = n!*Sum[a[[i+1]]*a[[n-i+1]]/i!/(n-i)!,{i,0,n}] + n!*Sum[a[[j+1]]/(j-1)!*Sum[a[[i+1]]*a[[n-j-i+1]]/i!/(n-j-i)!,{i,0,n}],{j,1,n}],{n,1,18}]; a (* _Vaclav Kotesovec_, Feb 19 2014 *)
%t A224788 FindRoot[ExpIntegralEi[1/Sqrt[r]] - ExpIntegralEi[1] + E == (r+Sqrt[r]) * E^(1/Sqrt[r]),{r,1/2},WorkingPrecision->50] (* program for numerical value of the radius of convergence r, _Vaclav Kotesovec_, Feb 19 2014 *)
%o A224788 (PARI) {a(n)=local(A=1+x);for(i=1,n,A=exp(intformal(A/(1-x*A^2 +x*O(x^n)))));n!*polcoeff(A,n)}
%o A224788 for(n=0,20,print1(a(n),", "))
%Y A224788 Cf. A225052, A091725.
%K A224788 nonn
%O A224788 0,3
%A A224788 _Paul D. Hanna_, Apr 28 2013