Let
Differentiate with respect to a :
Substitute y = √x :
Polynomial division yields
Computing I'(a) isn't so difficult from here. You'd find (assuming a ≥ 0)
Integrate both sides with respect to a. On the right side, substituting b = √a yields
Noting that a = 0 makes the integral I(a) vanish, we have
and so
We recover the integral we want with a = 1, which gives a value of