1. Consider the integral $$ \int\limits^\infty_0 \frac{\cos 2 mx \ dx}{\{1+x^2/a^2\} \{1+x^2/(a+1)^2\}\{1+x^2/(a+1)^2\}\cdots}, $$ where $m$ and $a$ are positive. It can be easily proved that \begin{eqnarray*} \left\{1-\left(\frac{t}{a}\right)^2\right\} \left\{1 - \left(\frac{t}{a+1}\right)^2\right\} \left\{1-\left(\frac{t}{a+2}\right)^2\right\} \cdots \left\{1-\left(\frac{t}{a+n-1}\right)^2\right\}\\ = \frac{\Gamma(a+n-t)\Gamma(a+n+t)\{\Gamma (a)\}^2}{\Gamma(a-t) \Gamma(a+t) \{\Gamma(a+n)\}^2 }, \end{eqnarray*} where $n$ is any positive integer. Hence, by splitting $$ \frac{1}{\{1+x^2/a^2\} \{1+x^2/(a+1)^2\}\cdots \{1+x^2/(a+n-1)^2\}}$$ into partial fractions, we see that it is equal to \begin{eqnarray*} \frac{2 \Gamma(2a) \{\Gamma (a+n)\}^2}{\{\Gamma(a)\}^2 \Gamma(n) \Gamma(2a+n)} && \left\{\frac{a}{a^2+x^2} - \frac{2a}{1!} \frac{n-1}{n+2a} \frac{a+1}{(a+1)^2 + x^2}\right. \\ &+&\left. \frac{2a(2a+1)}{2!} \frac{(n-1)(n-2)}{(n+2a)(n+2a+1)} \frac{a+2}{(a+2)^2+x^2} - \cdots \right\}. \end{eqnarray*} Multiplying both sides by $\cos 2 mx$ and integrating from 0 to $\infty$ with respect to $x$, we have \begin{eqnarray*} \int\limits^\infty_0 && \frac{\cos 2 mx \ dx}{\{1+x^2/a^2\} \{1+x^2/(a+1)^2\} \cdots \{1+x^2/(a+n-1)^2\}} \\ &=& \frac{\pi \Gamma (2a) \{\Gamma(a+n)\}^2}{\{\Gamma (a)\}^2 \Gamma (n) \Gamma (2a+n)} \left\{e^{-2a m} - \frac{2a}{1!}\frac{n-1}{n+2a} e^{-2(a+1)m} + \cdots \right\}. \end{eqnarray*} The limit of the right-hand side, as $n \to \infty$, is \begin{eqnarray*} \frac{\pi \Gamma (2a)}{\{\Gamma (a)\}^2} && \left\{e^{-2am} - \frac{2a}{1!} e^{-2(a+1)m} + \frac{2a(2a+1)}{2!} e^{-2(a+2)m} -\cdots \right\}\\ &=& \half\sqrt{\pi} \frac{\Gamma (a+\half )}{\Gamma(a)} \mbox{sech}^{2a} m. \end{eqnarray*} Hence

Since \begin{equation*} \left\{ 1+ \left(\frac{x}{a}\right)^2\right\} \left\{1+ \left(\frac{x}{a+1} \right)^2\right\} \left\{1+\left(\frac{x}{a+2}\right)^2\right\} \cdots = \frac{\{\Gamma(a)\}^2}{\Gamma(a+ix) \Gamma(a-ix)}, \end{equation*} the formula (1) is equivalent to

2. In a similar manner we can prove that \begin{eqnarray*} \int\limits^\infty_0 && \left(\frac{1+x^2/b^2}{1+x^2/a^2}\right) \left(\frac{1+x^2/(b+1)^2}{1+x^2/(a+1)^2}\right) \left(\frac{1+x^2/(b+2)^2}{1+x^2/(a+2)^2}\right) \cdots \cos mx \ dx\\ &=& \frac{\pi \Gamma (2a) \{\Gamma(b)\}^2}{\{\Gamma (a)\}^2 \Gamma(b+a) \Gamma(b-a)} \left\{e^{-am} - \frac{2a}{1!}\frac{b - a- 1}{b+a} e^{-(a+2)m}\right.\\ &+& \left.\frac{2a(2a+1)}{2!} \frac{(b-a-1)(b-a-2)}{(b+a)(b+a+1)} e^{-(a+2)m} - \cdots \right\}, \end{eqnarray*} where $m$ is positive and $0 \lt a \lt b$. When $0 \lt a \lt b-\frac{1}{2}$, the integral and the series remain convergent for $m=0$, and we obtain the formulæ

If $a_1, a_2, a_3, \ldots, a_n$ be $n$ positive numbers in arithmetical progression, then $$\int\limits^\infty_0 \frac{dx}{(a_1^2 + x^2)(a_2^2+x^2)(a_3^2+x^2)\cdots(a_n^2+x^2)}$$ is a particular case of the above integral, and its value can be written down at once. Thus, for example, by putting $a=\frac{11}{10}$ and $b=\frac{61}{10}$, we obtain \begin{eqnarray*} \int\limits^\infty_0 \frac{dx}{(x^2+11^2)(x^2+21^2)(x^2+31^2)(x^2+41^2)(x^2+51^2)}\\ =\frac{5 \pi}{12\cdot 13\cdot 16\cdot 17\cdot 18\cdot 22\cdot 23\cdot 24\cdot 31\cdot 32\cdot 41}. \end{eqnarray*}

3. It follows at once from equation (1), by applying Fourier's theorem $$\int\limits^\infty_0 \cos ny dy \int\limits^\infty_0 \phi (x) \cos xy \ dx = \half\pi \phi (n) ,$$ that, when $a$ and $n$ are positive,

Hence the function $$ \phi(a) =\int\limits^\infty_0\mbox{sech}^a x \cos nx \ dx (0 \lt a \lt 2)$$ satisfies the functional equation $$\phi (a) \phi(2-a) = \frac{\pi \sin \pi a}{2(1-a)(\cosh \pi n -\cos \pi a)}.$$

4. Let $$ \int\limits^b_a f(x) F(nx) \ dx = \Psi (n),$$ and $$ \int\limits^\beta_\alpha \phi(x) F(nx) \ dx =\chi (n).$$ Then, if we suppose the functions $f, \phi,$ and $F$ to be such that the order of integration is indifferent, we have

A number of curious relations between definite integrals may be deduced from this result. We have, for example, the formulæ By applying the general result (6) to the integrals (7) and (8), we obtain $$\sqrt{3} \int\limits^\infty_0 \frac{dx}{\cosh \pi x(1+2 \cosh 2 nx)} =\int\limits^\infty_0 \frac{dx}{\cosh nx(1+2 \cosh \frac{2}{3} \pi x)};$$ or, in other words, if $\alpha \beta = \frac{3}{4} \pi^2$, then In the same way, from (8) and (9), we obtain with the condition $\alpha \beta = \frac{4}{3} \pi;$ and, from (7) and (9), with the conditions $\alpha \beta=\pi$. ¹ Similarly, by taking the two integrals $$ \int\limits^\infty_0 \frac{\sin nx}{e^{2 \pi x}-1} \ dx= \half \left(\frac{1}{e^n-1} + \half - \frac{1}{n}\right),$$ and $$\int\limits^\infty_0 xe^{-x^2} \sin nx \ dx = \frac{1}{4} \sqrt{\pi} ne^{-\frac{1}{4} n^2}, $$ we can prove that, if $\alpha \beta = \pi^2$, then and so on.

5. Suppose now that $a, b$ and $n$ are positive, and

Then, if the conditions of Fourier's double integral theorem are satisfied, we have Applying the formula (6) to (14) and (15), we obtain Thus, when $a=b$, we have the formula $$ \half \pi \int\limits^\infty_0 \phi (x) \phi (nx) \ dx = \int\limits^\infty_0 \Psi (x) \Psi (nx) \ dx, $$ where $$\psi (t) = \int\limits^\infty_0 \phi(x) \begin{array}{cc}\cos \\ \sin \end{array} tx \ dx; $$ and, in particular, if $n=1$, then $$\half \pi \int\limits^\infty_0 \{\phi (x)\}^2 \ dx = \int\limits^\infty_0 \{\Psi (x) \}^2 \ dx.$$ if $$ \phi (a,x) =\frac{1}{\{1+x^2/a^2\}\{1+x^2/(a+1)^2\}\cdots}\quad (a>0),$$ then, by (1), $$ \Psi (a, x) = \half \sqrt{\pi} \frac{\Gamma (a+\half)}{\Gamma(a)} \mbox{sech}^{2a} \half x. $$ Hence, by (16), $$\int\limits^\infty_0 \phi(a,x) \phi(b.x) \ dx = \frac{\Gamma(a+\half )\Gamma(b+\half)}{2\Gamma(a) \Gamma(b)} \int\limits^\infty_0 \mbox{sech}^{2a+2b} \half x \ dx; $$ and so $a$ and $b$ being positive: or As particular cases of the above result, we have, when $b=1$, $$\int\limits^\infty_0 \frac{x}{\sin h \pi x} \frac{dx}{\{1+x^2/a^2\}\{1+x^2/(a+1)^2\}\cdots} = \frac{a}{2(1+2a)} ; $$ when $b=2,$ $$ \int\limits^\infty_0 \frac{x^3}{\sin h \pi x} \frac{dx}{\{1+x^2/a^2\}\{1+x^2/(a+1)^2\}\cdots} = \frac{a^2}{2(1+2a)(3+2a)};$$ and so on. Since $\Pi\{1+x^2/(a+n)^2\}$ can be expressed in finite terms by means of hyperbolic functions when $2a$ is an integer, we can deduce a large number of special formulæ from the preceding results.

6. Another curious formula is the following. If $0 \lt r \lt 1, n \gt 0,$ and $0 \lt a \lt r^{n-1},$ then

unless $n$ is an integer or $a$ is of the form $r^p$, where $p$ is a positive integer. If $a=r^p$, the formula reduces to If $n$ is an integer, the value of the integral is in any case $$ -\frac{\log r}{1-a} \frac{(1-r)(1-r^2)\cdots(1-r^{n-1})}{(r-a)(r^2-a)\cdots(r^{n-1}-a)}. $$ My own proofs of the above results make use of a general formula, the truth of which depends on conditions which I have not yet investigated completely. A direct proof depending on Cauchy's theorem will be found in Mr Hardy's note which follows this paper. The final formula used in Mr Hardy's proof can be proved as follows. Let $$ f(t) = \prod^{m=\infty}_{m=0} \left(\frac{1-btx^m}{1-atx^m}\right) = A_0 + A_1 t+ A_2t^2+ \cdots . $$ Then it is evident that $$ (1-at) f(t) = (1-bt) f(tx).$$ That is $$ (1-at)(A_0 + A_1 t + A_2 t^2 + \cdots) = (1-bt)(A_0 + A_1 tx + A_2 t^2 x^2 + \cdots). $$ Equating the coefficients of $t^n$, we obtain $$ A_n = A_{n-1} \frac{1-bx^{n-1}}{1-x^n}; $$ and $A_0$ is evidently 1. Hence we have

7. As a particular case of (19), we have, when $a=0$,

When $n$ is an integer, the value of the integral reduces to $$ - r^{-\frac{1}{2} n(n-1)} (1-r)(1-r^2)\cdots (1-r^{n-1}) \log r. $$ When we put $n = \half$ in (19), we have If, in particular, $n = \half$ in (22), or $a=0$ in (23), then the $n$th term in the denominator being $r^{\frac{1}{2} n(n-1)}$. Thus, for example, when $r=e^{-5 \pi}$, we have \begin{eqnarray*} \int\limits^\infty_0 && \frac{dx}{(1+x^2)(1+e^{- 10 \pi} x^2) (1+ e^{-20 x^2})\cdots} \nonumber \\ &=& \frac{\pi}{2(1+e^{-5\pi}+e^{-15\pi}+e^{-30\pi}+\cdots )} \nonumber\\ &=& \pi^{\frac{3}{4}} \Gamma \left(\mbox{$\frac{3}{4}$}\right) \sqrt{5}\sqrt[8]{2} \half (1+ \sqrt[4]{5}) \{\half (1+\sqrt{5})\}^{\frac{1}{2}} e^{-5 \pi/8}. \end{eqnarray*} Similarly \begin{eqnarray*} \int\limits^\infty_0 && \frac{dx}{(1+x^2)(1+e^{-20 \pi} x^2)(1+e^{-40 \pi} x^2)\cdots} \nonumber\\ &=& \pi^{\frac{3}{4}} \Gamma \left(\mbox{$\frac{3}{4}$}\right) \sqrt{5} \sqrt[4]{2} \half (1+ \sqrt[4]{5})^2 \{\half (1+\sqrt{5})\} \mbox{$\frac{5}{2}$} e^{-5 \pi/4}; \end{eqnarray*} and \begin{eqnarray*} \int\limits^\infty_0 \frac{dx}{(1+x^2)(1+.001 x^2)(1+.00001 x^2)\cdots} &=& \half \pi \frac{10}{11} \frac{1110}{1111} \frac{111110}{111111}\cdots \\ &=& \frac{\pi}{2.202~002~000~200~002~000~002\ldots}. \end{eqnarray*}