OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..50
FORMULA
a(n) = [y^n] ( sqrt(1+y^2) + y )^(2^n), since log(sqrt(1+y^2) + y) = arcsinh(y); [y^n] F(y) denotes the coefficient of y^n in F(y).
EXAMPLE
G.f.: A(x) = 1 + 2*x + 8*x^2 + 84*x^3 + 2688*x^4 + 276892*x^5 + ...
This is a special application of the following identity.
Let F(x), G(x), be power series in x such that F(0)=1, G(0)=1, then
Sum_{n>=0} m^n * H(q^n*x) * log( F(q^n*x)*G(x) )^n / n! =
Sum_{n>=0} x^n * G(x)^(m*q^n) * [y^n] H(y)*F(y)^(m*q^n).
MAPLE
m:=30;
S:=series( add( arcsinh(2^j*x)^j/j! , j=0..m+2), x, m+1):
seq(coeff(S, x, j), j=0..m); # G. C. Greubel, Mar 15 2021
MATHEMATICA
With[{m=30}, CoefficientList[Series[Sum[ArcSinh[2^j*x]^j/j!, {j, 0, m+2}], {x, 0, m}], x]] (* G. C. Greubel, Mar 15 2021 *)
PROG
(PARI) {a(n)=polcoeff(sum(k=0, n, asinh(2^k*x +x*O(x^n))^k/k!), n)}
(PARI) {a(n)=polcoeff((sqrt(1+x^2)+x+x*O(x^n))^(2^n), n)}
(Magma)
m:=30;
R<x>:=PowerSeriesRing(Rationals(), 30);
Coefficients(R!( (&+[Argsinh(2^j*x)^j/Factorial(j): j in [0..m+2]]) )); // G. C. Greubel, Mar 15 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jan 20 2008
STATUS
approved