%I #25 Aug 08 2015 05:05:11
%S 5,6,9,5,9,6,1,5,8,1,8,3,6,1,4,5,0,6,2,3,6,4,5,5,5,3,6,7,2,7,1,7,4,6,
%T 9,0,1,0,7,8,7,6,1,2,6,8,2,1,2,2,8,7,8,3,6,8,2,8,1,8,4,0,8,1,2,4,8,5,
%U 2,3,0,0,2,5,0,2,9,9,1,8,1,1,6,1,4,0,5,6,5,7,4,2,2,2,7,2,4,5,8,6,8
%N Decimal expansion of the integral of log(2+x^2+y^2)/((1+x^2)*(1+y^2)) dx dy over the square [0,1]x[0,1].
%C The computation of this integral is given by Bailey & Borwein as an example of the use of CAS packages (and additional tools) to simplify large symbolic expressions.
%H Vincenzo Librandi, <a href="/A244843/b244843.txt">Table of n, a(n) for n = 0..1000</a>
%H D. H. Bailey and J. M. Borwein, <a href="http://moodle.thecarma.net/jon/ontology.pdf">Experimental computation as an ontological game changer</a>, 2014, see p. 5.
%H D. H. Bailey, J. M. Borwein and A. D. Kaiser, <a href="http://www.carma.newcastle.edu.au/jon/auto.pdf">Automated Simplification of Large Symbolic Expressions</a>
%H Eric Weisstein's MathWorld, <a href="http://mathworld.wolfram.com/ClausensIntegral.html">Clausen's Integral</a>.
%H Eric Weisstein's MathWorld, <a href="http://mathworld.wolfram.com/Polylogarithm.html">Polylogarithm</a>.
%F Pi^2/8*log(2) - 7/48*zeta(3) + 11/24*Pi*Cl2(Pi/6) - 29/24*Pi*Cl2(5*Pi/6), where Cl2 is the Clausen function Cl2(t) = Sum_{n>0} sin(n*t)/n^2.
%e 0.56959615818361450623645553672717469010787612682122878368281840812485230025...
%t Clausen2[x_] := Im[PolyLog[2, Exp[x*I]]]; Pi^2/8*Log[2] - 7/48*Zeta[3] + 11/24*Pi*Clausen2[Pi/6] - 29/24*Pi*Clausen2[5*Pi/6] // RealDigits[#, 10, 101]& // First
%o (PARI) Cl2(x)=imag(polylog(2,exp(x*I)));
%o Pi^2/8*log(2) - 7/48*zeta(3) + 11/24*Pi*Cl2(Pi/6) - 29/24*Pi*Cl2(5*Pi/6) \\ _Charles R Greathouse IV_, Aug 27 2014
%Y Cf. A261027 (Cl_2(Pi/6)), A261028 (Cl_2(5*Pi/6)).
%K cons,nonn
%O 0,1
%A _Jean-François Alcover_, Jul 07 2014
|