OFFSET
0,1
COMMENTS
Let N be a positive integer divisible by 4. We have the asymptotic expansion 2*((Pi/2 - log(1 + sqrt(2)))/(2*sqrt(2)) - Sum_{k = 0..N/4 - 1} (-1)^k/(4*k + 3)) ~ 1/N - 1/N^2 - 3/N^3 + 11/N^4 + 57/N^5 - ..., where the sequence of coefficients [1, -1, -3, 11, 57, ...] is A212435. This follows from Borwein et al., Lemma 2 with f(x) = 1/x and then set x = N/4 and h = 3/4. An example is given below. Cf. A181048. - Peter Bala, Sep 23 2016
LINKS
Gheorghe Coserea, Table of n, a(n) for n = 0..2015
J. M. Borwein, P. B. Borwein, and K. Dilcher, Pi, Euler numbers and asymptotic expansions, Amer. Math. Monthly, 96 (1989), 681-687.
Eric W. Weisstein, Euler's Series Transformation.
FORMULA
Equals Integral_{x=0..1} (x^2 dx)/(1+x^4).
Equals (1/2) * Integral_{x = 0..Pi/4} sqrt(tan(x)) dx. Cf. A247719. - Peter Bala, Sep 23 2016
Equals Sum_{n >= 0} 2^(n-1)*n!/(Product_{k = 0..n} 4*k + 3) = Sum_{n >= 0} 2^(n-1)*n!/A008545(n+1) (apply Euler's series transformation to Sum_{k >= 0} (-1)^k/(4*k+3)). - Peter Bala, Dec 01 2021
From Peter Bala, Mar 03 2024: (Start)
Continued fraction: 1/(3 + 3^2/(4 + 7^2/(4 + 11^2/(4 + 15^2/(4 + ... ))))) due to Euler.
Equals (1/3)*hypergeom([3/4, 1], [7/4], -1).
Gauss's continued fraction: 1/(3 + 3^2/(7 + 4^2/(11 + 7^2/(15 + 8^2/(19 + 11^2/(23 + 12^2/(27 + 15^2/(31 + 16^2/(35 + 19^2/(39 + ... )))))))))). (End)
Equals Integral_{x=1..oo, y=1..oo} 1/(x^4 + y^4) dx. - Vaclav Kotesovec, Jun 13 2024
EXAMPLE
0.2437477471996805241799750836323027110...
From Peter Bala, Sep 23 2016: (Start)
At N = 100000 the truncated series 2*Sum_{k = 0..N/4 - 1} (-1)^k/(4*k + 3) = 0.4874(8)5494(4)9936(4)048(24)99(444)67(625)6... to 32 digits. The bracketed numbers show where this decimal expansion differs from that of 2*A181049. The numbers 1, -1, -3, 11, 57, -361 must be added to the bracketed numbers to give the correct decimal expansion to 32 digits: 2*( Pi/2 - log(1+sqrt(2)))/(2*sqrt(2) ) = 0.4874(9)5494(3)9936(1)048(35)99(501)67(264)6.... (End)
MATHEMATICA
First@ RealDigits[N[(Pi/2 - Log[1 + Sqrt@ 2])/(2 Sqrt@ 2), 105]] (* Michael De Vlieger, Oct 07 2015 *)
PROG
(PARI)
default(realprecision, 106);
eval(vecextract(Vec(Str(sumalt(n=0, (-1)^(n)/(4*n+3)))), "3..-2")) \\ Gheorghe Coserea, Oct 06 2015
(PARI) (Pi/2 - log(1+sqrt(2)))/(2*sqrt(2)) \\ G. C. Greubel, Nov 28 2017
(Magma) C<i> := ComplexField(); [(Pi(C)/2 - Log(1+Sqrt(2)))/(2*Sqrt(2))]; // G. C. Greubel, Nov 28 2017
CROSSREFS
KEYWORD
cons,nonn
AUTHOR
Jonathan D. B. Hodgson, Oct 01 2010, Oct 05 2010
STATUS
approved