|
FORMULA
|
E.g.f. A(x) satisfies:
(1) A(x) = x + A(x)^2*exp(A(x)^2).
(2) A(x) = x*Catalan( x*exp(A(x)^2) ) where Catalan(x) = (1-sqrt(1-4*x))/(2*x).
(3) A(x) = x*Sum_{n>=0} binomial(2*n+1,n)/(2*n+1) * x^n * exp(n*A(x)^2).
(4) A(x) = x*exp( Sum_{n>=1} binomial(2*n-1,n) * x^n/n * exp(n*A(x)^2) ).
(5) A(x) = x + Sum_{n>=1} d^(n-1)/dx^(n-1) x^(2*n) * exp(n*x^2) / n!.
(6) A(x) = x*exp( Sum_{n>=1} d^(n-1)/dx^(n-1) x^(2*n-1) * exp(n*x^2) / n! ).
a(n) ~ n^(n-1) * s * sqrt((1+2*s^2)/(2+10*s^2+4*s^4)) / (exp(n) * ((s*(1+2*s^2))/(2*(1+s^2)))^n), where s = 0.3788063540000847107637564... is the root of the equation 2*s*(1+s^2)*exp(s^2) = 1. - Vaclav Kotesovec, Jan 07 2014
|
|
PROG
|
(PARI) {a(n)=n!*polcoeff(serreverse(x-x^2*exp(x^2 +x*O(x^n))), n)}
for(n=1, 25, print1(a(n), ", "))
(PARI) /* E.g.f. A(x) = x*Catalan( x*exp(A(x)^2) ): */
{a(n)=local(A=x); for(i=1, n, A=(1-sqrt(1-4*x*exp(A^2 +x^2*O(x^n)) ))/2*exp(-A^2 +x*O(x^n)) ); n!*polcoeff(A, n)}
for(n=1, 25, print1(a(n), ", "))
(PARI) {Dx(n, F)=local(D=F); for(i=1, n, D=deriv(D)); D}
{a(n)=local(A=x); A=x+sum(m=1, n, Dx(m-1, x^(2*m)*exp(x^2+x*O(x^n))^m/m!)); n!*polcoeff(A, n)}
for(n=1, 25, print1(a(n), ", "))
(PARI) {Dx(n, F)=local(D=F); for(i=1, n, D=deriv(D)); D}
{a(n)=local(A=x+x^2+x*O(x^n)); A=x*exp(sum(m=1, n, Dx(m-1, x^(2*m-1)*exp(x^2+x*O(x^n))^m/m!)+x*O(x^n))); n!*polcoeff(A, n)}
for(n=1, 25, print1(a(n), ", "))
|