OFFSET
0,1
COMMENTS
According to Steven Finch, no exact expression of this moment is known.
REFERENCES
Steven R. Finch, Mathematical Constants, Cambridge University Press, 2003, Section 5.16 Extreme value constants, p. 365.
FORMULA
integral_(-infinity..infinity) 8*x^2*F(x)^7*f(x) dx - mu(8)^2, where f(x) is the normal (0,1) density and F(x) its cumulative distribution, mu(8) being the moment A243961.
EXAMPLE
0.3728971432867289942202112287621146...
MATHEMATICA
digits = 101; m0 = 5; dm = 5; f[x_] := 1/ Sqrt[2*Pi]*Exp[-x^2/2]; F[x_] := 1/2*Erf[x/Sqrt[2]] + 1/2; Clear[mu8]; mu8[m_] := mu8[m] = 8*NIntegrate[x*F[x]^7*f[x], {x, -m , m}, WorkingPrecision -> digits+5, MaxRecursion -> 20]; mu8[m0]; mu8[m = m0+dm]; While[RealDigits[mu8[m]] != RealDigits[mu8[m-dm]], Print["m1 = ", m]; m = m+dm]; m8 = mu8[m]; Clear[v, m]; v[m_] := v[m] = 8*NIntegrate[x^2*F[x]^7*f[x], {x, -m , m}, WorkingPrecision -> digits+5, MaxRecursion -> 20]; v[m0]; v[m = m0+dm]; While[RealDigits[v[m]] != RealDigits[v[m-dm]], Print["m2 = ", m]; m = m+dm]; v8 = v[m]-m8^2; RealDigits[v8, 10, digits] // First
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Jean-François Alcover, Jun 16 2014
STATUS
approved