Por ahora ignoramos la justificación para cambiar el orden de integración y resolver un problema más fácil con la ayuda de Mathematica. La región de integración es el triángulo delimitado por la $x$ $z$ ejes y la línea de $x + z = 1$. Cuando se cambie el orden de integración tenemos la integral
$$J(a,b) = \int_0^1 \int_0^{1-z} \frac{dx dz}{az(z-1)+axz+x(1-b)}$$
El denominador es $0$ siempre
$$x_0 = \frac{az(1-z)}{az+1-b}$$
así que tenemos que calcular una integral impropia de primera. Tenemos
$$\int_0^{x_0-\epsilon} \frac{dx}{az(z-1)+axz+x(1-b)} = \frac{\log((-1+b-az)\epsilon)-\log(az(z-1))}{az+1-b}$$
y como esto no es finito como $\epsilon$$0$, tenemos que admitir que la integral impropia no converge. Del mismo modo
$$\int_{x_0+\epsilon}^{1-z} \frac{dx}{az(z-1)+axz+x(1-b)} = \frac{\log((1-b)(1-z))-\log((1-b+az)\epsilon)}{az+1-b}$$
no es finito como $\epsilon$$0$. Sin embargo, podemos calcular el valor del capital por poner las dos integrales juntos y dejar que el límite de enfoque de integración de $x_0$ a la misma velocidad en ambos lados
$$P.V.\int_0^{1-z} = \frac{\log((1-b)(1-z))-\log(az(1-z))}{1-b+az}$$
Integrando de nuevo con respecto a z tenemos
$$P.V.J(a,b) = \int_0^1\frac{\log((1-b)(1-z))-\log(az(1-z))}{1-b+az} dz = \frac{1}{a}\left(\log\left(\frac{a}{1-b}\right)\log\left(\frac{1-b}{1+a-b}\right)-Li_2\left(-\frac{a}{1-b}\right)\right)$$
Podemos ampliar esta como
$$-\frac{1}{a}(\log(a)\log(1+a)+Li_2(-a))-\frac{\log(a)}{1+a}b-\frac{1+a+(2+a)\log(a)}{2(1+a)^2}b^2+O(b^3)$$
Al $b$ es pequeña, nos quedamos con
$$-\frac{1}{a}(\log(a)\log(1+a)+Li_2(-a))$$
La pregunta es - fuimos justificados en intercambiar el orden de integración en el primer lugar? No estoy seguro de si esto le ayuda.