I ran across this integral I get no where with. Can someone suggest a method of attack?.

$$\int_0^{\infty}\frac{\sin(\pi x^2)}{\sinh^2 (\pi x)}\mathrm dx=\frac{2-\sqrt{2}}{4}$$

I tried series, imaginary parts, and so forth, but have made no progress.

Thanks very much.

  • $\begingroup$ If it's any consolation, Mathematica can't solve "Integrate[Sin[Pix^2]/Sinh[Pix]^2, {x,0,Infinity}]", which means it's not easy. $\endgroup$ – barrycarter Sep 3 '11 at 16:53
  • $\begingroup$ @barrycarter: Almost two years later, Mathematica still can't get this one. $\endgroup$ – robjohn Jul 30 '13 at 22:04
  • 3
    $\begingroup$ I love it when a brilliant human mind evaluates an integral the fancy math engines can not.:) $\endgroup$ – Cody Jul 31 '13 at 11:15
  • $\begingroup$ @robjohn $\texttt{Mathematica}$ claims that $\texttt{Log[...]}$ assumes a branch-cut along $\left(-\infty,0\right]$. However, it returns $\texttt{Log[-1]} = \pi\,\mathrm{i}$ $\endgroup$ – Felix Marin Apr 6 '17 at 22:21
  • $\begingroup$ @robjohn $\texttt{Mathematica}$ is unable to evaluate ( at least version ) the simple integral $$ \int_{-\infty}^{\infty}{g/\pi \over \left(\omega - \varepsilon\right)^2 + g^{2}}\, {1 \over \mathrm{e}^{\omega/t} + 1}\,\mathrm{d}\omega = {1 \over 2} - {1 \over \pi}\,\Im\Psi\left({1 \over 2} + { g + \varepsilon\,\mathrm{i}\over 2\pi t}\right)\,,\quad g > 0\,,\ \varepsilon \in \mathbb{R}\,,\ t > 0 $$ $\endgroup$ – Felix Marin Apr 6 '17 at 23:37

Although this question is two years old, the integral was mentioned in chat recently, I evaluated it, and then found this question. Since there is no complete solution, although Hans Lundmark's suggestion is excellent and similar in nature, I am posting what I have done.


Since the integrand is even, $$ \begin{align} \int_0^\infty\frac{\sin(\pi x^2)}{\sinh^2(\pi x)}\,\mathrm{d}x &=\frac12\int_{-\infty}^\infty\frac{\sin(\pi x^2)}{\sinh^2(\pi x)}\,\mathrm{d}x \end{align} $$ Define $$ f(z)=\frac{\cos\left(\pi z^2\right)}{\sinh(2\pi z)\sinh^2(\pi z)} $$ Note that because $$ f(x\pm i) =\frac{-\cos\left(\pi x^2\right)\cosh(2\pi x)\pm i\sin\left(\pi x^2\right)\sinh(2\pi x)}{\sinh(2\pi x)\sinh^2(\pi x)}\\ $$ we have $$ \begin{align} \int_\gamma f(z)\,\mathrm{d}z &=\int_{-\infty}^\infty\big[f(x-i)-f(x+i)\big]\,\mathrm{d}x\\ &=-2i\int_{-\infty}^\infty\frac{\sin(\pi x^2)}{\sinh^2(\pi x)}\,\mathrm{d}x\\ &=2\pi i\times\begin{array}{}\text{the sum of the residues}\\\text{inside the contour}\end{array} \end{align} $$ where $\gamma$ is the contour

$\hspace{3.2cm}$enter image description here

Therefore, $$ \int_0^\infty\frac{\sin(\pi x^2)}{\sinh^2(\pi x)}\,\mathrm{d}x =-\frac\pi2\times\begin{array}{}\text{the sum of the residues}\\\text{inside the contour}\end{array} $$ Residues

near $0$ : $$ \begin{align} f(z) &=\frac{\cos\left(\pi z^2\right)}{\sinh(2\pi z)\sinh^2(\pi z)}\\ &=\frac{1-\frac12\pi^2z^4+O(z^8)}{2\pi z\left(1+\frac23\pi^2z^2+O(z^4)\right)\pi^2 z^2\left(1+\frac13\pi^2z^2+O(z^4)\right)}\\ &=\frac{1-\pi^2z^2}{2\pi^3z^3}+O(z)\\[10pt] &\implies\text{residue}=-\frac1{2\pi} \end{align} $$ at $\pm i/2$, use L'Hosptal : $$ \begin{align} \text{residue} &=\lim_{z\to\pm i/2}\frac{(z\mp i/2)\cos\left(\pi z^2\right)}{\sinh(2\pi z)\sinh^2(\pi z)}\\ &=\frac1{2\pi\cosh(\pm\pi i)}\frac{\cos(-\pi/4)}{\sinh^2(\pm\pi i/2)}\\ &=\frac1{2\pi\cos(\pm\pi)}\frac{\sqrt2/2}{-\sin^2(\pm\pi/2)}\\[4pt] &=\frac{\sqrt2}{4\pi} \end{align} $$ near $\pm i$ : $$ \begin{align} f(z\pm i) &=\frac{-\cos\left(\pi z^2\right)\cosh(2\pi z)\pm i\sin\left(\pi z^2\right)\sinh(2\pi z)}{\sinh(2\pi z)\sinh^2(\pi z)}\\ &=\frac{-\left(1-\frac12\pi^2z^4+O(z^8)\right)\left(1+2\pi^2z^2+O(z^4)\right)+O(z^3)}{2\pi z\left(1+\frac23\pi^2z^2+O(z^4)\right)\pi^2 z^2\left(1+\frac13\pi^2z^2+O(z^4)\right)}\\ &=-\frac{1+\pi^2z^2}{2\pi^3z^3}+O(1)\\[10pt] &\implies\text{residue}=-\frac1{2\pi} \end{align} $$ Result

Thus, $$ \begin{align} \int_0^\infty\frac{\sin(\pi x^2)}{\sinh^2(\pi x)}\,\mathrm{d}x &=-\frac\pi2\left(-\frac1{2\pi}-\frac1{2\pi}+\frac{\sqrt2}{4\pi}+\frac{\sqrt2}{4\pi}\right)\\[6pt] &=\frac{2-\sqrt2}{4} \end{align} $$

  • $\begingroup$ Wow, that's nice robjohn. It sure doesn't bother me that it is two years old. I always enjoy seeing your clever solutions. I can give an upvote. I hate to take Hand greenie away :) $\endgroup$ – Cody Jul 30 '13 at 21:17
  • $\begingroup$ No, I wouldn't want to do that. Now, I may change my mind if I find a real-only solution :-) $\endgroup$ – robjohn Jul 30 '13 at 22:02
  • 3
    $\begingroup$ While it is true that $\sinh(z)\approx z$, one has to be careful since $\frac1{z^3}$ has a residue of $0$ at $z=0$, but $\frac1{\sinh^3(z)}$ has a residue of $-\frac12$ at $z=0$. Rather than looking at $\epsilon e^{it}$, look at the power series near the singularity; e.g. near $0$, $$ \begin{align} \frac{e^{i\pi z^2}e^{\pi z}}{\sinh^2(\pi z)\cosh(\pi z)} &=\frac{\left(1+O\left(z^2\right)\right) \left(1+\pi z+O\left(z^2\right)\right)}{\pi^2z^2\left(1+O\left(z^2\right)\right) \left(1+O\left(z^2\right)\right)}\\ &=\frac1{\pi^2z^2}+\frac1{\pi z}+O(1) \end{align} $$ $\endgroup$ – robjohn Aug 2 '13 at 18:49
  • 1
    $\begingroup$ Thus, that function has a residue of $\frac1\pi$ at $z=0$. If the contour looks similar to mine, but encompasses $0$ and misses $i$, then we only need worry about the sum of the residues at $0$ and $i/2$. However, you are right, and he shouldn't be considering half the residue at points where the singularity is greater than degree 1. $\endgroup$ – robjohn Aug 2 '13 at 18:51
  • 2
    $\begingroup$ $$ \begin{align} &\lim_{\epsilon\to0}\mathrm{Im}\int_\pi^0\left(\color{#C00000} {\frac1{\pi^2\epsilon^2} e^{-2it}}\color{#00A000} {+\frac1{\pi\epsilon} e^{-it}}\right)\epsilon ie^{it}\,\mathrm{d}t\\ =&\lim_{\epsilon\to0}\mathrm{Im}\int_\pi^0\left(\color{#C00000} {\frac1{\pi^2\epsilon} e^{-it}}\color{#00A000} {+\frac1{\pi}}\right)i\,\mathrm{d}t\\ =&\lim_{\epsilon\to0}\mathrm{Im}\left(\color{#C00000}{\frac{-2}{\pi^2\epsilon}} \color{#00A000}{-i}\right)\\ =&\color{#C00000}{0}\color{#00A000}{-i} \end{align} $$ $\endgroup$ – robjohn Aug 2 '13 at 20:11

It can be done using contour integration and the calculus of residues.

Sketch: Integrate $$ f(z) = \frac{e^{i\pi z^2} e^{\pi z}}{\sinh^2 (\pi z) \cosh(\pi z)} $$ around a rectangular contour with corners at $\pm R$ and $\pm R + i$ and with semicircular indentations of radius $\epsilon$ to avoid the poles at $0$ and $i$, take imaginary parts and let $R\to\infty$, $\epsilon\to 0^+$.

You'll need to use $$ f(x)-f(x+i)=\frac{2 e^{i \pi x^2}}{\sinh^2(\pi x)} $$ together with $$ \operatorname*{res}_{z=0} \, f(z) = \operatorname*{res}_{z=i} \, f(z) = \frac{1}{\pi} $$ (since these will each contribute $-i \pi$ times the residue in the limit $\epsilon \to 0^+$) and $$ \operatorname*{res}_{z=i/2} \, f(z) = \frac{-1+i}{\pi\sqrt{2}}. $$

  • $\begingroup$ Thanks everyone for your contributions and help. $\endgroup$ – Cody Sep 4 '11 at 11:42
  • $\begingroup$ I wonder if something involving the differentiation of an integral could be done... $\endgroup$ – J. M. is a poor mathematician Sep 4 '11 at 14:55
  • 2
    $\begingroup$ @J. M.: I half-heartedly tried something like that first, but I didn't get anywhere. I'm not Feynman... ;-) $\endgroup$ – Hans Lundmark Sep 4 '11 at 17:30

I would write the $\sin(x^2)$ as $(e^{ix^2}-e^{-ix^2})/2i$ and the sinh as $(e^{ x}-e^{-x})/2$. Then I'd maybe put the integrand in the form of $(e^{p_1(x)}+e^{p_2(x)}+\cdots)^{-1}+(e^{p_3(x)}+e^{p_4(x)}+\cdots)^{-1}+\cdots$ where $p_i(x)$ are polynomes with complex coefficients. I have no clue if that helps, to be honest.

Another idea would be partial integration after multiplying with 1, like: $\int\mathrm dx 1\cdot f(x)= xf(x)-\int\mathrm dx \; x\cdot f'(x)$ Sometimes this helps to handle a $x^2$ in the argument of a complicated function.


Numeric answer would be possible to get using the following tools:

  1. Riemann integration
  2. root-finding algorithm for the equation
  3. some limit sequence for the infinity giving better and better approximations

Riemann integration is needed to calculate F(x). Basically you'll need a root-finding algorithm that works with G : R->R functions, and gives a single x as solution. Just move the constant to the other side to get F(x)-F(0)-c=0. with G(x)=F(x)-F(0)-c. The infinity will break the riemann integration, so you'll need a sequence like { G(a_1)=0, G(a_2)=0, G(a_3)=0, ... } to get better and better approximations with a_1,a_2,a_3, ... sequence increasing towards infinity. The result then looks like {x_1, x_2,x_3,...} sequence which contains the values of x coming from root-finding algorithm.

But there could be better ways to solve this problem...

EDIT: there is problems with this solution. Namely, the a_i is a constant, not a variable, so root-finding might not be needed after all. All I get is approximation of 0=0.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.