Como mencioné anteriormente, usemos análisis complejo. Creo que es más fácil hacer primero la segunda integral. Para empezar, consideremos la integral compleja
$$\oint_Cdz \, \frac{z^{-2/3} (z-1)^{-1/3}}{z^2+1} $$
donde $C$ es el siguiente contorno:
donde el radio del arco circular grande es $R$ y el radio de los arcos circulares pequeños es $\epsilon$.
Para evaluar la integral de contorno, dividimos el contorno en varias piezas y evaluamos cada pieza parametrizando según corresponda. En este caso, la integral de contorno es igual a
$$e^{i \pi} \int_R^{\epsilon} dx \, \frac{e^{-i 2 \pi/3} x^{-2/3} e^{-i \pi/3} (x+1)^{-1/3}}{x^2+1} +i \epsilon \int_{\pi}^0 d\phi \, e^{i \phi} \frac{\epsilon^{-2/3} e^{-i 2 \phi/3} \left ( \epsilon e^{i \phi} -1 \right)^{-1/3}}{\epsilon^2 e^{i 2 \phi}+1} \\ + \int_{\epsilon}^{1-\epsilon} dx \, \frac{x^{-2/3} e^{-i \pi/3} (1-x)^{-1/3}}{x^2+1} + i \epsilon \int_{\pi}^{-\pi} d\phi \, e^{i \phi} \frac{\left (1+\epsilon e^{i \phi} \right)^{-2/3} \epsilon^{-1/3} e^{-i \phi/3}}{\left (1+\epsilon e^{i \phi}\right )^2+1} \\ + \int_{1-\epsilon}^{\epsilon} dx \, \frac{x^{-2/3} e^{i \pi/3} (1-x)^{-1/3}}{x^2+1} + i \epsilon \int_0^{-\pi} d\phi \, e^{i \phi} \frac{\epsilon^{-2/3} e^{-i 2 \phi/3} \left ( \epsilon e^{i \phi} -1 \right)^{-1/3}}{\epsilon^2 e^{i 2 \phi}+1} \\ +e^{-i \pi} \int_{\epsilon}^R dx \, \frac{e^{i 2 \pi/3} x^{-2/3} e^{i \pi/3} (x+1)^{-1/3}}{x^2+1} + i R \int_{-\pi}^{\pi} d\theta \, e^{i \theta} \frac{R^{-2/3} e^{-i 2 \theta/3} \left ( R e^{i \theta}-1 \right )^{-1/3}}{R^2 e^{i 2 \theta}+1}$$
Consideremos los límites cuando $\epsilon \to 0$ y $R \to \infty$. Independientemente de estos límites, las primeras y séptimas integrales se cancelan. En estos límites, el lector puede verificar que las segundas, cuartas, sextas y octavas integrales tienden a cero. Por lo tanto, nos quedan las terceras y quintas integrales. La integral de contorno en los límites anteriores se convierte entonces en
$$-i 2 \sin{\frac{\pi}{3}} \int_0^1 dx \, \frac{x^{-2/3} (1-x)^{-1/3}}{1+x^2}$$
Por el teorema de los residuos, la integral de contorno también es igual a $i 2 \pi$ veces la suma de los residuos en los polos $z=\pm i = e^{\pm i\pi/2}$. Tenga en cuenta que la elección del argumento de los polos no es realmente una elección, ya que ya hemos definido esta elección con nuestro corte de rama. En consecuencia, la integral de contorno también es igual a
$$i 2 \pi \left [\frac{e^{-i \pi/3} \left (e^{i \pi/2}-1 \right )^{-1/3} }{i 2} + \frac{e^{i \pi/3} \left (e^{-i \pi/2}-1 \right )^{-1/3} }{-i 2} \right ] $$
Con el corte de rama que hemos definido, tenga en cuenta que $e^{i \pi/2}-1 = 2^{1/2} e^{i 3 \pi/4}$ y $e^{-i \pi/2}-1 = 2^{1/2} e^{-i 3 \pi/4}$. Haciendo el resto de la aritmética y estableciendo el resultado igual a la integral de contorno anterior, concluimos que
$$\int_0^1 dx \, \frac{x^{-2/3} (1-x)^{-1/3}}{1+x^2} = 2^{-1/6} \pi \frac{\sin{\frac{5 \pi}{12}}}{\sin{\frac{\pi}{3}}} = \frac{\sqrt{2+\sqrt{3}}}{2^{1/6} \sqrt{3}} \pi$$
Dejo como ejercicio para el lector demostrar que el resultado es igual al producido por Wolfram Alpha/Mathematica como lo muestra Robert Israel.