Limit matters in improper integral

Recently, one of my best friends, Joking, asked me to solve an integral \begin{equation}\int_0^1 \frac{x-1}{\log x}dx\,.\end{equation}

His intention is to show me a trick as in this video: Let $I(b)\equiv \int_0^1 \frac{x^b-1}{\log x}dx$. By taking derivative to $b$, we obtain an simple differential equation $I'(b)=\frac{1}{b+1}$. With $I(0)=0$, we solve $I(b)$ and thus have $I(1)=\log 2$.

However, my first attempt gave the answer of zero: Let $y=x^2$, $\int_0^1 \frac{x}{\log x}dx=\int_0^1 \frac{dx^2}{\log x^2}=\int_0^1 \frac{dy}{\log y}$. As a result, $\int_0^1 \frac{x-1}{\log x}dx=\int_0^1 \frac{x}{\log x}dx-\int_0^1 \frac{dx}{\log x}=\int_0^1 \frac{dy}{\log y}-\int_0^1 \frac{dx}{\log x}=0$. Clearly, the integral cannot be zero. So what's wrong here?

The reason is that this is an improper integral. Strictly speaking, we should solve \begin{equation} \lim_{\epsilon\rightarrow 0} \int_0^{1-\epsilon} \frac{x-1}{\log x}dx\,.\end{equation} The limit here matters. With this limit, we end with \begin{equation} \int_0^{1-\epsilon} \frac{x-1}{\log x}dx= \int_0^{(1-\epsilon)^2} \frac{dy}{\log y}-\int_0^{(1-\epsilon)} \frac{dx}{\log x}=\int_{(1-\epsilon)}^{(1-\epsilon)^2}\frac{dx}{\log x}\,,\tag{1}\end{equation} which is not necessarily zero.

To solve (1), we define $f(x)\equiv \int_{(1-\epsilon)}^x \frac{dt}{\log t}$. The question is to solve the value $f((1-\epsilon)^2)$. By Taylor expansion at $1-\epsilon$, we have
\begin{equation} f((1-\epsilon)^2)=\sum_{k=0}^{\infty} \frac{1}{k!} f^{(k)}(1-\epsilon)\left[(1-\epsilon)^2-(1-\epsilon)\right]^{k} = \sum_{k=0}^2 \frac{1}{k!} f^{(k)}(1-\epsilon) (-\epsilon)^k\,. \end{equation}
Because of the existence of $\epsilon^k$, we only have to compute $f^{(k)}(1-\epsilon)$ to the leading order $1/\epsilon^k$. Keeping this in mind, we can compute $f'(x)=\frac{1}{\log x}$, $f''(x)=-\frac{1}{(\log x)^2 x}$, $f^{(3)}(x)\approx\frac{2}{(\log x)^3 x^2}$, etc. So $f^{(k)}(x)\approx-\frac{(k-1)!}{(-\log x)^k x^{k-1}}$, and thus $f^{(k)}(1-\epsilon)\approx-\frac{(k-1)!}{(-\log(1-\epsilon))^k}=-\frac{(k-1)!}{\epsilon^k}$. Put it back to the Taylor expansion, we finally have \begin{equation} f((1-\epsilon)^2)=-\sum_{k=0}^{\infty}\frac{1}{k!} \frac{(k-1)!}{\epsilon^k} (-\epsilon)^k=-\sum_{k=0}^{\infty}\frac{(-1)^k}{k}=\log 2\,,\end{equation} where we used the series $\sum_{k=0}^{\infty}{x^k}/{k}=-\log(1-x)$ in the last equal sign. The final answer is the same as that from the trick. Cheers!

Comments

Post a Comment

Popular posts from this blog

529 Plan

How to offset W2 tax

Retirement Accounts