OFFSET
0,2
COMMENTS
Let ((1 + k*x)/(1 - k*x))^(1/k) = a(0) + a(1)*x + a(2)*x^2 + ...
Then n*a(n) = 2*a(n-1) + k^2*(n-2)*a(n-2) for n > 1.
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..1000
FORMULA
a(n) ~ 2^(2*n + 1/4) / (Gamma(1/4) * n^(3/4)). - Vaclav Kotesovec, Apr 26 2018
n*a(n) = 2*a(n-1) + 4^2*(n-2)*a(n-2) for n > 1.
G.f.: A(x)=F(x*G(x^2)), where F(x) is the g.f. for A063886, and G(x) is the g.f. for A151403. - Alexander Burstein, Nov 13 2023
MATHEMATICA
CoefficientList[Series[Surd[(1+4x)/(1-4x), 4], {x, 0, 40}], x] (* Harvey P. Dale, Jul 25 2021 *)
PROG
(PARI) N=66; x='x+O('x^N); Vec(((1+4*x)/(1-4*x))^(1/4))
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Apr 25 2018
STATUS
approved