Tenemos
$$ f(a,b,c) = \int_{0}^{\frac{\pi}{2}} \frac{\arctan(a \tan^2\theta)}{b \sin^2\theta + c\cos\theta^2} \, d\theta = \frac{\pi}{\sqrt{bc}} \left( \arctan\left(1 + \sqrt{\frac{2ac}{b}}\right) - \frac{\pi}{4} \right). $$
Aquí es un boceto de la computación.
Paso 1. Reducción del problema (no imprescindible)
Considere la curva $\gamma(\theta) = (\sqrt{c}\cos\theta, \sqrt{b}\sin\theta)$$\theta \in [0,\pi/2]$. Esto es parte de una elipse en el primer cuadrante. Observando que en el $ \sqrt{bc} d\theta = x dy - y dx$, se deduce que
$$ f(a,b,c) = \int_{\gamma} \omega, \quad \text{where} \quad \omega = \frac{1}{\sqrt{bc}} \arctan\left(\frac{ac}{b}\cdot\frac{y^2}{x^2}\right) \frac{x dy - y dx}{x^2 + y^2}. $$
Escribir $g = ac/b$ por la simplicidad. A continuación, la observación crucial es que
$$ d\omega = \frac{1}{\sqrt{bc}} \left( \frac{\partial}{\partial x} \frac{x \arctan(gy^2/x^2)}{x^2 + y^2} + \frac{\partial}{\partial y} \frac{y \arctan(gy^2/x^2)}{x^2 + y^2} \right) dx \wedge dy = 0 $$
en el primer cuadrante. Así que por el Verde del teorema, podemos substituir $\gamma$ por cualquier definidas a trozos $C^1$-curva en el primer cuadrante que se une a de$P = (\sqrt{c}, 0)$$Q = (0, \sqrt{b})$. Puesto que la integral de $\omega$ a lo largo de cualquier segmento de línea en cada eje es cero, podemos también sustituir a $P$ por cualquier punto positivo de la $x$-eje y de la misma manera para $Q$. Poner en conjunto, que puede sustituir a $\gamma$ por el arco circular $(x,y) = (\cos\theta, \sin\theta)$ obtener
$$ f(a,b,c)
= \frac{1}{\sqrt{bc}} \int_{0}^{\frac{\pi}{2}} \arctan(g\tan^2\theta) \, d\theta
= \frac{1}{2\sqrt{bc}} \int_{0}^{\infty} \frac{\arctan(g u)}{\sqrt{u}(1+u)} \, du. $$
Paso 2. Cálculo
Ahora el resto de cálculo se puede hacer con la ayuda de residuos de la computación. En efecto, considerar la última integral como una función de la $g$. A continuación, el residuo de cálculo muestra que
$$ \frac{\partial f}{\partial g}
= \frac{1}{2\sqrt{bc}} \int_{0}^{\infty} \frac{\sqrt{u}}{(1+u)(1+ g^2 u^2)} \, du
= \frac{\pi}{2\sqrt{bc}} \cdot \frac{1}{\sqrt{2 g}(1 + \sqrt{2 g} + g)}. $$
La integración de ambos lados y usando el valor inicial $f|_{g=0} = 0$ da la respuesta deseada.
Observación. Me di cuenta de que podemos omitir el Paso 1 y aplicar técnicas en el Paso 2 directamente a
$$ f(a,b,c) = \frac{1}{2} \int_{0}^{\infty} \frac{\arctan (a u)}{\sqrt{u}(c + bu)} \, du. $$
Aún así, el Paso 1 se ve agradable a mí como se desinfecta parámetros y muestra un poco de la técnica.