|
%I #5 Jul 19 2013 05:06:46
%S 1,1,5,61,1481,66361,5390285,803252341,224927827601,121129543555441,
%T 127545238071714965,265238370995975176621,1095520296374502654008921,
%U 9015241470782090221556516521,148067303294213271502974778276445
%N E.g.f. A(x) satisfies: d/dx log(A(x)) = A(2*x)^2.
%C Row 2 of array A159314.
%F E.g.f. satisfies: A'(x) = A(x)*A(2*x)^2.
%F a(n) = Sum_{i=0..n-1} C(n-1,i)*4^i*A126444(i)*a(n-1-i) for n>0 with a(0)=1.
%F E.g.f.: A(x) = G(2*x)^(1/2) where G(x) = e.g.f. of A126444.
%F E.g.f.: A(x) = F(4*x)^(1/4) where F(x) = e.g.f. of A159315.
%e E.g.f.: A(x) = 1 +x +5*x^2/2! +61*x^3/3!+1481*x^4/4!+66361*x^5/5! +...
%e Related expansions:
%e log(A(x)) = x + 4*x^2/2! + 48*x^3/3! + 1216*x^4/4! + 57600*x^5/5! +...
%e A(2*x)^2 = 1 + 4*x + 48*x^2/2! + 1216*x^3/3! + 57600*x^4/4! +...
%e A(x)*A(2*x)^2 = 1 + 5*x +61*x^2/2! +1481*x^3/3! +66361*x^4/4! +...
%o (PARI) {a(n)=local(A=vector(n+4, j, 1+j*x)); for(i=0, n+3, for(j=0, n+2, m=n+3-j; A[m]=exp(intformal((A[m+1]+x*O(x^n))^(2^(m-1)))))); n!*polcoeff(A[3], n, x)}
%Y Cf. A159314, A159315, A126444.
%K nonn
%O 0,3
%A _Paul D. Hanna_, Apr 19 2009
| 