OFFSET
0,3
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 0..402
FORMULA
E.g.f.: arcsin(2*x)/(2*(1-2*x)^(3/2)).
a(n) = sum(binomial(n,2*k+1)*(4^k)*((2*k-1)!!)^2*(2*(n-2*k)-1)!!,k=0..floor(n/2)), with (2*n-1)!!:= A001147(n) (double factorials).
a(n) ~ Pi * 2^(n-1/2) * n^(n+1) / exp(n) * (1 - sqrt(2/(Pi*n))). - Vaclav Kotesovec, Mar 18 2014
2*(n+1)*(3+2*n)^2*a(n)-(4*n^2+8*n+1)*a(n+1)-(2*(n+4))*a(n+2)+a(n+3)=0. - Robert Israel, Feb 07 2018
MAPLE
f:= gfun:-rectoproc({2*(n+1)*(3+2*n)^2*a(n)-(4*n^2+8*n+1)*a(n+1)-(2*(n+4))*a(n+2)+a(n+3)=0, a(0)=0, a(1)=1, a(2)=6}, a(n), remember):
map(f, [$0 .. 30]); # Robert Israel, Feb 07 2018
MATHEMATICA
With[{nn=20}, CoefficientList[Series[ArcSin[2x]/(2(1-2x)^(3/2)), {x, 0, nn}], x] Range[0, nn]!] (* Harvey P. Dale, Mar 18 2019 *)
PROG
(PARI) x = 'x + O('x^40); concat(0, Vec(serlaplace(asin(2*x)/(2*(1-2*x)^(3/2))))) \\ Michel Marcus, Jun 18 2017
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Wolfdieter Lang, Sep 15 2008
STATUS
approved