Figure 1 shows a graphic in a tweet by Cliff Pickover that caught my attention recently:
 |
| Figure 1 |
It turns out that the denominator of \(x^2+1)^2\) gives the same result for \(x^2+1\) so in this video I'll just be considering the integral:$$\text{I}=\int_{-\infty} ^{\infty} \frac{\cos x}{x^2+1} \text{ d}x$$It's such an aesthetic result that I thought I'd investigate how the result was arrived at. Flammable Maths YouTube channel solves this definite integral in three different ways, although in this post I'll only be looking at the Feynman method. I'll be simply reproducing the steps outlined in the video that is shown below. It's a good opportunity to practise my LaTeX.
The first step is to introduce a variable \(t\) into the integral:$$\text{I}(t)=\int_{-\infty} ^{\infty} \frac{\cos (t \, x)}{x^2+1} \text{ d}x$$When \(t=1\), we have \(\text{I}(1)=\text{I}\). Now we'll determine the first and second derivatives. Let's start with the first derivative:$$\text{I'}(t)=\int_{-\infty} ^{\infty} \partial_t \,\frac{\cos (t \, x)}{x^2+1} \text{ d}x=\int_{-\infty} ^{\infty} \frac{1}{x^2+1} \, -x \sin (t \, x) \text{ d}x$$Now we multiply top and bottom by \(x\) and add \(0\) in the form of \(+1 + -1\):$$\text{I'}(t)=-\int_{-\infty} ^{\infty} \frac{(x^2+1-1) \, \sin (t \, x)}{x \,(x^2+1)} \text{ d}x$$We can now split the integral into two parts thanks to our \(+1 + -1\) trick:$$\text{I'}(t)=-\int_{-\infty} ^{\infty} \frac{(x^2+1) \, \sin (t \, x)}{x \,(x^2+1)} \text{ d}x+\int_{-\infty} ^{\infty} \frac{ \sin (t \, x)}{x \,(x^2+1)} \text{ d}x$$The first part of the integral simplifies to:$$-\int_{-\infty} ^{\infty} \frac{ \sin (t \, x)}{x} \text{ d}x=-\pi $$We'll just accept that result for the moment and so we have the following result for the first derivative:$$\text{I'}(t)=-\pi+\int_{-\infty} ^{\infty} \frac{ \sin (t \, x)}{x \,(x^2+1)} \text{ d}x$$Now we'll find the second derivative:$$\text{I''}(t)=\int_{-\infty} ^{\infty} \partial_t \frac{ \sin (t \, x)}{x \,(x^2+1)} \text{ d}x=\int_{-\infty} ^{\infty} \frac{ \cos (t \, x)}{(x^2+1} \text{ d}x=\text{I}(t) $$The fact the the second derivative of the function is equal to the original function means that we have a second order linear differential equation of the form:$$\text{I''}(t)-\text{I}(t)=0 \text{ where we'll assume that I}(t)=c \, e^{\, \lambda \,t}$$This means that \(\text{I''}(t)=\lambda^2 \,c\,e^{\, \lambda \, t}\) and substituting this into the differential equation we get:$$\lambda^2 \,c\,e^{\, \lambda \, t}-c \, e^{\, \lambda \,t}=0 \implies c \, e^{\, \lambda \,t}(\lambda^2-1)=0 \text{ and }\lambda=\pm1$$Thus we have: $$ \text{I}(t)=c_1 \,e^{\, t}+c_2 \,e^{-t} \text{ and }\text{I'}(t)=c_1 \,e^{\, t}-c_2 \,e^{-t}$$$$\text{I}(0)=\pi=c_1+c_2 \text{ and } \text{I'}(0)=-\pi =c_1-c_2 \implies c_1=0 \text{ and }c_2=\pi$$ $$ \text{I}(t)=\pi \,e^{-t} \text{ and thus }\text{I}(1)=\text{I} =\frac{\pi}{e}$$.Figure 2 shows a graph of the function:
 |
| Figure 2 |
Clearly there are negative areas under the curve that will cancel with the positive areas but the bulk of the area lies between \( \frac{-\pi}{2} \) and \(\frac{\pi}{2}\).
Here are some links that Flammable Maths provides for some of the techniques used in this video:
There are lots of interesting integrals like the one we have just dealt with and I am certainly out of practice in dealing with them. I should make posts like this more frequently. Looking back over my previous posts, I notice that I've made the following integration-related posts:
- The Wallis Formula for Pi and the Dawson Function (October 22nd 2016): this post deals with the following integral:$$ \int e^{x^2} \, \text{d}x$$
- Novel Integration Technique (September 29th 2017): this post also deals with Feynman's technique as it applies to the following integral:$$ \int_0^{\infty} \frac{x}{\sin x}\, \text{d}x$$
- Average Distance Between Two Points in a Square (April 28th 2019): this entails modifying an integral by changing into polar coordinates so that it is easier to deal with.
