It's a very Challenging question perhaps a kind of olympiad question :
Prove that :
$$\int_{0}^{\sqrt{5}}x^{x}dx>4$$
It's pretty sharp since the left hand side is almost $4.0005$ .
I recall the Sophomore dream :
$$\int_{0}^{1}x^{x}dx=-\sum_{n=1}^{\infty}\left(-n\right)^{-n}$$
For the other part we have :
$$\int_{1}^{\sqrt{5}}x^{x}dx\geq \int_{1}^{\sqrt{5}}\left(1-\sqrt{x}+x\right)^{1+x}dx$$
Because I show in another question the inequality $x>0$:
$$\left(1-\sqrt{x}+x\right)^{1+x}\leq x^x$$
But we are far from the goal .
We have numerically :
$$\int_{1}^{1.5}x^{x}dx>\int_{1}^{1.5}\left(e^{\left(-1+\frac{\left(1+x+x^{2}\right)}{3}\right)}+\frac{1}{6}\left(x-1\right)^{2}+\frac{1}{26}\left(x-1\right)^{4}\right)dx$$
How to show it by hand or make a proof where the final steps is the calculus of some constant ?
$\color{green}{\textbf{Splitting.}}$
Let $\;L=\dfrac{582651}{582601},\quad q=\dfrac{44}{43},\;$ then $\;L<q^{-35}\sqrt 5,$ $$I>\int\limits_0^{\sqrt5} x^x\,\text{dx} >\int\limits_0^1 x^x\,\text dx+\int\limits_1^L x^x\,\text dx+\sum_{k=0}^{34}\int\limits_{q^{k-35}\sqrt5}^{q^{k-34}\sqrt5} x^x\,\text dx.\tag1$$
$\color{green}{\textbf{First integral.}}$
In accordance with OP, $$I_0=\int_0^1 x^x\,\text dx = \sum_{k=1}^\infty (-1)^{k+1} k^{-k},\tag{$\Diamond$}$$ so $$I_0 > \sum_{k=1}^8 (-1)^{k+1} k^{-k} >\dfrac1{8^8}\sum_{k=1}^8 \genfrac\lfloor\rfloor{}{}{(-1)^{k+1}8^8}{k^k}=\dfrac1{16777216}(16777216-4194304+621378-262144+5368-359+20-1),$$ $$I_0=\dfrac{13143781}{16777216}>0.783430.\tag2$$
$\color{green}{\textbf{Second interal.}}$
$$I_1=\int\limits_1^L x^x\,\text dx = \int\limits_0^{L-1} (1+y)^{1+y}\,\text dy \ge \int_0^{L-1} (1+y(1+y)) \,\text dy = \dfrac{(L-1)(6+3(L-1)+2(L-1)^2)}6,$$ $$I_1=\dfrac{50(6\cdot 582601^2+3\cdot 50\cdot 582601+2\cdot 50^2)}{6\cdot582601^3} > 0.0000858.\tag3$$
$\color{green}{\textbf{Integrals under the sum.}}$
Let $x\in\big(a,a+\delta a\big),\;\delta=q-1=\dfrac1{43}\ll1,\;$ then $$x=a(1+\delta z)=a +\delta a z,\qquad \big(z\in(0,1),\;a\in(1,\infty)\big),$$ $$x^{x}=a^{\large a+\delta a z}\;\left(1+\delta z \right)^{\large a(1+\delta z)},$$
Taking in account the inequality $$(1+t)^{\large a(1+t)}- e^{at} \ge \dfrac a2 t^2\left(1+\dfrac{3a-1}3 t\right)\tag{4}$$ and identities
$$\int\limits_0^d t^k a^{a t}\,\text dt = \dfrac{\Gamma(1+k)-\Gamma(1+k,-a d\ln(a))}{(-a \ln(a))^{1+k}},$$ $$\dfrac{\Gamma(k+1)-\Gamma(k+1,-y)}{(-y)^{k+1}} =\dfrac1{k+1}+\dfrac y{k+2}+\dfrac{y^2}{2!(k+3)}+\dots,$$ one can get $$\int\limits_0^d t^k a^{a t}\,\text dt = \dfrac{\Gamma(1+k)-\Gamma(1+k,-ad\ln(a))}{(-a d\ln(a))^{k+1}} \dfrac{(ad\ln a)^{k+1}}{(-a \ln(a))^{k+1}},$$ $$\int\limits_0^d t^k a^{a t}\,\text dt =d^{k+1}\left(\dfrac1{k+1}+\dfrac {a d\ln(a)}{k+2}+\dfrac{(ad\ln(a))^2}{2!(k+3)}+\dots\right),\tag5$$ $$\begin{align} &J(a, \delta)=\int\limits_a^{a+\delta a} x^x\,\text dx = \delta a\cdot a^a \int\limits_0^1 a^{\delta az}\left(1+\delta z\right)^{\large a(1+\delta z)}\,\text dz\\[4pt] &= a^{1+a}\int\limits_0^\delta a^{at}\left(1+t\right)^{\large a(1+t)}\,\text dt \ge a^{1+a} \int\limits_0^\delta a^{\large at} \left(e^{at}+\frac a2 t^2 + \frac {a(3a-1)}6 t^3\right)\,\text dt\\[4pt] &=a^{1+a} \int\limits_0^\delta (ea)^{\large a t}\,\text dt + \dfrac12a^{2+a} \int\limits_0^\delta t^2 a^{\large a t}\,\text dt +\dfrac {3a-1}6 a^{2+a} \int\limits_0^\delta t^3 a^{\large a t} \,\text dt,\\[4pt] &J(a,\delta)\ge a^a\dfrac{(ea)^{\large \delta a}-1}{\ln(ea)} +\dfrac12a^{a+2}\delta^3\left(\dfrac13+\dfrac {a\delta\ln(a)}4 +\dfrac{(a\delta\ln(a))^2}{10}+\dots\right)\\[4pt] &+\dfrac{3a-1}6 a^{a+2}\delta^4\left(\dfrac14+\dfrac {a\delta\ln(a)}5+\dfrac{(a\delta\ln(a))^2}{12}+\dots\right), \end{align}$$ $$J\left(a,\dfrac1{43}\right) > a^a\left(\dfrac{(ea)^{\large \frac a{43}}-1}{\ln(ea)}+\dfrac{a^2}{6\cdot 43^3}\left(1+\dfrac{3a-1}{196}\right)\right).\tag6$$ Therefore, $$I_2>\sum_{k=1}^{35} J\left(\left(\dfrac{43}{44}\right)^k\,\sqrt5, \dfrac1{43}\right) > 3.216991.\tag7$$
$\color{green}{\textbf{Final summation.}}$
Taking in account $(1),$ $(2),$ $(3),$ $(7),$ easily to get
$$\color{brown}{\mathbf{\int\limits_0^{\sqrt5} x^x\,\text dx> I_0+I_1+I_2 > 0.783430+ 0.0000858+3.216991 = 4.000506>4.}}$$
Note, that calculations contain only elementary functions and can be executed on the calculator.