$\newcommand{\+}{^{\dagger}}%
\newcommand{\angles}[1]{\left\langle #1 \right\rangle}%
\newcommand{\braces}[1]{\left\lbrace #1 \right\rbrace}%
\newcommand{\bracks}[1]{\left\lbrack #1 \right\rbrack}%
\newcommand{\dd}{{\rm d}}%
\newcommand{\isdiv}{\,\left.\right\vert\,}%
\newcommand{\ds}[1]{\displaystyle{#1}}%
\newcommand{\equalby}[1]{{#1 \atop {= \atop \vphantom{\huge A}}}}%
\newcommand{\expo}[1]{\,{\rm e}^{#1}\,}%
\newcommand{\floor}[1]{\,\left\lfloor #1 \right\rfloor\,}%
\newcommand{\ic}{{\rm i}}%
\newcommand{\imp}{\Longrightarrow}%
\newcommand{\ket}[1]{\left\vert #1\right\rangle}%
\newcommand{\pars}[1]{\left( #1 \right)}%
\newcommand{\partiald}[3][]{\frac{\partial^{#1} #2}{\partial #3^{#1}}}
\newcommand{\pp}{{\cal P}}%
\newcommand{\root}[2][]{\,\sqrt[#1]{\,#2\,}\,}%
\newcommand{\sech}{\,{\rm sech}}%
\newcommand{\sgn}{\,{\rm sgn}}%
\newcommand{\totald}[3][]{\frac{{\rm d}^{#1} #2}{{\rm d} #3^{#1}}}
\newcommand{\ul}[1]{\underline{#1}}%
\newcommand{\verts}[1]{\left\vert #1 \right\vert}%
\newcommand{\yy}{\Longleftrightarrow}$
$\ds{\pp\pars{\mu} \equiv \int_{0}^{1}{x^{\mu} - 1 \over \ln\pars{x}}\,\dd x}$
$$
\pp'\pars{\mu} \equiv \int_{0}^{1}{x^{\mu}\ln\pars{x} \over \ln\pars{x}}\,\dd x
=
\int_{0}^{1}x^{\mu}\,\dd x = {1 \over \mu + 1}
\quad\imp\quad
\pp\pars{\mu} - \overbrace{\pp\pars{0}}^{=\ 0} = \ln\pars{\mu + 1}
$$
$$
\pp\pars{7} = \color{#0000ff}{\large\int_{0}^{1}{x^{7} - 1 \over \ln\pars{x}}
\,\dd x}
=
\ln\pars{7 + 1} = \ln\pars{8} = \color{#0000ff}{\large 3\ln\pars{2}}
$$
To formally justify this integration, we can do the following. Put $f_1(x) = \sqrt x$ ,and for each $n\ge 1$, put $f_{n+1}(x) = \sqrt{x + f_{n}(x)}$. We want to show that $f(x) = \lim_{n\to\infty}f_n(x)$ exists and is integrable on $[0,2]$. We will show that $f(x)$ is continuous on $(0,2]$, and has a discontinuity at $x = 0$, where $f(0) = 0$. This is sufficient because a function $f\colon[0,2]\to\mathbf R$ is Riemann integrable if and only if it is bounded and continuous almost everywhere. Finally, by "correcting" the value of $f$ at $x=0$ to produce a continuous function $\tilde f$ on $[0,2]$, we can evaluate the integral $\int_0^2f(x)\,\mathrm dx$ by applying the fundamental theorem of calculus to integrate $\int_0^2\tilde f(x)\,\mathrm dx$.
For each $x\in (0,2]$, put $x_n = f_n(x)$. We claim that $x_n\to \frac{1}{2} + \frac{1}{2}\sqrt{4x+1}$. First we will justify the convergence of the sequence $(x_n)$ using the monotone convergence theorem and then use algebraic limit laws to deduce its limit rigorously.
Consider the sequence $b_{n+1} = \sqrt{2 + b_n}$, where $b_1 = \sqrt 2$. Clearly $b_1 \leqslant 2$. Suppose that $b_n\leqslant 2$; then $\sqrt{2 + b_n} \leqslant \sqrt{2 + 2} = 2$. By induction, $b_n \leqslant 2$ for every $n$. Thus the sequence $(b_n)$ is bounded above by $2$. The sequence $(b_n)$ is clearly monotone increasing. Hence $(b_n)$ converges.
Since $(x_n)$ is monotone and bounded above by $\lim_{n\to\infty}b_n$, it converges, so set $L = \lim_{n\to\infty}x_n$. Since $(x_n)$ converges to $L$, every subsequence converges to $L$, so by continuity of the map $y\mapsto y^2$, we have
$$
L^2 = \big(\lim_{n\to\infty}x_{n+1}\big)^2 = \lim_{n\to\infty}x_{n+1}^2 = \lim_{n\to\infty}x + x_n = x + L.
$$
Therefore $L$ is a root of the polynomial $y\mapsto y^2 - y - x$, and we can conclude that
$$
L = \frac{1}{2} + \frac{1}{2}\sqrt{4x+1}.
$$
(Note that $L$ could not be the negative root since $x_n \geqslant 0$ for every $n\geqslant 1$ and every $x\in(0,2]$.)
If $x = 0$, then $x_n = 0$ for each $n$, so all told,
$$
f(x) = \begin{cases}
\frac{1}{2} + \frac{1}{2}\sqrt{4x+1}, &\text{if $x\in (0,2]$,}\\
0, &\text{if $x = 0$.}
\end{cases}
$$
Thus we have proved that the integrand $f(x)=\sqrt{x+\sqrt{x+\dotsb}}$ is bounded and continuous almost everywhere on $[0,2]$, so it is Riemann integrable on $[0,2]$. The function $\tilde f\colon [0,2]\to\mathbf R$ defined by
$$
\tilde f(x) = \frac{1}{2} + \frac{1}{2}\sqrt{4x+1},\quad\text{for each $x\in[0,2]$,}
$$
is continuous on $[0,2]$ and agrees with $f$ almost everywhere, so $\int_0^2 f(x)\,\mathrm dx = \int_0^2\tilde f(x)\,\mathrm dx$. Since $\tilde f(x)$ is continuous on the entire closed interval $[0,2]$, its integral can be evaluated via the fundamental theorem of calculus. Since
$$
F(x) = \frac{x}{2} + \frac{(4x+1)^{3/2}}{12}
$$
satisfies $F'(x) = \tilde f(x)$ for each $x\in [0,2]$, by the fundamental theorem of calculus,
\begin{align*}
\int_0^2\tilde f(x)\,\mathrm dx = F(2)-F(0) = \bigg[\frac{2}{2} + \frac{(4\cdot 2+1)^{3/2}}{12}\bigg] - \bigg[\frac{0}{2}+\frac{(4\cdot0+1)^{3/2}}{12}\bigg]=\frac{19}{6},
\end{align*}
as desired.
Best Answer
The form of the integrand suggests writing $$(1 + \log x)x^x = (1+\log x)e^{x \log x},$$ then observing that by the product rule, $$\frac{d}{dx}\bigl[x \log x\bigr] = x \cdot \frac{1}{x} + 1 \cdot \log x = 1 + \log x.$$ Consequently, the integrand is of the form $f'(x) e^{f(x)}$, and its antiderivative is simply $$e^{f(x)} = e^{x \log x} = x^x.$$