login
a(n) = (2n-1)!! * Sum_{k=0..n-1}(-1)^k/(2k+1).
20

%I #79 Sep 08 2022 08:44:48

%S 0,1,2,13,76,789,7734,110937,1528920,28018665,497895210,11110528485,

%T 241792844580,6361055257725,163842638377950,4964894559637425,

%U 147721447995130800,5066706567801827025,171002070002301095250,6548719685561840296125,247199273204273879989500

%N a(n) = (2n-1)!! * Sum_{k=0..n-1}(-1)^k/(2k+1).

%C (2*n + 1)!!/a(n+1), n>=0, is the n-th approximant for William Brouncker's continued fraction for 4/Pi = 1 + 1^2/(2 + 3^2/(2 + 5^2/(2 + ... ))) See the C. Brezinski and J.-P. Delahaye references given under A142969 and A142970, respectively. The double factorials (2*n + 1)!! = A001147(n+1) enter. - _Wolfdieter Lang_, Oct 06 2008

%D A. E. Jolliffe, Continued Fractions, in Encyclopaedia Britannica, 11th ed., pp. 30-33; see p. 31.

%H Vincenzo Librandi, <a href="/A024199/b024199.txt">Table of n, a(n) for n = 0..200</a>

%H Peter Luschny, <a href="http://www.luschny.de/math/factorial/hadamard/HadamardsGammaFunctionMJ.html">Is the Gamma function misdefined? Hadamard versus Euler - Who found the better Gamma function?</a>

%H Wikipedia, <a href="http://en.wikipedia.org/wiki/William_Brouncker,_2nd_Viscount_Brouncker">William Brouncker, 2nd Viscount Brouncker</a>

%F a(n) = s(1)s(2)...s(n)(1/s(1) - 1/s(2) + ... + c/s(n)) where c=(-1)^(n+1) and s(k) = 2k-1 for k = 1, 2, 3, ... (was previous definition). - _Clark Kimberling_

%F D-finite with recurrence a(0) = 0, a(1) = 1, a(n+1) = 2*a(n) + (2*n-1)^2*a(n-1). - _N. J. A. Sloane_, Jul 19 2002

%F a(n) + A024200(n) = A001147(n) = (2n-1)!!. - _Max Alekseyev_, Sep 23 2007

%F a(n)/A024200(n) -> Pi/(4-Pi) as n -> oo. - _Max Alekseyev_, Sep 23 2007

%F From _Wolfdieter Lang_, Oct 06 2008: (Start)

%F E.g.f. for a(n+1), n>=0: (sqrt(1-2*x)+arcsin(2*x)*sqrt(1+2*x)/2)/((1-4*x^2)^(1/2)*(1-2*x)). From the recurrence, solving (1-4*x^2)y''(x)-2*(8*x+1)*y'(x)-9*y=0 with inputs y(0)=1, y'(0)=2.

%F a(n+1) = A003148(n) + A143165(n), n>=0 (from the two terms of the e.g.f.). (End)

%F From _Johannes W. Meijer_, Nov 12 2009: (Start)

%F a(n) = (-1)^(n-1)*(2*n-3)!! + (2*n-1)*a(n-1) with a(0) = 0.

%F a(n) = (2*n-1)!!*sum((-1)^(k)/(2*k+1), k=0..n-1)

%F (End)

%F E.g.f.: Pi/4/sqrt(1-2*x) - 1/2*log(2*x+sqrt(4*x^2-1))/sqrt(2*x-1). - _Vaclav Kotesovec_, Mar 18 2014

%F a(n) ~ Pi * 2^(n-3/2) * n^n / exp(n). - _Vaclav Kotesovec_, Mar 18 2014

%F a(n) = (2*H(n+1/2)-Gamma(n+1/2))*2^(n-2)*sqrt(Pi) with H(x) the Hadamard factorial (see the link section). - _Cyril Damamme_, Jul 19 2015

%F a(n) = A135457(n) - (-1)^n A001147(n-1). - _Cyril Damamme_, Jul 19 2015

%F a(n) = (Pi + (-1)^n*(Psi(n/2 + 1/4) - Psi(n/2 + 3/4)))*Gamma(n+1/2)*2^(n-2)/sqrt(Pi). - _Robert Israel_, Jul 20 2015

%F a(n) = A167576(n) - A135457(n). - _Cyril Damamme_, Jul 22 2015

%F a(n)/A001147(n) -> Pi/4 as n -> oo. - _Daniel Suteu_, Jul 21 2016

%F From _Peter Bala_, Nov 15 2016: (Start)

%F Conjecture: a(n) = 1/2*Sum_{k = 0..2*n-1} (-1)^(n-k+1)*k!*(2*n - 2*k - 3)!!, where the double factorial of an odd integer (positive or negative) may be defined in terms of the gamma function as (2*N - 1)!! = 2^((N+1)/2)*Gamma(N/2 + 1)/sqrt(Pi).

%F E.g.f. 1/2*arcsin(2*x)/sqrt(1 - 2*x) = x + 2*x^2/2! + 13*x^3/3! + 76*x^4/4! + .... (End)

%e a(3) = (2*3 - 1)!! * Sum_{k=0..2} (-1)^k/(2k + 1) = 5!! * (1/(2*0 + 1) - 1/(2*1 + 1) + 1/(2*2 + 1)) = 5*3*1*(1/1 - 1/3 + 1/5) = 15 - 5 + 3 = 13. Notice that the first factor always cancels the common denominator of the sum. - _Michael B. Porter_, Jul 22 2016

%p a := proc(n) option remember; if n=0 then 0 elif n=1 then 1 else 2*a(n-1)+(2*n-3)^2* a(n-2) fi end: seq(a(n), n=0..20); # _Peter Luschny_, Nov 16 2016 after _N. J. A. Sloane_

%t f[k_] := (2 k - 1) (-1)^(k + 1)

%t t[n_] := Table[f[k], {k, 1, n}]

%t a[n_] := SymmetricPolynomial[n - 1, t[n]]

%t Table[a[n], {n, 1, 22}] (* A024199 signed *)

%t (* _Clark Kimberling_, Dec 30 2011 *)

%t RecurrenceTable[{a[n+1] == 2*a[n] + (2*n-1)^2*a[n-1],a[0] == 0, a[1] == 1},a,{n,0,20}] (* _Vaclav Kotesovec_, Mar 18 2014 *)

%t CoefficientList[Series[Pi/4/Sqrt[1-2*x] - 1/2*Log[2*x+Sqrt[4*x^2-1]]/Sqrt[2*x-1], {x, 0, 20}], x] * Range[0, 20]! (* _Vaclav Kotesovec_, Mar 18 2014 *)

%o (Magma) [0] cat [ n le 2 select (n) else 2*Self(n-1)+(2*n-3)^2*Self(n-2): n in [1..25] ]; // _Vincenzo Librandi_, Feb 17 2015

%Y Cf. A004041, A000407.

%Y From _Johannes W. Meijer_, Nov 12 2009: (Start)

%Y Cf. A007509 and A025547.

%Y Equals first column of A167584.

%Y Equals row sums of A167591.

%Y Equals first right hand column of A167594.

%Y (End)

%Y Cf. A167576 and A135457.

%K nonn,easy

%O 0,3

%A _Clark Kimberling_

%E Edited by _N. J. A. Sloane_, Jul 19 2002

%E New name from _Cyril Damamme_, Jul 19 2015