Sophomore's dream

In mathematics, sophomore's dream is a name occasionally used for the identities (especially the first)

$\begin{align} \int_0^1 x^{-x}\,\mathrm{d}x &= \sum_{n=1}^\infty n^{-n}&&(\scriptstyle{= 1.29128599706266354040728259059560054149861936827\dots)} \\ \int_0^1 x^x \,\mathrm{d}x &= \sum_{n=1}^\infty (-1)^{n+1}n^{-n} = - \sum_{n=1}^\infty (-n)^{-n} &&(\scriptstyle{= 0.78343051071213440705926438652697546940768199014\dots}) \end{align}$

discovered in 1697 by Johann Bernoulli.

The name "sophomore's dream", which appears in (Borwein, Bailey & Girgensohn 2004), is in contrast to the name "freshman's dream" which is given to the incorrect^{[note 1]} equation (x + y)ⁿ = xⁿ + yⁿ. The sophomore's dream has a similar too-good-to-be-true feel, but is true.

Proof

Graph of the functions y = x^x and y = x^−x on the interval x ∈ (0, 1].

We prove the second identity; the first is completely analogous.

The key ingredients of the proof are:

• Write x^x = exp(x log x).
• Expand exp(x log x) using the power series for exp.
• Integrate termwise.
• Integrate by substitution.

Expand x^x as

$x^x = \exp(x \log x) = \sum_{n=0}^\infty \frac{x^n(\log x)^n}{n!}.$

Therefore, we have : $\int_0^1 x^x\,\mathrm{d}x = \int_0^1 \sum_{n=0}^\infty \frac{x^n(\log x)^n}{n!} \,\mathrm{d}x.$

By uniform convergence of the power series, we may interchange summation and integration

$\int_0^1 x^x\,\mathrm{d}x = \sum_{n=0}^\infty \int_0^1 \frac{x^n(\log x)^n}{n!} \,\mathrm{d}x.$

To evaluate the above integrals we perform the change of variable in the integral $\scriptstyle x=\exp\, \left(-\frac{u}{n+1}\right)$, with $\scriptstyle 0 < u < \infty$, giving us

$\int_0^1 x^n(\log\, x)^n\,\mathrm{d}x = (-1)^n (n+1)^{-(n+1)} \int_0^\infty u^n e^{-u}\,\mathrm{d}u.$

By the well-known Euler's integral identity for the Gamma function

$\int_0^\infty u^n e^{-u}\,\mathrm{d}u=n!$

so that:

$\int_0^1 \frac{x^n (\log x)^n}{n!}\,\mathrm{d}x = (-1)^n (n+1)^{-(n+1)}.$

Summing these (and changing indexing so it starts at n = 1 instead of n = 0) yields the formula.

Historical proof

The original proof, given in Bernoulli (1697), and presented in modernized form in Dunham (2005), differs from the one above in how the termwise integral $\int_0^1 x^n(\log\, x)^n\,\mathrm{d}x$ is computed, but is otherwise the same, omitting technical details to justify steps (such as termwise integration). Rather than integrating by substitution, yielding the Gamma function (which was not yet known), Bernoulli used integration by parts to iteratively compute these terms.

The integration by parts proceeds as follows, varying the two exponents independently to obtain a recursion. An indefinite integral is computed initially, omitting the constant of integration $+ C$ both because this was done historically, and because it drops out when computing the definite integral. One may integrate $\scriptstyle \int x^m (\ln x)^n\,\mathrm{d}x$ by taking u = (ln x)ⁿ and dv = x^m dx, which yields:

$\begin{align} \int x^m (\ln x)^n\,\mathrm{d}x & = \frac{x^{m+1}(\ln x)^n}{m+1} - \frac{n}{m+1}\int x^{m+1} \frac{(\ln x)^{n-1}}{x}\,\mathrm{d}x \qquad\mbox{(for }m\neq -1\mbox{)}\\ & = \frac{x^{m+1}}{m+1}(\ln x)^n - \frac{n}{m+1}\int x^m (\ln x)^{n-1}\,\mathrm{d}x \qquad\mbox{(for }m\neq -1\mbox{)} \end{align}$

(also in the list of integrals of logarithmic functions). This reduces the power on the logarithm in the integrand by 1 (from $n$ to $n-1$) and thus one can compute the integral inductively, as

$\int x^m (\ln x)^n\,\mathrm{d}x = \frac{x^{m+1}}{m+1} \cdot \sum_{i=0}^n (-1)^i \frac{(n)_i}{(m+1)^i} (\ln x)^{n-i}$

where (n) _i denotes the falling factorial; there is a finite sum because the induction stops at 0, since n is an integer.

In this case m = n, and they are integers, so

$\int x^n (\ln x)^n\,\mathrm{d}x = \frac{x^{n+1}}{n+1} \cdot \sum_{i=0}^n (-1)^i \frac{(n)_i}{(n+1)^i} (\ln x)^{n-i}.$

Integrating from 0 to 1, all the terms vanish except the last term at 1,^{[note 2]} which yields:

$\int_0^1 \frac{x^n (\ln x)^n}{n!}\,\mathrm{d}x = \frac{1}{n!}\frac{1^{n+1}}{n+1} (-1)^n \frac{(n)_n}{(n+1)^n} = (-1)^n (n+1)^{-(n+1)}.$

From a modern point of view, this is (up to a scale factor) equivalent to computing Euler's integral identity for the Gamma function, $\Gamma(n+1) = n!\,$ on a different domain (corresponding to changing variables by substitution), as Euler's identity itself can also be computed via an analogous integration by parts.

See also

• Series (mathematics)

Notes

1. Incorrect unless one is working over a field or unital commutative ring of prime characteristic n or a factor of n. The correct result is given by the binomial theorem.
2. All the terms vanish at 0 because $\scriptstyle\lim_{x \to 0^+} x^m (\ln x)^n \, = \, 0$ by l'Hôpital's rule (Bernoulli omitted this technicality), and all but the last term vanish at 1 since ln(1) = 0.

References

Formula

• Johann Bernoulli, 1697, collected in Johannis Bernoulli, Opera omnia, vol. 3, pp. 376–381
• Borwein, Jonathan; Bailey, David H.; Girgensohn, Roland (2004), Experimentation in Mathematics: Computational Paths to Discovery, pp. 4, 44, ISBN 978-1-56881-136-9
• Dunham, William (2005), "3: The Bernoullis (Johann and $x^x$)", The Calculus Gallery, Masterpieces from Newton to Lebesgue, Princeton, NJ: Princeton University Press, pp. 46–51, ISBN 978-0-691-09565-3
• OEIS, (sequence A083648 in OEIS) and (sequence A073009 in OEIS)
• Pólya, George; Szegő, Gábor (1998), "part I, problem 160", Problems and Theorems in Analysis, p. 36, ISBN 978-3-54063640-3
• Weisstein, Eric W., "Sophomore's Dream", MathWorld.

• Max R. P. Grossmann (2013): Sophomore's dream. 200,000 digits of the first constant
• Danil Krotkov (2015): Generalized identity.

Function

• Literature for x^x and Sophomore's Dream, Tetration Forum, 03/02/2010
• The Coupled Exponential, Jay A. Fantini, Gilbert C. Kloepfer, 1998
• Sophomore's Dream Function, Jean Jacquelin, 2010, 13 pp.
• Lehmer, D. H. (1985). "Numbers associated with Stirling numbers and x^x". Rocky Mountain Journal of Mathematics 15: 461. doi:10.1216/RMJ-1985-15-2-461.
• Gould, H. W. (1996). "A Set of Polynomials Associated with the Higher Derivatives of y = x^x". Rocky Mountain Journal of Mathematics 26: 615. doi:10.1216/rmjm/1181072076.