Another integral that Mathematica cannot do

The integral to evaluate is

$$\int_0^{\infty} dx \frac{\sin{\left (\pi x^2 \right )}}{\sinh^2{(\pi x)}} \cosh{(\pi x)}$$

Given the trig functions in the integrand, it makes sense to use the residue theorem based on a complex integral around a rectangular contour. As has been my experience with these integrals, the integrand of the complex integral will not match the integrand of the real integral to be evaluated.

I am asked frequently how I determine the integrand of the complex integral when I am performing a residue theorem evaluation of a real integral. Of course, first one must determine a contour of integration in the complex plane. In general, this can be difficult but for problems involving trig functions or periodic functions in general, a rectangular contour is used. That said, once the contour is determined, if I do not know any better I experiment with the integrand. So perhaps I evaluate the integral around the contour using the integrand of the real integral to be evaluated. In most cases, that will not produce the real integral, but it will provide insight into what the integrand should be.

I encourage the reader to try integrating the original integrand about the rectangular contour described below to see what happens.

Without further ado, consider the complex integral

$$\oint_C dz \frac{\cos{\left (\pi z^2\right )}}{\sinh^3{\left (\pi z\right )}} $$

about the rectangle with vertices $\pm R \pm i$ with small semicircular detours around the poles at $z=\pm i$. The contour integral is then equal to

$$PV \int_{-R}^R dx \frac{\cos{[\pi (x-i)^2]}}{\sinh^3{[\pi (x-i)]}} + i \epsilon \int_{\pi}^0 d\phi \, e^{i \phi} \frac{\cos{\left [\pi \left (-i+\epsilon e^{i \phi} \right )^2 \right ]}}{\sinh^3{\left [\pi \left (-i+\epsilon e^{i \phi} \right ) \right ]}} \\ + PV \int_R^{-R} dx \frac{\cos{[\pi (x+i)^2]}}{\sinh^3{[\pi (x+i)]}} + i \epsilon \int_{2 \pi}^{\pi} d\phi \, e^{i \phi} \frac{\cos{\left [\pi \left (i+\epsilon e^{i \phi} \right )^2 \right ]}}{\sinh^3{\left [\pi \left (i+\epsilon e^{i \phi} \right ) \right ]}} \\ + i \int_{-1}^1 dy \frac{\cos{\left [\pi \left (R+i y \right )^2 \right ]}}{\sinh^3{\left [\pi \left (R+i y \right ) \right ]}}+i \int_1^{-1} dy \frac{\cos{\left [\pi \left (-R+i y \right )^2 \right ]}}{\sinh^3{\left [\pi \left (-R+i y \right ) \right ]}}$$

Note that the first and third integrals are actually expressed as Cauchy principal values because the individual integrals themselves do not converge. That said, when combined, the resulting integral does converge and we may remove the $PV$ label.

As $R \to \infty$, the last two integrals go to zero.

The second integral approaches, in the limit as $\epsilon \to 0$:

$$i \epsilon \int_{\pi}^0 d\phi \, e^{i \phi} \frac{-1-2 \pi^2 \epsilon^2 e^{i 2 \phi}+\cdots}{\left (\pi \epsilon e^{i \phi} \right )^3 \left (1 + \frac16 \pi^2 \epsilon^2 e^{i 2 \phi} + \cdots \right )^3} \\ = i \epsilon \int_{\pi}^0 d\phi \, e^{i \phi} \frac{-1-2 \pi^2 \epsilon^2 e^{i 2 \phi}}{\left (\pi \epsilon e^{i \phi} \right )^3} \left (1-\frac12 \pi^2 \epsilon^2 e^{i 2 \phi} \right ) \to -i \frac32$$

The fourth integral approaches an identical limit as $\epsilon \to 0$.

The first and third integrals combine to produce, as $R \to \infty$,

$$i 4 \int_{-\infty}^{\infty} dx \frac{\sin{\left ( \pi x^2 \right )}}{\sinh^2{(\pi x)}} \cosh{(\pi x)} $$

The contour integral is also equal to $i 2 \pi$ times the residue of the integrand at $z=0$. The residue may be computed by expanding the integrand in a Laurent series about $z=0$, which is

$$\frac1{(\pi z)^3} \left (1 – \frac12 \pi^2 z^4+\cdots \right ) \left (1 – \frac12 \pi^2 z^2+\cdots \right )$$

The residue is the coefficient of $z^{-1}$, or $-1/(2 \pi)$. Thus,

$$i 4 \int_{-\infty}^{\infty} dx \frac{\sin{\left ( \pi x^2 \right )}}{\sinh^2{(\pi x)}} \cosh{(\pi x)} – i 3 = -i$$

Rearranging things a bit, we find that the original integral is

$$\int_0^{\infty} dx \frac{\sin{\left ( \pi x^2 \right )}}{\sinh^2{(\pi x)}} \cosh{(\pi x)} = \frac14 $$

Computing an integral over an absolute value using Cauchy’s theorem

The problem is to compute the following integral: $$\int_{-1}^1 dx \frac{|x-y|^{\alpha}}{(1-x^2)^{(1+\alpha)/2}}$$ I will show how to compute this integral using Cauchy’s theorem. It was remarked that it should not be possible to use Cauchy’s theorem, as Cauchy’s theorem only applies to analytic functions, and an absolute value certainly does not qualify. True. Nevertheless, for the […]

Very nifty limit

Problem: Find the value of the limit $$\lim_{n \to \infty}n\left(\left(\int_0^1 \frac{1}{1+x^n}\,\mathrm{d}x\right)^n-\frac{1}{2}\right)$$ Solution: I chose this problem because the answer is highly nontrivial and just out of left field. But the process of getting there seems so straightforward; it is not really. Substitute $x=u^{1/n}$ in the integral and get $$I(n) = \int_0^1 \frac{dx}{1+x^n} = \frac1n \int_0^1 […]

Integral of polynomial times rational function of trig function over multiple periods

Problem: Compute $$\int_0^{12\pi} dx \frac{x}{6+\cos 8x}$$ It looks like the OP is trying to compute the antiderivative and use the fundamental theorem of calculus. With multiple periods, that approach is paved with all sorts of difficulty that is really an artifice related to the functional form of the antiderivative. In truth, there should be no […]

A method of evaluating a double integral that nobody taught you in school

Many times we are given integrals to evaluate. The standard way to evaluate is to find a series of transformations that will render the integral into something we know how to evaluate and then proceed. Examples of such transformations are substitutions, parts, replacement of an integrand with another integral, reversing order of integration, and so […]

Deceptively easy product

The problem is to evaluate $$\prod_{n=2}^{\infty} \left (1+\frac{(-1)^{n-1}}{a_n} \right ) $$ where $$a_n = n! \sum_{k=1}^{n-1} \frac{(-1)^{k-1}}{k!} $$ This is one of those cases where trying out a bunch of numbers really isn’t going to help all that much. The $a_n$ look kind of like $n!$, except off by some. Even so, I can find […]

A sum I can only imagine being evaluated using the residue theorem

The challenge this time is to evaluate the following sum: $$\sum_{n=1}^{\infty} \frac{1}{n^3\sin{\left (\sqrt{2} \, n\pi \right)}}$$ NB (20 Nov 2016) I just got word that Wolfram fixed the problem described below. See below for details. It should be clear that the sum converges…right? No? Then how do we show this? Numerical experiments are more or […]

A little perspective on reaching 100K Rep on Math.SE

I stumbled onto on December 16, 2012. Before that, I was occasionally on the prowl for Putnam exam prep questions, various university math contests, and stuff like the IBM monthly puzzle.  At this point 1132 days later, I forget what I was looking for in the first place – likely a new problem to […]

Systematic treatment of a deceptively messy Cauchy principal value integral

The problem here is to evaluate the following: $$PV \int_0^{\infty} dx \frac{\log^2{x}}{(x-1)^2(x-4) \sqrt{x}} $$ This can be done using complex analysis, but it is a pretty involved affair, deceptively so. After struggling with the problem of how to present the solution, I am going to lay out a systematic approach that ignores the motivation behind […]

Cauchy principal value of a convolution

The problem here is to compute the following convolution-type integral: $$\int_{-1}^1 dx \frac{\sqrt{1-x^2}}{\lambda-x}$$ When $-1 \lt \lambda \lt 1$, this integral is infinite, but its Cauchy principal value may be defined. This integral is interesting because of the branch points. And so, away we go… Consider the following contour integral: $$\oint_C dz \frac{\sqrt{z^2-1}}{\lambda-z} $$ where […]

Computing the Convolution of Two Pulses: Graphical vs Analytical

Recently, a user on Math.SE presented a problem of computing the convolution of two pulses: a triangular pulse (impulse response) $$h(t) = \begin{cases} t & 0 \lt t \lt 2 T \\ 0 & t \lt 0 \cup t \gt 2 T\end{cases} $$ and a rectangular pulse (input) $$x(t) = \begin{cases} 1 & 0 \lt […]

Integral with two branch cuts II

The problem here is to compute $$\int_0^\infty \log(1+tx)t^{-p-1}dt$$ where $p\in(0,1)$ and $x>0$. This is a great problem for contour integration. Just tricky enough to be really interesting. What makes it interesting is that there are two functions in the integrand needing their own separate branch cuts. One must keep in mind that each function only […]

Integral with two branch cuts

The problem is to evaluate the following integral: $$\int_{-1}^1 dx \frac{\log{(x+a)}}{(x+b) \sqrt{1-x^2}} $$ where $a \gt 1$ and $|b| \lt 1$. It should be obvious to those who spend time around these integrals that this integral does not converge as stated. However, we only have a simple pole at $x=-b$ so that we can compute […]

Fascinating Fourier Transform

Sometimes I come across a Fourier integral that I have no idea how to attack. And then I find that I can convert it to another, more familiar integral using complex analysis. So here’s an example of such a satisfying situation. The problem is to evaluate $$\int_{-\infty}^{\infty} dx \, (1+i a^2 x)^{-1/2} (1+i b^2 x)^{-1/2} […]

Cauchy Principal Value

On Math.SE, we frequently get integrals that simply do not exist in the usual sense because the integration path intersects a pole of the integrand. Occasionally, we come across such integrals in the course of evaluating integrals of functions with removable singularities using complex methods. However, sometimes such integrals are interesting in their own right. […]

Inverse Laplace Transforms and Delta Functions

Problem: find the inverse Laplace transform of $$F(s) = \frac{s}{s-1}$$ Solution: Well, this one should be easy and one wonders why we are even bothering.  Just split $F$ up as follows: $$F(s) = \frac1{s-1} + 1$$ The ILT of $1/(s-1)$ is simply $e^t$, and the ILT of $1$ is $\delta(t)$.  Done. Or are we?  Why […]

Adventures in integration, University Edition

I got a request from an old friend whom I didn’t even know existed. He is the son of my grand-advisor, if that makes any sense. He teaches, among others, a course in Real Analysis at a university in Australia. I know I must like this guy because he says stuff like this: Currently, within the School […]

Inverse Laplace Transform with a coinciding pole and branch point

Recently, the following Laplace transform was asked to be inverted: $$F(s) = s^{-a-1} e^{-s^a} $$ where $a \in (0,1)$. This is a tough problem for two reasons. One is that there is very little chance of there being an analytical result for arbitrary values of $a$. The other, however, is more subtle: there is a […]

The art of using the Residue Theorem in evaluating definite integrals

As many of you know, using the Residue Theorem to evaluate a definite integral involves not only choosing a contour over which to integrate a function, but also choosing a function as the integrand. Many times, this is an easy task when integrating, say, rational functions over the real line. Sometimes there are less trivial […]

A Tale of Two Integrals

Two integrals look almost the same, and even to those fairly well-versed in the art, are the same. But alas, as we shall see. Consider the integral $$I_1 = \int_0^{\pi} dx \frac{x \sin{x}}{1+\cos^2{x}} $$ This may be evaluated by subbing $x \mapsto \pi-x$ as follows: $$\begin{align} I_1 &= \int_0^{\pi} dx \frac{(\pi-x) \sin{x}}{1+\cos^2{x}} \\&= \pi \int_0^{\pi} […]

Weird integral whose simplicity is only apparent in the complex plane

Every so often we come across an integral that seems absolutely impossible on its face, but is easily attacked – in fact, is custom designed – for the residue theorem. I wonder why a first year complex analysis class doesn’t show off this case as Exhibit A in why the residue theorem is so useful. […]

Expansions of $e^x$

A very basic question was asked recently: What is a better approximation to $e^x$, the usual Taylor approximation, or a similar approximation involving $1/e^{-x}$? More precisely, given an integer $m$, which is a better approximation to $e^x$: $$f_1(x) = \sum_{k=0}^m \frac{x^k}{k!} $$ or $$f_2(x) = \frac1{\displaystyle \sum_{k=0}^m \frac{(-1)^k x^k}{k!}} $$ The answer is amazingly simple: […]