En busca de la integral:
Con $b>0, L>0, \alpha_0>b,\sigma>0, x>L$ ,
$$\phi(x;\alpha_0,\sigma_)=\int_b^\infty\alpha L^{\alpha } x^{-\alpha -1} \frac{e^{-\frac{\left(\log (\alpha -b)-\log (\alpha_0-b)+\frac{\sigma ^2}{2}\right)^2}{2 \sigma ^2}}}{\sqrt{2 \pi } \sigma (\alpha -b)}d\alpha$$
$\textbf{Background}$: Esta es la densidad de una distribución de Pareto $\alpha L^{\alpha } x^{-\alpha -1} $ con su cola exponente $\alpha$ distribuidos de acuerdo a un desplazado logarítmico-normal con una media de $\alpha_0$ y el límite inferior $b$. En otras palabras $\alpha + b$ sigue una $\text{LogNormal}\left(\log (\alpha_0-b)-\frac{\sigma ^2}{2},\sigma \right)$.