%I #17 Sep 08 2022 08:45:11
%S 29,2077,389467,23511309071,250074841297,217632439585619,
%T 2271157731457180823,39331108008268763851,152552947614179997630583,
%U 30344459362884140864563052777
%N Numerators of successive approximations to zeta(3) = Sum_{k>0} 1/k^3, using Zeilberger's formula with s=2.
%H G. C. Greubel, <a href="/A084223/b084223.txt">Table of n, a(n) for n = 1..329</a>
%H D. Zeilberger, <a href="https://arxiv.org/abs/math/9804126">Faster and Faster convergent series for zeta(3)</a>, arXiv:math/9804126 [math.CO], 1998.
%F a(n) = numerator( Sum_{k=1..n} (1/4)*(-1)^(k-1)*(56*k^2 - 32*k + 5)/((2*k - 1)^2* binomial(3*k, k)*binomial(2*k, k)*k^3). - _G. C. Greubel_, Oct 08 2018
%p a:=n->add((1/4)*(-1)^(k-1)*(56*k^2-32*k+5)/((2*k-1)^2*binomial(3*k,k)*binomial(2*k,k)*k^3),k=1..n): seq(numer(a(n)),n=1..10); # _Muniru A Asiru_, Oct 09 2018
%t Table[Numerator[Sum[(1/4)*(-1)^(k - 1)*(56*k^2 - 32*k + 5)/((2*k - 1)^2* Binomial[3*k, k]*Binomial[2*k, k]*k^3), {k, 1, n}]], {n, 1, 30}] (* _G. C. Greubel_, Oct 08 2018 *)
%o (PARI) for(n=1,15,print1(numerator(sum(k=1,n,1/4*(-1)^(k-1)*(56*k^2 -32*k+5)/(2*k-1)^2/binomial(3*k,k)/binomial(2*k,k)/k^3)), ","))
%o (Magma) [Numerator((&+[(1/4)*(-1)^(k-1)*(56*k^2-32*k+5)/((2*k-1)^2*Binomial(3*k,k)*Binomial(2*k,k)*k^3): k in [1..n]])): n in [1..30]]; // _G. C. Greubel_, Oct 08 2018
%o (GAP) List(List([1..10],n->Sum([1..n],k->(1/4)*(-1)^(k-1)*(56*k^2-32*k+5)/((2*k-1)^2*Binomial(3*k,k)*Binomial(2*k,k)*k^3))),NumeratorRat); # _Muniru A Asiru_, Oct 09 2018
%Y Denominators are in A084224, decimal expansion is in A002117.
%Y Cf. A084225 (s=3).
%K nonn,frac
%O 1,1
%A _Ralf Stephan_, May 19 2003
|