Ramanujan's master theorem

Template:Short description In mathematics, Ramanujan's master theorem, named after Srinivasa Ramanujan,Template:Sfn is a technique that provides an analytic expression for the Mellin transform of an analytic function.

Page from Ramanujan's notebook stating his Master theorem.

The result is stated as follows:

If a complex-valued function $f(x)$ has an expansion of the form $$f(x)=\sum _{k=0}^{\mathrm{\infty }}\frac{\phi (k)}{k!}(-x{)}^k$$

where $\phi (s)$ is an analytic function, then the Mellin transform of $f(x)$ is given by

$${\int }_0^{\mathrm{\infty }}{x}^{s-1}f(x)dx=\mathrm{\Gamma }(s)\phi (-s)$$

where $\mathrm{\Gamma }(s)$ is the gamma function.

It was widely used by Ramanujan to calculate definite integrals and infinite series.

Higher-dimensional versions of this theorem also appear in quantum physics through Feynman diagrams.Template:Sfn

A similar result was also obtained by Glaisher.Template:Sfn

Alternative formalism

An alternative formulation of Ramanujan's master theorem is as follows:

$${\int }_0^{\mathrm{\infty }}{x}^{s-1}(\lambda (0)-x\lambda (1)+x^2\lambda (2)-\cdots )dx=\frac{\pi }{\sin (\pi s)}\lambda (-s)$$

which gets converted to the above form after substituting $\lambda (n)\equiv \frac{\phi (n)}{\mathrm{\Gamma }(1+n)}$ and using the functional equation for the gamma function.

The integral above is convergent for $0<ℛℯ(s)<1$ subject to growth conditions on $\phi$.Template:Sfn

Proof

A proof subject to "natural" assumptions (though not the weakest necessary conditions) to Ramanujan's master theorem was provided by G. H. HardyTemplate:Sfn(chapter XI) employing the residue theorem and the well-known Mellin inversion theorem. Template:Collapse top

Recall Euler's representation of the Gamma function ${\int }_0^{\mathrm{\infty }}{x}^{n-1}{e}^{-mx}dx=m^{-n}\mathrm{\Gamma }(n)$

choosing $m=r^k$, multiplying both sides by $\frac{f^k(a)h^k}{k!}$ and then summing over $0\le k$ to obtain: $\sum _{k=0}^{\mathrm{\infty }}\frac{f^k(a)h^k}{k!}{\int }_0^{\mathrm{\infty }}{x}^{n-1}{e}^{-r^{k}x}dx=\mathrm{\Gamma }(n)\sum _{k=0}^{\mathrm{\infty }}\frac{f^k(a)(h{r}^{-n}{)}^k}{k!}$

Observing the sum in the RHS is a Taylor Series and writing $e^{-r^{k}x}$ in its series expansion:

$\sum _{k=0}^{\mathrm{\infty }}\frac{f^k(a)h^k}{k!}{\int }_0^{\mathrm{\infty }}{x}^{n-1}\sum _{j=0}^{\mathrm{\infty }}\frac{(-r^{k}x{)}^j}{j!}dx=\mathrm{\Gamma }(n)f(h{r}^{-n}+a)$

Rewriting the LHS:

$\sum _{j=0}^{\mathrm{\infty }}\frac{(-1{)}^j}{j!}{\int }_0^{\mathrm{\infty }}{x}^{n-1+j}\sum _{k=0}^{\mathrm{\infty }}\frac{f^k(a)(h{r}^j{)}^k}{k!}dx=\mathrm{\Gamma }(n)f(h{r}^{-n}+a)$

Then once again observing the sum over $k$ is a Taylor series: ${\int }_0^{\mathrm{\infty }}{x}^{n-1}\sum _{j=0}^{\mathrm{\infty }}\frac{f^k(h{r}^j+a)(-x{)}^j}{j!}dx=\mathrm{\Gamma }(n)f(h{r}^{-n}+a)$

Finally defining $f(h{r}^N+a)=\varphi (N)$ and letting $F(x)=\sum _{j=0}^{\mathrm{\infty }}\frac{\varphi (j)(-x{)}^j}{j!}$ we gain the master theorem: ${\int }_0^{\mathrm{\infty }}{x}^{n-1}F(x)dx=\mathrm{\Gamma }(n)\varphi (-n)$ Template:Collapse bottom

Template:Collapse top Let $\lambda (z)$ be an analytic single-valued function defined on a half plane

$H(\delta )=\{z\in ℂ:\Re (z)\ge -\delta \}$

for some $0<\delta <1$. Suppose that, for some $A<\pi$, $\lambda$ satisfies the growth condition

$|\lambda (v+iw)|<C{e}^{Pv+A|w|}$

for all $z=v+iw\in H(\delta )$. Let $0<x<e^{-P}$. This implies the series $F(x)=\sum _{n=0}^{\mathrm{\infty }}(-x{)}^n\lambda (n)$ converges. Observing $\frac{\pi }{\sin (\pi s)}$ has poles at $s\in {\mathbb{Z}}^{-}$ with residue $(-1{)}^{-s}$, application of the residue theorem yields

$F(x)=\frac{1}{2\pi i}{\int }_{c-i\mathrm{\infty }}^{c+i\mathrm{\infty }}\frac{\pi }{\sin (\pi s)}\lambda (-s)x^{-s}ds$

for any $0<c<\delta$. This integral converges absolutely and uniformly for $c\in (a,b)$ for any $0<a,b<\delta$. Applying the Mellin inversion theorem yields

${\int }_0^{\mathrm{\infty }}{x}^{s-1}F(x)dx=\frac{\pi }{\sin (\pi s)}\lambda (-s)$ Template:Collapse bottom

Application to Bernoulli polynomials

The generating function of the Bernoulli polynomials $B_k(x)$ is given by:

$$\frac{z{e}^{xz}}{e^z-1}=\sum _{k=0}^{\mathrm{\infty }}{B}_k(x)\frac{z^k}{k!}$$

These polynomials are given in terms of the Hurwitz zeta function:

$$\zeta (s,a)=\sum _{n=0}^{\mathrm{\infty }}\frac{1}{(n+a{)}^s}$$

by $\zeta (1-n,a)=-\frac{B_n(a)}{n}$ for $n\ge 1$. Using the Ramanujan master theorem and the generating function of Bernoulli polynomials one has the following integral representation:Template:Sfn

$${\int }_0^{\mathrm{\infty }}{x}^{s-1}(\frac{e^{-ax}}{1-e^{-x}}-\frac{1}{x})dx=\mathrm{\Gamma }(s)\zeta (s,a)$$

which is valid for $0<ℛℯ(s)<1$.

Application to the gamma function

Weierstrass's definition of the gamma function

$$\mathrm{\Gamma }(x)=\frac{e^{-\gamma x}}{x}\prod _{n=1}^{\mathrm{\infty }}{(1+\frac{x}{n})}^{-1}{e}^{x/n}$$

is equivalent to expression

$$\log \mathrm{\Gamma }(1+x)=-\gamma x+\sum _{k=2}^{\mathrm{\infty }}\frac{\zeta (k)}{k}(-x{)}^k$$

where $\zeta (k)$ is the Riemann zeta function.

Then applying Ramanujan master theorem we have:

$${\int }_0^{\mathrm{\infty }}{x}^{s-1}\frac{\gamma x+\log \mathrm{\Gamma }(1+x)}{x^2}\mathrm{d}x=\frac{\pi }{\sin (\pi s)}\frac{\zeta (2-s)}{2-s}$$

valid for $0<ℛℯ(s)<1$.

Special cases of $s=\frac{1}{2}$ and $s=\frac{3}{4}$ are

$${\int }_0^{\mathrm{\infty }}\frac{\gamma x+\log \mathrm{\Gamma }(1+x)}{x^{5/2}}\mathrm{d}x=\frac{2\pi }{3}\zeta \left(\frac{3}{2}\right)$$

$${\int }_0^{\mathrm{\infty }}\frac{\gamma x+\log \mathrm{\Gamma }(1+x)}{x^{9/4}}\mathrm{d}x=\sqrt{2}\frac{4\pi }{5}\zeta \left(\frac{5}{4}\right)$$

Application to Bessel functions

The Bessel function of the first kind has the power series $$J_{\nu }(z)=\sum _{k=0}^{\mathrm{\infty }}\frac{(-1{)}^k}{\mathrm{\Gamma }(k+\nu +1)k!}(\frac{z}{2}{)}^{2k+\nu }$$

By Ramanujan's master theorem, together with some identities for the gamma function and rearranging, we can evaluate the integral

$$\frac{2^{\nu -2s}\pi }{\sin (\pi (s-\nu ))}{\int }_0^{\mathrm{\infty }}{z}^{s-1-\nu /2}{J}_{\nu }(\sqrt{z})dz=\mathrm{\Gamma }(s)\mathrm{\Gamma }(s-\nu )$$

valid for $0<2ℛℯ(s)<ℛℯ(\nu )+\frac{3}{2}$.

Equivalently, if the spherical Bessel function $j_{\nu }(z)$ is preferred, the formula becomes

$$\frac{2^{\nu -2s}\sqrt{\pi }(1-2s+2\nu )}{\cos (\pi (s-\nu ))}{\int }_0^{\mathrm{\infty }}{z}^{s-1-\nu /2}{j}_{\nu }(\sqrt{z})dz=\mathrm{\Gamma }(s)\mathrm{\Gamma }(\frac{1}{2}+s-\nu )$$

valid for $0<2ℛℯ(s)<ℛℯ(\nu )+2$.

The solution is remarkable in that it is able to interpolate across the major identities for the gamma function. In particular, the choice of $J_{1/2}(\sqrt{z})$ gives the square of the gamma function, $j_0(\sqrt{z})$ gives the duplication formula, $z^{-1/2}{J}_1(\sqrt{z})$ gives the reflection formula, and fixing to the evaluable $s=\frac{1}{2}$ or $s=1$ gives the gamma function by itself, up to reflection and scaling.

Bracket integration method

The bracket integration method (method of brackets) applies Ramanujan's master theorem to a broad range of integrals.Template:Sfn The bracket integration method generates the integrand's series expansion, creates a bracket series, identifies the series coefficient and formula parameters and computes the integral.Template:Sfn

Integration formulas

This section identifies the integration formulas for integrands with and without consecutive integer exponents and for single and double integrals. The integration formula for double integrals may be generalized to any multiple integral. In all cases, there is a parameter value $n^{\ast }$ or array of parameter values $N^{\ast }$ that solves one or more linear equations derived from the exponent terms of the integrand's series expansion.

Consecutive integer exponents, 1 variable

This is the function series expansion, integral and integration formula for an integral whose integrand's series expansion contains consecutive integer exponents.Template:Sfn $$\begin{array}{cc}& f(y)=\sum _{n=0}^{\mathrm{\infty }}\frac{(-1{)}^n}{n!}\phi (n)y^n\\ & {\int }_0^{\mathrm{\infty }}{y}^{c-1}f(y)dy\\ & ={\int }_0^{\mathrm{\infty }}\sum _{n=0}^{\mathrm{\infty }}\frac{(-1{)}^n}{n!}\phi (n)y^{n+c-1}dy\\ & =\mathrm{\Gamma }(-n^{\ast })\phi (n^{\ast }).\end{array}$$ The parameter $n^{\ast }$ is a solution to this linear equation. $$n^{\ast }+c=0,n^{\ast }=-c$$

General exponents, 1 variable

Applying the substitution $y=x^a$ generates the function series expansion, integral and integration formula for an integral whose integrand's series expansion may not contain consecutive integer exponents.Template:Sfn $$\begin{array}{cc}& f(x)=\sum _{n=0}^{\mathrm{\infty }}\frac{(-1{)}^n}{n!}\phi (n)x^{an}\\ & {\int }_0^{\mathrm{\infty }}{x}^{c-1}f(x)dx\\ & ={\int }_0^{\mathrm{\infty }}\sum _{n=0}^{\mathrm{\infty }}\frac{(-1{)}^n}{n!}\phi (n)x^{an+c-1}dx\\ & =a^{-1}\mathrm{\Gamma }(-n^{\ast })\phi (n^{\ast }).\end{array}$$ The parameter $n^{\ast }$ is a solution to this linear equation. $$a{n}^{\ast }+c=0,n^{\ast }=-a^{-1}c$$

Consecutive integer exponents, double integral

This is the function series expansion, integral and integration formula for a double integral whose integrand's series expansion contains consecutive integer exponents.Template:Sfn $$\begin{array}{cc}& f(y_1,y_2)=\sum _{n=0}^{\mathrm{\infty }}\frac{(-1{)}^{n_1}}{n_1!}\frac{(-1{)}^{n_2}}{n_2!}\phi (n_1,n_2)y_1^{n_1}{y}_2^{n_2}\\ & {\int }_0^{\mathrm{\infty }}{y}_1^{c_1-1}{y}_2^{c_2-1}f(y_1,y_2)d{y}_{1}d{y}_2\\ & ={\int }_0^{\mathrm{\infty }}{\int }_0^{\mathrm{\infty }}\sum _{n_1=0}^{\mathrm{\infty }}\sum _{n_2=0}^{\mathrm{\infty }}\frac{(-1{)}^{n_1}}{n_1!}\frac{(-1{)}^{n_2}}{n_2!}\phi (n_1,n_2)y_1^{n_1+c_1-1}{y}_2^{n_2+c_2-1}d{y}_{1}d{y}_2\\ & =\mathrm{\Gamma }(-n_1^{\ast })\mathrm{\Gamma }(-n_2^{\ast })\phi (n_1^{\ast },n_2^{\ast }).\end{array}$$ The parameters $n_1^{\ast }$ and $n_2^{\ast }$ are solutions to these linear equations. $$n_1^{\ast }+c_1=0,n_2^{\ast }+c_2=0,n_1^{\ast }=-c_1,n_2^{\ast }=-c_2$$

General exponents, double integral

This section describes the integration formula for a double integral whose integrand's series expansion may not contain consecutive integer exponents. Matrices contain the parameters needed to express the exponents in a series expansion of the integrand, and the determinant of invertible matrix $A$ is $det|A|$.Template:Sfn $$A=\left|\begin{array}{cc}{a}_{11}& a_{12}\\ a_{21}& a_{22}\end{array}\right|,C=\left|\begin{array}{c}{c}_1\\ c_2\end{array}\right|,N^{\ast }=\left|\begin{array}{c}{n}_1^{\ast }\\ n_2^{\ast }\end{array}\right|$$ Applying the substitution $$y_1=x_1^{a_{11}}{x}_2^{a_{21}},y_2=x_1^{a_{12}}{x}_2^{a_{22}}$$ generates the function series expansion, integral and integration formula for a double integral whose integrand's series expansion may not contain consecutive integer exponents.Template:Sfn The integral and integration formula areTemplate:SfnTemplate:Sfn $$\begin{array}{cc}& {\int }_0^{\mathrm{\infty }}{\int }_0^{\mathrm{\infty }}\sum _{n_1=0}^{\mathrm{\infty }}\sum _{n_2=0}^{\mathrm{\infty }}\frac{(-1{)}^{n_1}}{n_1!}\frac{(-1{)}^{n_2}}{n_2!}\phi (n_1,n_2)x_1^{n_1{a}_{11}+n_2{a}_{12}+c_1-1}{x}_2^{n_1{a}_{21}+n_2{a}_{22}+c_2-1}d{x}_{1}d{x}_2\\ & =det|A{|}^{-1}\mathrm{\Gamma }(-n_1^{\ast })\mathrm{\Gamma }(-n_2^{\ast })\phi (n_1^{\ast },n_2^{\ast }).\end{array}$$ The parameter matrix $N^{\ast }$ is a solution to this linear equation.Template:Sfn $$A{N}^{\ast }+C=0,N^{\ast }=-A^{-1}C$$.

Positive complexity index

In some cases, there may be more sums than variables. For example, if the integrand is a product of 3 functions of a common single variable, and each function is converted to a series expansion sum, the integrand is now a product of 3 sums, each sum corresponding to a distinct series expansion.

• The number of brackets is the number of linear equations associated with an integral. This term reflects the common practice of bracketing each linear equation.Template:Sfn
• The complexity index is the number of integrand sums minus the number of brackets (linear equations). Each series expansion of the integrand contributes one sum.Template:Sfn
• The summation indices (variables) are the indices that index terms in a series expansion. In the example, there are 3 summation indices $n_1,n_2$ and $n_3$ because the integrand is a product of 3 series expansions.Template:Sfn
• The free summation indices (variables) are the summation indices that remain after completing all integrations. Integration reduces the number of sums in the integrand by replacing the series expansions (sums) with an integration formula. Therefore, there are fewer summation indices after integration. The number of chosen free summation indices equals the complexity index.Template:Sfn

Integrals with a positive complexity index

The free summation indices ${\overline{n}}_1,\dots ,{\overline{n}}_f$ are elements of set $F$. The matrix of free summation indices is $\overline{N}$ and the coefficients of the free summation indices is matrix $\overline{A}$. $$\overline{A}=\left|\begin{array}{ccc}{\overline{a}}_{11}& \dots & {\overline{a}}_{1f}\\ ⋮& & ⋮\\ {\overline{a}}_{b1}& \dots & {\overline{a}}_{bf}\end{array}\right|,\overline{N}=\left|\begin{array}{c}{\overline{n}}_1\\ ⋮\\ {\overline{n}}_f\end{array}\right|$$ The remaining indices are set $B$ containing indices $n_1,\dots ,n_b$. Matrices $A,C$ and $N^{\ast }$ contain matrix elements that multiply or sum with the non-summation indices. The selected free summation indices must leave matrix $A$ non-singular. $$A=\left|\begin{array}{ccc}{a}_{11}& \dots & a_{1b}\\ ⋮& & ⋮\\ a_{b1}& \dots & a_{bb}\end{array}\right|,C=\left|\begin{array}{c}{c}_1\\ ⋮\\ c_b\end{array}\right|,N^{\ast }=\left|\begin{array}{c}{n}_1^{\ast }\\ ⋮\\ n_b^{\ast }\end{array}\right|$$. This is the function's series expansion, integral and integration formula.Template:Sfn $$\begin{array}{cc}& f(x_1,\dots ,x_b)\\ & =\sum _{n\in B}^{\mathrm{\infty }}\sum _{\overline{n}\in F}^{\mathrm{\infty }}\frac{(-1{)}^{n_1}}{n_1!}\frac{(-1{)}^{{\overline{n}}_1}}{{\overline{n}}_1!}\dots \frac{(-1{)}^{n_b}}{n_b!}\frac{(-1{)}^{\overline{n_f}}}{{\overline{n}}_f!}\phi (n_1,\dots ,n_b,{\overline{n}}_1,\dots ,{\overline{n}}_f)\prod _{x_k\in B}{x}_k^{n_1{a}_{k1}+\dots +{\overline{n}}_1{\overline{a}}_{k1}+\dots +c_k-1}\\ & {\int }_0^{\mathrm{\infty }}\dots {\int }_0^{\mathrm{\infty }}{x}^{c_1-1}\dots x^{c_b-1}f(x_1,\dots ,x_b)d{x}_1\dots d{x}_b\\ & =det|A{|}^{-1}\sum _{\overline{n}\in F}^{\mathrm{\infty }}\frac{(-1{)}^{{\overline{n}}_1}}{{\overline{n}}_1!}\dots \frac{(-1{)}^{{\overline{n}}_f}}{{\overline{n}}_f!}\mathrm{\Gamma }(-n_1^{\ast })\dots \mathrm{\Gamma }(-n_b^{\ast })\phi (n_1^{\ast },\dots ,n_b^{\ast },{\overline{n}}_1,\dots ,{\overline{n}}_f).\end{array}$$ The parameters $n_1^{\ast },\dots ,n_b^{\ast }$ are linear functions of the parameters ${\overline{n}}_1^{\ast },\dots ,{\overline{n}}_f^{\ast }$.Template:Sfn $$A{N}^{\ast }+\overline{A}\overline{N}+C=0,N^{\ast }=-A^{-1}(\overline{A}\overline{N}+C)$$.

Bracket series

Table 1. Bracket series notations

Notation type	Power series notation	Bracket series notation
Indicator	$$\frac{(-1{)}^n}{n!}$$	$${\varphi }_n$$
Multi-indicator	$$\prod _{j=1}^N\left(\frac{(-1{)}^{n_j}}{n_j!}\right)$$	$${\varphi }_{n_1,\dots ,n_N}$$
Bracket	$${\int }_0^{\mathrm{\infty }}dx{x}^{a_1{n}_1+\dots +a_m{n}_m+c-1}$$	$$⟨a_1{n}_1+\dots +a_m{n}_m+c⟩$$

Bracket series notations are notations that substitute for common power series notations (Table 1).Template:Sfn Replacing power series notations with bracket series notations transforms the power series to a bracket series. A bracket series facilitates identifying the formula parameters needed for integration. It is also recommended to replace a sum raised to a power:Template:Sfn $$\frac{1}{(x_1+\dots +x_b{)}^{\alpha }}$$ with this bracket series expression:$$\sum _{m_1=0}^{\mathrm{\infty }}\dots \sum _{m_b=0}^{\mathrm{\infty }}{\varphi }_{m_1,\dots ,m_b}{x}_1^{m_1}\dots x_b^{m_b}\frac{⟨\alpha +m_1+\dots +m_b⟩}{\mathrm{\Gamma }(\alpha )}.$$

Algorithm

This algorithm describes how to apply the integral formulas.Template:SfnTemplate:SfnTemplate:Sfn

Table 2. Integral formulas

Complexity index	Integral formula
Zero, single integral	$$a^{-1}\mathrm{\Gamma }(-n^{\ast })\phi (n^{\ast })$$
Zero, multiple integral	$$det|A{|}^{-1}\mathrm{\Gamma }(-n_1^{\ast })\dots \mathrm{\Gamma }(-n_b^{\ast })\phi (n_1^{\ast },\dots ,n_b^{\ast })$$
Positive	$$det|A{|}^{-1}\sum _{\overline{n}\in F}^{\mathrm{\infty }}\frac{(-1{)}^{{\overline{n}}_1}}{{\overline{n}}_1!}\dots \frac{(-1{)}^{{\overline{n}}_f}}{{\overline{n}}_f!}\mathrm{\Gamma }(-n_1^{\ast })\dots \mathrm{\Gamma }(-n_b^{\ast })\phi (n_1^{\ast },\dots ,n_b^{\ast },{\overline{n}}_1,\dots ,{\overline{n}}_f)$$

Input Integral expression

Output Integral value or integral cannot be assigned a value

1. Express the integrand as a power series.
2. Transform the integrand's power series to a bracket series.
3. Obtain the complexity index, formula parameters and series coefficient function.
   1. Complexity index is the number of integrand sums minus number of brackets.
   2. Parameters $n^{\ast }$ or array $N^{\ast }$ are solutions to linear equations $a{n}^{\ast }+c=0$ (zero complexity index, single integral), $A{N}^{\ast }+C=0$ (zero complexity index, single integral) or $A{N}^{\ast }+\overline{A}\overline{N}+C=0$ (positive complexity index).
   3. Identify parameter $a$ or (zero complexity index, single integral) or compute $det|A|$ (all other cases) from the associated linear equations.
   4. Identify the series coefficient function $\phi ()$ of the bracket series.
4. If the complexity index is negative, return integral cannot be assigned a value.
5. If the complexity index is zero, select the formula from table 2 for zero complexity index, single or multiple integral, compute the integral value with this formula, and return this integral value.
6. If the complexity index is positive, select the formula from table 2 for positive complexity index, and compute the integral value as a series expansion with this formula for all possible choices of the free summation indices. Select the lowest complexity index, convergent series expansion, adding series that converge in the same region.
   1. If all series expansions are divergent series or null series (all series terms zero), then return integral cannot be assigned a value.
   2. If the series expansion is non-null and non-divergent, return this series expansion as the integral value.

Examples

Zero complexity index

The bracket method will integrate this integral. $${\int }_0^{\mathrm{\infty }}{x}^{3/2}{e}^{-x^3/2}dx$$

1. Express the integrand as a power series.$${\int }_0^{\mathrm{\infty }}\sum _{n=0}^{\mathrm{\infty }}{2}^{-n}\frac{(-1{)}^n}{n!}{x}^{(3\cdot n+5/2)-1}dx$$
2. Transform the power series to a bracket series. $$\sum _{n=0}^{\mathrm{\infty }}{2}^{-n}\varphi (n)\cdot ⟨3n+\frac{5}{2}⟩$$
3. Obtain the complexity index, formula parameters and series coefficient function.

Complexity index is zero.

$3n^{\ast }+5/2=0$

$n^{\ast }=-5/6,a=3$

$\phi (n)=2^{-n}$.

4. Use table 2 to compute the integral.

$${\int }_0^{\mathrm{\infty }}{x}^{3/2}\cdot e^{-x^3/2}dx$$ $$=a^{-1}\mathrm{\Gamma }(-n^{\ast })\phi (n^{\ast })$$ $$=\frac{\mathrm{\Gamma }\left(\frac{5}{6}\right)2^{5/6}}{3}$$

Positive complexity index

The bracket method will integrate this integral. $${\int }_0^{\mathrm{\infty }}\frac{1}{(1+x^3+x^5{)}^{1/2}}dx$$ 1. Express the integrand as a power series. Use the sum raised to a power formula. $${\int }_0^{\mathrm{\infty }}\sum _{n_1,n_2,n_3}\frac{1}{\sqrt{\mathrm{\Gamma }(1/2)}}{\varphi }_{123}{1}^{n_1}{x}^{5n_2+3n_3}⟨n_1+n_2+n_3+1/2⟩dx$$ 2. Transform the power series to a bracket series. $${\int }_0^{\mathrm{\infty }}\sum _{n_1,n_n,n_3}\frac{1}{\sqrt{\mathrm{\Gamma }(1/2)}}{\varphi }_{123}⟨5n_2+3n_3+1⟩⟨n_1+n_2+n_3+1/2⟩$$ 3. Obtain the complexity index, formula parameters and series coefficient function.

Complexity index is 1 as 3 sums and 2 brackets.

Select $n_3$ as the free index, ${\overline{n}}_3$. The linear equations, solutions, determinant and series coefficient are

$$5n_2^{\ast }+3{\overline{n}}_3+1=0,n_1^{\ast }+n_2^{\ast }+{\overline{n}}_3+1/2=0$$ $$\left|\begin{array}{cc}1& 1\\ 0& 5\end{array}\right|\left|\begin{array}{c}{n}_1^{\ast }\\ n_2^{\ast }\end{array}\right|+\left|\begin{array}{c}1\\ 3\end{array}\right|\left|\begin{array}{c}{\overline{n}}_3\end{array}\right|+\left|\begin{array}{c}1/2\\ 1\end{array}\right|=0$$ $$A{N}^{\ast }+\overline{A}\overline{N}+C=0$$ $$det|A|=5$$ $$n_1^{\ast }=-\frac{2}{5}{\overline{n}}_3-\frac{3}{10},n_2^{\ast }=-\frac{3}{5}{\overline{n}}_3-\frac{1}{5}.$$ $$\phi (n_1^{\ast },n_2^{\ast },{\overline{n}}_3)=\frac{1}{\sqrt{\mathrm{\Gamma }(1/2)}}=\frac{1}{\sqrt{\pi }}$$ 4. Use table 2 to compute the integral $$\begin{array}{cc}& {\int }_0^{\mathrm{\infty }}\frac{1}{(1+x^3+x^5{)}^{1/2}}dx\\ & =\sum _{{\overline{n}}_3=0}^{\mathrm{\infty }}\frac{(-1{)}^{{\overline{n}}_3}}{{\overline{n}}_3!}det|A{|}^{-1}\mathrm{\Gamma }(-n_1^{\ast })\mathrm{\Gamma }(-n_2^{\ast })\phi (n_1^{\ast },n_2^{\ast },{\overline{n}}_3)\\ & =\sum _{{\overline{n}}_3=0}^{\mathrm{\infty }}\frac{(-1{)}^{{\overline{n}}_3}}{{\overline{n}}_3!}\frac{\mathrm{\Gamma }(\frac{2}{5}{\overline{n}}_3+\frac{3}{10})\mathrm{\Gamma }(\frac{3}{5}{\overline{n}}_3+\frac{1}{5})}{5\sqrt{\pi }}\end{array}$$

References

• Template:Cite journal
• Template:Cite journal
• Template:Cite book
• Template:Cite journal
• Template:Cite journal
• Template:Cite arXiv
• Template:Cite journal
• Template:Cite journal
• Template:Cite book
• Template:Cite journal
• Template:Cite journal
• Template:Cite book

External links

• Template:MathWorld