OFFSET
0,2
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 0..200
FORMULA
G.f.: A(x) = (1-(1-16*x)^(1/8))/(2*x).
G.f.: A(x) = (1/x)*(series reversion of x/A097182(x)).
a(n) = A097183(n)/(n+1).
D-finite with recurrence: (n+1)*a(n) +2*(-8*n+1)*a(n-1)=0. - R. J. Mathar, Nov 16 2012
a(n) = 16^n * Gamma(n+7/8) / (Gamma(7/8) * Gamma(n+2)). - Vaclav Kotesovec, Feb 09 2014
a(n) ~ 16^n / (Gamma(7/8) * n^(9/8)). - Vaclav Kotesovec, Feb 09 2014
MAPLE
seq(coeff(series((1-(1-16*x)^(1/8))/(2*x), x, n+2), x, n), n = 0..20); # G. C. Greubel, Sep 17 2019
MATHEMATICA
CoefficientList[Series[(1-(1-16*x)^(1/8))/(2*x), {x, 0, 20}], x] (* Vaclav Kotesovec, Feb 09 2014 *)
Table[FullSimplify[16^n*Gamma[n+7/8]/(Gamma[7/8]*Gamma[n+2])], {n, 0, 20}] (* Vaclav Kotesovec, Feb 09 2014 *)
PROG
(PARI) a(n)=polcoeff((1-(1-16*x+x^2*O(x^n))^(1/8))/(2*x), n, x)
(Magma) R<x>:=PowerSeriesRing(Rationals(), 20); Coefficients(R!( (1-(1-16*x)^(1/8))/(2*x) )); // G. C. Greubel, Sep 17 2019
(Sage)
def A097184_list(prec):
P.<x> = PowerSeriesRing(QQ, prec)
return P((1-(1-16*x)^(1/8))/(2*x)).list()
A097184_list(20) # G. C. Greubel, Sep 17 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Aug 03 2004
EXTENSIONS
More terms from Vincenzo Librandi, Feb 10 2014
STATUS
approved