70

I'm trying to find $\int x^x \, dx$, but the only thing I know how to do is this:

Let $u=x^x$.

$$\begin{align} \int x^x \, dx&=\int u \, du\\[6pt] &=\frac{u^2}{2}\\[6pt] &=\dfrac{\left(x^x\right)^2}{2}\\[6pt] &=\frac{x^{2x}}{2} \end{align}$$

But it's certain that this isn't the correct way to evaluate that, and the answer must be wrong.

Garmen1778
  • 2,338

9 Answers9

76

As noted in the comments, your derivation contains a mistake.

To answer the question, this function can not be integrated in terms of elementary functions. So there is no "simple" answer to your question, unless you are willing to consider a series approximation, obtained by expanding the exponential as a series:

$$\int{x^xdx} = \int{e^{\ln x^x}dx} = \int{\sum_{k=0}^{\infty}\frac{x^k\ln^k x}{k!}}dx$$

38

If you are willing to put bounds on your integral, it is possible to compute that $$\int_0^1 x^x\,dx = \sum_{n=1}^\infty \frac{(-1)^{n-1}}{n^n}.$$ Indeed, if you start like nbubis suggests, and make the substitution $u = -\log x$, you get that $$\int_0^1 x^x\,dx = \sum_{k=0}^\infty \frac{1}{k!}\int_0^1x^k(\log x)^k\,dx = \sum_{k=0}^\infty \frac{(-1)^k}{k!}\int_0^\infty e^{u(k+1)}u^k\,du$$$$ = \sum_{k=0}^\infty \frac{(-1)^k}{k!}\frac{1}{(k+1)^k}\int_0^\infty e^{u(k+1)}[(k+1)u]^k\,du.$$ If you then make the substitution $t = (k+1)u$ this becomes $$\sum_{k=0}^\infty \frac{(-1)^k}{k!}\frac{1}{(k+1)^k}\int_0^\infty e^tt^k\,dt = \sum_{k=0}^\infty \frac{(-1)^k}{(k+1)!}\frac{1}{(k+1)^k}\Gamma(k+1),$$ where $\Gamma$ is the usual Gamma function. Since $\Gamma(k+1) = k!$, the final expression is $$ \sum_{k=0}^\infty \frac{(-1)^k}{(k+1)^{k+1}} = \sum_{n=1}^\infty \frac{(-1)^{n-1}}{n^n}.$$ Similarly you can derive $\int_0^1 x^{-x}\,dx = \sum_{n=1}^\infty n^{-n}$. In don't think any further simplification is possible.

froggie
  • 9,941
  • 1
    This is nice. I like these kinds of answers. This gives, essentially, $$\int_0^1 x^{-zx},dx = \sum_{n=1}^{\infty} \frac{z^{n-1}}{n^n}.$$ – Antonio Vargas May 09 '12 at 06:51
  • 7
    And slightly more generally, $\int_0^1 x^{r-zx}\ dx = \sum_{n=1}^\infty \dfrac{z^{n-1}}{(r+n)^n}$ for $r > -1$. – Robert Israel May 09 '12 at 07:59
  • 7
    These identities for $\int_0^1 x^{-x}\ dx$ and $\int_0^1 x^x\ dx$ are sometimes called the "sophomore's dream". Look that up on Wikipedia. – Robert Israel May 09 '12 at 08:04
  • How do you know it's legal to switch the summation and the integral? I know you can do it with finite sums but I thought there were certain conditions under which it invalid to switch them. – sam-pyt Nov 20 '19 at 13:00
  • @sam-pyt late reply but you determine if they are switchable by testing if they satisfy Fubini's theorem and/or the dominated convergence theorem – Max0815 Nov 29 '21 at 14:28
22

let ${x}^{x} = {\left({e}^{\ln {x}} \right)}^{x} = {e}^{x \ln {x}}. $

By the series expansion of ${e}^{x}$: $${e}^{x \ln {x}} = \sum _{ n=0 }^{ \infty }{ \frac { { \left( x \ln{x} \right) }^{ n } }{ n! } }$$

Thus $$\int _{ 0 }^{ 1 }{ { x }^{ x } } dx=\sum _{ n=0 }^{ \infty }{ \int _{ 0 }^{ 1 }{ \frac { { { x }^{ n }\left( \ln {x} \right) }^{ n } }{ n! } } }=\sum _{ n=0 }^{ \infty }{ \frac { 1 }{ n! } \int _{ 0 }^{ 1 }{ { x }^{ n } } { \left( \ln { x } \right) }^{ n }dx }$$

Let $u = {\left(\ln {x} \right)}^{n} $, $dv = {x}^{n} dx $, $du = \frac{{n \left(\ln {x} \right)}^{n-1}}{x} dx$ and $v=\frac{{x}^{n+1}}{n+1}$, then using integration by parts, we arrive at

$$\lim _{ a\rightarrow 0 }{ \int _{ a }^{ 1 }{ { x }^{ n } } { \left( \ln { x } \right) }^{ n }dx } =\lim _{ a\rightarrow 0 }{ { \left[ \frac { { x }^{ n+1 } }{ n+1 } { \left( \ln { x } \right) }^{ n } \right] }_{ a }^{ 1 } } -\lim _{ a\rightarrow 0 }{ \int _{ a }^{ 1 }{ { \frac { n }{ n+1 } x }^{ n } } { \left( \ln { x } \right) }^{ n-1 } } dx$$

which becomes $$\lim _{ a\rightarrow 0 }{ \int _{ a }^{ 1 }{ { x }^{ n } } { \left( \ln { x } \right) }^{ n }dx } =-\int _{ 0 }^{ 1 }{ { \frac { n }{ n+1 } x }^{ n } } { \left( \ln { x } \right) }^{ n-1 }dx = \frac{{(-1)}^{n}n!}{{(n+1)}^{n+1}}$$

Therefore $$\int _{ 0 }^{ 1 }{ { x }^{ x } } dx=\sum _{ n=1 }^{ \infty }{ \frac { { \left( -1 \right) }^{ n-1 } }{ { n }^{ n } } }$$

Lord_Farin
  • 17,743
  • 3
    As to the content, you didn't answer the question, for it was asking about the indefinite, not the definite, integral. – Lord_Farin Jan 03 '15 at 10:09
8

On can find a compendium of properties of the special function : $$\text{Sphd}(\alpha\:;\:x)=\int_0^x t^{\alpha\:t}dt$$ and the particular case : $$\int x^x dx = \text{Sphd}(1\:;\:x) +\text{constant}$$ in https://fr.scribd.com/doc/34977341/Sophomore-s-Dream-Function

JJacquelin
  • 66,221
  • 3
  • 37
  • 87
4

The integral $\int{x^x}{dx}$ can be expressed as a double series. I asked about this series form here and the answers there show it is correct and my own answer there shows you can differentiate this back to get a power series for $x^x$: $$ \int{x^x}{dx}=\sum _{n=1}^{\infty } \sum _{k=0}^{n-1} \frac{x^n \log ^k(x) (-1)^{1+n+k}}{n^{n-k}\ k!} $$

flinty
  • 13
  • Note this is valid for $x>0 \land x\neq1$. The antiderivative $\int{x^{-x}dx}$ is very similar: $$\int{x^{-x}}{dx}=\sum _{n=1}^{\infty } \sum _{k=0}^{n-1} \frac{x^n \log ^k(x) (-1)^{k}}{n^{n-k}\ k!}$$ – flinty Oct 10 '20 at 15:21
3

Many people have pointed out that the integral you are looking for is equivalent to,

$$\sum_{n}^{\infty} \frac{1}{n!} \int_{0}^{x}x^{n}\ln(x)^ndx$$

But the integral within this equation can be simplified to

$$ \int_0^x x^n \ln(x)^ndx = (-1)^n (n+1)^{-1-n} \Gamma \left(n+1,(n+1)\ln \left(\frac{1}{x}\right)\right) $$

Where

$$ \Gamma \left(n+1,\left(n+1\right)\ln \left(\frac{1}{x}\right)\right) = n! e^{(n+1)\ln(x)} \sum_{k=0}^{n} \frac{(n+1)^k \ln \left(\frac{1}{x}\right)^k}{k!} $$

Is the incomplete Gamma function. Simplify the first equation, and you will get,

$$ \int_0^x x^x dx = \sum_{n=0}^\infty \frac{(-1)^n \Gamma \left(n+1,-(n+1) \ln({x})\right)} {n!(n+1)^{n+1}} $$

A demonstration of this function may be found on the desmos graphing calculator: https://www.desmos.com/calculator/2nfxrv0iba

gt6989b
  • 54,422
1

Start from the opposite task.

If $\displaystyle \int x^x \, dx=F(x)$ then $\displaystyle F'(x)=x^x$

First we need to find asymptotic evaluation of the integral. Let us take it in the form

$$F(x)=x^xg(x)$$

So it has to be:

$$F'(x)=x^x((1+\ln(x))g(x)+g'(x))=x^x$$

From there it is sufficient to take $g(x) \sim \frac{1}{1+\ln(x)}$

So we can start our journey:

$$F(x)=x^x(\frac{1}{1+\ln(x)}+f(g(x)))$$

If you calculate the derivative of this you have

$$g'(x)f'(g(x))-\frac{1}{x(\ln(x)+1)^2}+(\ln(x)+1)f(g(x))=0$$

For the purpose of cancellation if it best to take

$$g'(x)=\ln(x)+1$$

meaning

$$g(x)=x\ln(x)$$

Now we continue using the steps that are revealing the integral structure.

$$F(x)=x^x(\frac{1}{\ln(x)+1}+f(x\ln(x)))$$

Take derivative once more and you have got $$f(x\ln(x))=\frac{1}{x(1+\ln(x))^3}-f'(x\ln(x))$$ or $$F(x)=x^x(\frac{1}{\ln(x)+1}+\frac{1}{x(\ln(x)+1)^3}-f'(x\ln(x)))$$

We can then write:

$$F(x)=x^x(\frac{1}{\ln(x)+1}+\sum_{n=1}^{\infty}f_n(x))$$

where

$$\displaystyle f_{n}=-\frac{f_{n-1}'}{1+\ln(x)},\,f_0=\frac{1}{1+\ln(x)}$$

From $x=0$ to $1$ you can use probably more suitably $F(x)=xg(\ln(x))$

The derivation is similar to the one given above.

1

I will prove that $x^x$ has no elementary antiderivative, as many have claimed. I will use Risch's decision semi-algorithm for exponential polynomials in trascendental exponential extensions. Note that $g=x^x=\exp(x\log(x))$. We can then express our integrand as an element of the differential field $$\mathbb{Q}(x,\log(x),g)=\mathbb{Q}(x,\log(x))(g),$$ where we take $g$ as the last extension. This is, as we call it, a field of trascendental elementary functions. Then, $g$ can also be seen as a polynomial in $\mathbb{Q}(x,\log(x))[g]$, and if $g$ has an elementary antiderivative, then it has the form $\int g=q_1g$, where $q_1\in \mathbb{Q}(x,\log(x))$, and $q_1$ satisfies the Risch differential equation $$q_1'+(1+\log(x))q_1=1.$$ In this equation, $1+\log(x)$ is the derivative of the function inside the exponential $g=\exp(x\log(x))$, and the 1 on the right is the coefficient of our "monomial" $g$. To solve this equation, due to $\log(x)$ being trascendental over $\mathbb{Q}(x)$, we can equate coefficients over the polynomial ring $\mathbb{Q}(x)[\log(x)]$. Doing this we obtain $q_1'+q_1=1$ for degree 0 and $q_1=0$ for degree 1 in $\log(x)$. This is impossible, as $q_1=0$ implies $q_1'+q_1=0$. This means that $\int x^x dx$ has no elementary antiderivative.

The specific reason to take $\mathbb{Q}(x,\log(x))(g)$ as our differential field is that it is a field of trascendental elementary functions. A field of trascendental elementary functions takes as its base a field of rational functions, which could be $\mathbb{Q}(x)$, $\mathbb{R}(x)$ or $\mathbb{C}(x)$, for example. Then, we can attach a new function to this field, and if it is trascendental logarithmic or exponential over $\mathbb{Q}(x)$, it is said to be a trascendental elementary extension. For Risch's decision algorithm to work in a trascendental extension, we need that every extension of the base field be trascendental (exponential or logarithmic). As a note, the other type of elementary extension that could be presented is the algebraic extension (which is the opposite of a trascendental extension, and includes functions such as radicals).

In general, an elementary function is any function that can be seen as an element of a field of elementary functions. That is, an extension field with finite new functions over $\mathbb{Q}(x)$ such that each extension is exponential, logarithmic or algebraic over it.

  • A function $g$ is logarithmic over $F$ if there exists $f\in F$ such that $g'=f'/f$. Notation: $g=\log(f)$.
  • A function $g$ is exponential over $F$ if there exists $f\in F$ such that $g'=gf'$. Notation: $g=\exp(f)$.
  • A function $g$ is algebraic over $F$ if there exists a polynomial $p(z)\in F[z]$ such that $p(g)=0$. For example, $g=\sqrt{x}$ is algebraic over $\mathbb{Q}(x)$ because the polynomial $p(z)=z^2-x\in\mathbb{Q}(x)[z]$ is such that $p(g)=0$.

What we really proved is that $\int x^x dx$ can never be represented as a member of any field of elementary functions. Thus, it cannot be integrated in finite elementary symbols, which is why any answer to the integral requires infinitely many symbols, or a new notation. However, I obviously used some heavy results. If you want to see the proof of these results, I recommend researching differential algebra. There is a book by the name of "Algorithms for Computer Algebra" by Keith Geddes, which is very complete in regards to integration of trascendental elementary functions, and presents the proof in a very comprehensible level. Most of it is in chapters 11 and 12, but other important notions can be found in previous chapters.

  • Why should $\int g=q_1g$? This seems like a statement requiring more work to show. Note that $\exp(\cdots)$ is not a polynomial, so, there are many aspects of your answer that are unclear to me – FShrike Oct 20 '23 at 00:01
  • 1
    @FShrike $\exp(x\log(x))$ is not a polynomial in $\mathbb{Q}[x]$, which means that it is not a polynomial over the variable $x$. It is, however, a polynomial over the variable $g$. The reason why $\int g=q_1g$ is a result of the way "polynomials" in an exponential variable behave, which is in turn a consequence of a big result called "Liouville's principle". As I mentioned at the end of my answer, there are a lot of details that I left out because a whole book is needed to explain them. That is the reason why it's hard to come across of a proof that an integral is non-elementary. – Math Enjoyer Oct 20 '23 at 00:17
0

Here's an addition. As found on my site, I'm mostly surprised to not find an answer that was derived by simple methods (integration by parts using my calculator and for the series my simple derivation of the one found on Wolfram Alpha):

$$\int x^x\ dx = x^{x+1}-\int \ln{(x)x^{x+1}}\ dx-\int x^{x+1}\ dx+C =$$ $$C-\sum_{n=0}^{\infty}\frac{1}{n!(-n-1)^{n+1}}\int_{(-n-1)\ln{(x)}}^{\infty}e^{-t}t^n\ dt$$