A363516 Denominator of log(2) + (-1/4)^n*Integral_{x=0..1} (1 - x)^(4*n+2)/(1 + x^2) dx.
1, 120, 20160, 5765760, 313657344, 14898723840, 27413651865600, 4769975424614400, 4731815621217484800, 175077177985046937600, 28712657189547697766400, 2469288518301102007910400, 12998334760337000969640345600, 86113967787232631423867289600, 5511293938382888411127506534400
Offset: 0
Examples
n A363515(n)/a(n) approximate value - ------------------- ----------------- 0 1 1 1 79/120 0.6583333333... 2 14087/20160 0.6987599206... 3 3990557/5765760 0.6921129218... 4 217474889/313657344 0.6933518158... ... From _M. F. Hasler_, Jul 07 2023: (Start) Let f[n] = (-1/4)^n*(1 - x)^(4*n+2)/(1 + x^2), the rational fraction to be integrated from 0 to 1. We have: f[0] = 1 - 2*x/(1 + x^2), with primitive F[0] = x/2 - log(1 + x^2), whence an integral equal to 1/2 - log(2), and a(0) = 2 (denominator). f[1] = -x^4/4 + (3/2)*x^3 - (7/2)*x^2 + (7/2)*x - 1/4 - 2*x/(1 + x^2), and the term-wise integration of the polynomial part gives -1/20 + 3/8 - 7/6 + 7/4 - 1/4 = 79/120, whence A363515(1) = 79 and a(1) = 120. f[2] = (1/16)*x^8 - (5/8)*x^7 + (11/4)*x^6 - (55/8)*x^5 + (83/8)*x^4 - (71/8)*x^3 + (11/4)*x^2 + (11/8)*x + 1/16 - 2*x/(1 + x^2), so the integration gives 1/144 - 5/64 + 11/28 - 55/48 + 83/40 - 71/32 + 11/12 + 11/16 + 1/16 - log(2) = 14087/20160 - log(2), whence A363515(2) = 14087 and a(2) = 20160, etc. (End)
Links
- Mathematica Stack Exchange, How to prove using Mathematica that the sequence converges to log(2)?
Programs
-
Mathematica
Denominator[Simplify[Table[Log[2]+(-1)^n Integrate[(1-x)^(4n+2)/(1+x^2),{x,0,1}]/4^n,{n,0,14}]]]
Formula
Denominator of log(2) + HypergeometricPFQ([1/2, 1, 1], [2*(1 + n), 5/2 + 2*n], -1)/((3 + 4*n)*(-4)^n).
Limit_{n->oo} A363515(n)/a(n) = log(2).
Comments