%I #15 Sep 28 2019 06:07:41
%S 0,16,576,25600,1254400,65028096,3497066496,192980975616,
%T 10855179878400,619683355033600,35792910586740736,2087229562810269696,
%U 122682715414070296576,7259332273021911040000,432004345063916175360000,25835779854133582469529600
%N a(n) = A069722(n+1)^2.
%F G.f.: 1/AGM(1, (1-64*x)^(1/2)).
%F E.g.f.: 1 + Sum[n>=0, a(n)*x^(2n)/(2n)! ] = BesselI(0, 4x)^2. - _Ralf Stephan_, Jan 11 2005
%F From _Vaclav Kotesovec_, Sep 28 2019: (Start)
%F For n > 0, a(n) = 2^(2*n) * binomial(2*n, n)^2.
%F a(n) ~ 2^(6*n) / (Pi*n). (End)
%t Flatten[{0, Table[2^(2*n) * Binomial[2*n, n]^2, {n, 1, 20}]}] (* _Vaclav Kotesovec_, Sep 28 2019 *)
%t CoefficientList[Series[-1 + 2*EllipticK[1 - 1/(1 - 64*x)] / (Pi*Sqrt[1 - 64*x]), {x, 0, 20}], x] (* _Vaclav Kotesovec_, Sep 28 2019 *)
%Y Cf. A069722.
%K nonn
%O 0,2
%A _Benoit Cloitre_, Jan 03 2004
|