Primero, hacer cumplir la sustitución$x\mapsto x^{1/n}$ revela
$$ \begin{align}
\int_0^1 \frac{f(ax^n)-f(bx^n)}{x}\,dx=\frac1n \int_0^1 \frac{f(ax)-f(bx)}{x}\,dx
\end {align} $$
En segundo lugar, es sencillo mostrar que
$$ \begin{align}
\int_\epsilon^1 \frac{f(ax)}{x}\,dx-\int_\epsilon^1 \frac{f(bx)}{x}\,dx&=\int_{a\epsilon}^a \frac{f(x)}{x}\,dx-\int_{b\epsilon}^b \frac{f(x)}{x}\,dx\\\\
&=\int_{a\epsilon}^{b\epsilon}\frac{f(x)}{x}\,dx-\int_a^b \frac{f(x)}{x}\,dx\\\\
&=f(0)\log(b/a) +\int_{a\epsilon}^{b\epsilon}\frac{f(x)-f(0)}{x}\,dx-\int_a^b \frac{f(x)}{x}\,dx\\\\
\end {align} $$
En tercer lugar, tenemos la estimación
ps
que se acerca a$$\left|\int_{a\epsilon}^{b\epsilon}\frac{f(x)-f(0)}{x}\,dx\right|\le \int_{a\epsilon}^{b\epsilon} \frac{|f(x)-f(0)|}{x}\,dx\le \sup_{x\in[a\epsilon,b\epsilon]}|f(x)-f(0)| \log(b/a)$ como$0$ ya que$\epsilon \to 0$ es continuo.
Finalmente, armarlo revela
ps
como se iba a mostrar!