I am trying to solve the following improper integral $I$.
$$I = \int\limits_0^\infty\frac{1}{\left(x+1\right)\left(ax+1\right)\left(bx+1\right)} dx$$
However, I myself get the divergent answer. When I tried it through an online-calculator, it gives me something else. Can someone guide me through the following? I will really appreciate.
\begin{align*} I &=\lim_{x\to \infty}\frac{\left(a-1\right)b\ln\left(\left|bx+1\right|\right)+\left(a-ab\right)\ln\left(\left|ax+1\right|\right)+\left(b-a\right)\ln\left(\left|x+1\right|\right)}{\left(a-1\right)\left(b-1\right)\left(b-a\right)} \\ &\quad - \lim_{x\to 0}\dfrac{\left(a-1\right)b\ln\left(\left|bx+1\right|\right)+\left(a-ab\right)\ln\left(\left|ax+1\right|\right)+\left(b-a\right)\ln\left(\left|x+1\right|\right)}{\left(a-1\right)\left(b-1\right)\left(b-a\right)}. \end{align*}
While solving the part for $\lim_{x\to0}$, I get $0$ as the output, but if I solve $\lim_{x\to\infty}$, I get a divergent term. While in actual if I solve it through an online calculator, I get the following answer.
$$I = \dfrac{\left(a-1\right)b\ln\left(b\right)-a\ln\left(a\right)b+a\ln\left(a\right)}{\left(a-1\right)\left(b-1\right)\left(b-a\right)}$$
$\newcommand{\bbx}[1]{\,\bbox[15px,border:1px groove navy]{\displaystyle{#1}}\,} \newcommand{\braces}[1]{\left\lbrace\,{#1}\,\right\rbrace} \newcommand{\bracks}[1]{\left\lbrack\,{#1}\,\right\rbrack} \newcommand{\dd}{\mathrm{d}} \newcommand{\ds}[1]{\displaystyle{#1}} \newcommand{\expo}[1]{\,\mathrm{e}^{#1}\,} \newcommand{\ic}{\mathrm{i}} \newcommand{\mc}[1]{\mathcal{#1}} \newcommand{\mrm}[1]{\mathrm{#1}} \newcommand{\pars}[1]{\left(\,{#1}\,\right)} \newcommand{\partiald}[3][]{\frac{\partial^{#1} #2}{\partial #3^{#1}}} \newcommand{\root}[2][]{\,\sqrt[#1]{\,{#2}\,}\,} \newcommand{\totald}[3][]{\frac{\mathrm{d}^{#1} #2}{\mathrm{d} #3^{#1}}} \newcommand{\verts}[1]{\left\vert\,{#1}\,\right\vert}$ \begin{align} I & = \int_{0}^{\infty}{\dd x \over \pars{x + 1}\pars{ax + 1}\pars{bx + 1}} \\[5mm] & = \int_{0}^{\infty}\braces{2\int_{0}^{1}\dd u_{1}\int_{0}^{u_{1}}\dd u_{2}\,{1 \over \bracks{x + 1 + u_{1}\pars{a - 1}x + u_{2}\pars{b - a}x}^{\, 3}}}\dd x \\[5mm] & = 2\int_{0}^{1}\dd u_{1}\int_{0}^{u_{1}}\dd u_{2}\int_{0}^{\infty} {\dd x \over \braces{\bracks{\pars{a - 1}u_{1} + \pars{b - a}u_{2} + 1}x + 1}^{\, 3}} \\[5mm] & = \int_{0}^{1}\dd u_{1}\int_{0}^{u_{1}}\dd u_{2}\, {1 \over \pars{a - 1}u_{1} + \pars{b - a}u_{2} + 1} \\[5mm] & = \left.{1 \over b - a}\int_{0}^{1}\dd u_{1}\, \ln\pars{\bracks{a - 1}u_{1} + \bracks{b - a}u_{2} + 1} \,\right\vert_{\ u_{2}\ =\ 0}^{\ u_{2} =\ u_{1}} \\[5mm] & = {1 \over b - a}\int_{0}^{1}\dd u_{1}\braces{\vphantom{\Large A}% \ln\pars{\bracks{b - 1}u_{1} + 1} - \ln\pars{\bracks{a - 1}u_{1} + 1}} \\[5mm] & = {1 \over b - a}\braces{% \bracks{-1 + {b\ln\pars{b} \over b - 1}} - \bracks{-1 + {a\ln\pars{a} \over a - 1}}} \\[5mm] & = \bbx{{\pars{b - 1}a\ln\pars{a} - \pars{a - 1}b\ln\pars{b} \over \pars{a - 1}\pars{a - b}\pars{b - 1}}} \end{align}