Digamma function

The color representation of the digamma function,

\psi (z)

, in a rectangular region of the complex plane

In mathematics, the digamma function is defined as the logarithmic derivative of the gamma function:^[1]^[2]

\psi (x)={\frac {d}{dx}}\ln {\big (}\Gamma (x){\big )}={\frac {\Gamma '(x)}{\Gamma (x)}}.

It is the first of the polygamma functions.

The digamma function is often denoted as $ψ 0 (x)$ , $ψ (0) (x)$ or $Ϝ$ (the uppercase form of the archaic Greek consonant digamma meaning double-gamma).

Relation to harmonic numbers

The gamma function obeys the equation

\Gamma (z+1)=z\Gamma (z).\,

Taking the derivative with respect to $z$ gives:

\Gamma '(z+1)=z\Gamma '(z)+\Gamma (z)\,

Dividing by $Γ(z + 1)$ or the equivalent $z Γ(z)$ gives:

{\frac {\Gamma '(z+1)}{\Gamma (z+1)}}={\frac {\Gamma '(z)}{\Gamma (z)}}+{\frac {1}{z}}

or:

\psi (z+1)=\psi (z)+{\frac {1}{z}}

Since the harmonic numbers are defined as

H_{n}=\sum _{k=1}^{n}{\frac {1}{k}}

the digamma function is related to it by:

\psi (n)=H_{n-1}-\gamma

where $H n$ is the $n$ th harmonic number, and $γ$ is the Euler–Mascheroni constant. For half-integer values, it may be expressed as

\psi \left(n+{\tfrac {1}{2}}\right)=-\gamma -2\ln 2+\sum _{k=1}^{n}{\frac {2}{2k-1}}

Harmonic Mean Value Inequality

Inspired by the harmonic mean value inequality for the classical gamma function, Horzt Alzer and Graham Jameson proved, among other things, a harmonic mean-value inequality for the digamma function:

$-\gamma \leq {\frac {2\psi (x)\psi ({\frac {1}{x}})}{\psi (x)+\psi ({\frac {1}{x}})}}$ for $x>0$

Equality holds if and only if $x=1$ ^[3]

Integral representations

If the real part of $x$ is positive then the digamma function has the following integral representation

\psi (x)=\int _{0}^{\infty }\left({\frac {e^{-t}}{t}}-{\frac {e^{-xt}}{1-e^{-t}}}\right)\,dt.

This may be written as

\psi (s+1)=-\gamma +\int _{0}^{1}\left({\frac {1-x^{s}}{1-x}}\right)\,dx

which follows from Leonhard Euler's integral formula for the harmonic numbers.

Infinite product representation

The function $\psi (z)/\Gamma (z)$ is an entire function,^[4] and it can be represented by the infinite product

{\frac {\psi (z)}{\Gamma (z)}}=-e^{2\gamma z}\prod _{k=0}^{\infty }\left(1-{\frac {z}{x_{k}}}\right)e^{\frac {z}{x_{k}}}.

Here $x_{k}$ is the kth zero of $\psi$ (see below), and $\gamma$ is the Euler–Mascheroni constant.

Series formula

The digamma function can be computed in the complex plane outside negative integers (Abramowitz and Stegun 6.3.16),^[1] using

\psi (z+1)=-\gamma +\sum _{n=1}^{\infty }{\frac {z}{n(n+z)}}\qquad z\neq -1,-2,-3,\ldots

or

\psi (z)=-\gamma +\sum _{n=0}^{\infty }{\frac {z-1}{(n+1)(n+z)}}=-\gamma +\sum _{n=0}^{\infty }\left({\frac {1}{n+1}}-{\frac {1}{n+z}}\right)\qquad z\neq 0,-1,-2,-3,\ldots

This can be utilized to evaluate infinite sums of rational functions, i.e.,

\sum _{n=0}^{\infty }u_{n}=\sum _{n=0}^{\infty }{\frac {p(n)}{q(n)}},

where $p (n)$ and $q (n)$ are polynomials of $n$ .

Performing partial fraction on $u n$ in the complex field, in the case when all roots of $q (n)$ are simple roots,

u_{n}={\frac {p(n)}{q(n)}}=\sum _{k=1}^{m}{\frac {a_{k}}{n+b_{k}}}.

For the series to converge,

\lim _{n\to \infty }nu_{n}=0,

otherwise the series will be greater than the harmonic series and thus diverge. Hence

\sum _{k=1}^{m}a_{k}=0,

and

{\begin{aligned}\sum _{n=0}^{\infty }u_{n}&=\sum _{n=0}^{\infty }\sum _{k=1}^{m}{\frac {a_{k}}{n+b_{k}}}\\&=\sum _{n=0}^{\infty }\sum _{k=1}^{m}a_{k}\left({\frac {1}{n+b_{k}}}-{\frac {1}{n+1}}\right)\\&=\sum _{k=1}^{m}\left(a_{k}\sum _{n=0}^{\infty }\left({\frac {1}{n+b_{k}}}-{\frac {1}{n+1}}\right)\right)\\&=-\sum _{k=1}^{m}a_{k}{\big (}\psi (b_{k})+\gamma {\big )}\\&=-\sum _{k=1}^{m}a_{k}\psi (b_{k}).\end{aligned}}

With the series expansion of higher rank polygamma function a generalized formula can be given as

\sum _{n=0}^{\infty }u_{n}=\sum _{n=0}^{\infty }\sum _{k=1}^{m}{\frac {a_{k}}{(n+b_{k})^{r_{k}}}}=\sum _{k=1}^{m}{\frac {(-1)^{r_{k}}}{(r_{k}-1)!}}a_{k}\psi ^{r_{k}-1}(b_{k}),

provided the series on the left converges.

Taylor series

The digamma has a rational zeta series, given by the Taylor series at $z = 1$ . This is

\psi (z+1)=-\gamma -\sum _{k=1}^{\infty }\zeta (k+1)(-z)^{k},

which converges for $| z | < 1$ . Here, $ζ (n)$ is the Riemann zeta function. This series is easily derived from the corresponding Taylor's series for the Hurwitz zeta function.

Newton series

The Newton series for the digamma, sometimes referred to as Stern series,^[5]^[6] reads

\psi (s+1)=-\gamma -\sum _{k=1}^{\infty }{\frac {(-1)^{k}}{k}}{\binom {s}{k}}

where $(s k)$ is the binomial coefficient. It may also be generalized to

\Psi (s+1)=-\gamma -{\frac {1}{m}}\sum _{k=1}^{m-1}{\frac {m-k}{s+k}}-{\frac {1}{m}}\sum _{k=1}^{\infty }{\frac {(-1)^{k}}{k}}\left\{{\binom {s+m}{k+1}}-{\binom {s}{k+1}}\right\},\qquad \Re (s)>-1,

where $m = 2,3,4,...$ ^[6]

Series with Gregory's coefficients, Cauchy numbers and Bernoulli polynomials of the second kind

There exist various series for the digamma containing rational coefficients only for the rational arguments. In particular, the series with Gregory's coefficients $G n$ is

\Psi (v)=\ln v-\sum _{n=1}^{\infty }{\frac {{\big |}G_{n}{\big |}(n-1)!}{(v)_{n}}},\qquad \Re (v)>0,

\Psi (v)=2\ln \Gamma (v)-2v\ln v+2v+2\ln v-\ln 2\pi -2\sum _{n=1}^{\infty }{\frac {{\big |}G_{n}(2){\big |}}{(v)_{n}}}\,(n-1)!,\qquad \Re (v)>0,

\Psi (v)=3\ln \Gamma (v)-6\zeta '(-1,v)+3v^{2}\ln {v}-{\frac {3}{2}}v^{2}-6v\ln(v)+3v+3\ln {v}-{\frac {3}{2}}\ln 2\pi +{\frac {1}{2}}-3\sum _{n=1}^{\infty }{\frac {{\big |}G_{n}(3){\big |}}{(v)_{n}}}\,(n-1)!,\qquad \Re (v)>0,

where $(v) n$ is the rising factorial $(v) n = v (v +1)(v +2) ... (v + n -1)$ , $G n (k)$ are the Gregory coefficients of higher order with $G n (1) = G n$ , $Γ$ is the gamma function and $ζ$ is the Hurwitz zeta function.^[7]^[6] Similar series with the Cauchy numbers of the second kind $C n$ reads^[7]^[6]

\Psi (v)=\ln(v-1)+\sum _{n=1}^{\infty }{\frac {C_{n}(n-1)!}{(v)_{n}}},\qquad \Re (v)>1,

A series with the Bernoulli polynomials of the second kind has the following form^[6]

\Psi (v)=\ln(v+a)+\sum _{n=1}^{\infty }{\frac {(-1)^{n}\psi _{n}(a)\,(n-1)!}{(v)_{n}}},\qquad \Re (v)>-a,

where $ψ n (a)$ are the Bernoulli polynomials of the second kind defined by the generating equation

{\frac {z(1+z)^{a}}{\ln(1+z)}}=\sum _{n=0}^{\infty }z^{n}\psi _{n}(a)\,,\qquad |z|<1\,,

It may be generalized to

\Psi (v)={\frac {1}{r}}\sum _{l=0}^{r-1}\ln(v+a+l)+{\frac {1}{r}}\sum _{n=1}^{\infty }{\frac {(-1)^{n}N_{n,r}(a)(n-1)!}{(v)_{n}}},\qquad \Re (v)>-a,\quad r=1,2,3,\ldots

where the polynomials $N n,r (a)$ are given by the following generating equation

{\frac {(1+z)^{a+m}-(1+z)^{a}}{\ln(1+z)}}=\sum _{n=0}^{\infty }N_{n,m}(a)z^{n},\qquad |z|<1,

so that $N n,1 (a) = ψ n (a)$ .^[6] Similar expressions with the logarithm of the gamma function involve these formulas^[6]

\Psi (v)={\frac {1}{v+a-{\tfrac {1}{2}}}}\left\{\ln \Gamma (v+a)+v-{\frac {1}{2}}\ln 2\pi -{\frac {1}{2}}+\sum _{n=1}^{\infty }{\frac {(-1)^{n}\psi _{n+1}(a)}{(v)_{n}}}(n-1)!\right\},\qquad \Re (v)>-a,

and

\Psi (v)={\frac {1}{{\tfrac {1}{2}}r+v+a-1}}\left\{\ln \Gamma (v+a)+v-{\frac {1}{2}}\ln 2\pi -{\frac {1}{2}}+{\frac {1}{r}}\sum _{n=0}^{r-2}(r-n-1)\ln(v+a+n)+{\frac {1}{r}}\sum _{n=1}^{\infty }{\frac {(-1)^{n}N_{n+1,r}(a)}{(v)_{n}}}(n-1)!\right\},\qquad \Re (v)>-a,\quad r=2,3,4,\ldots

Reflection formula

The digamma function satisfies a reflection formula similar to that of the gamma function:

\psi (1-x)-\psi (x)=\pi \cot \pi x

Recurrence formula and characterization

The digamma function satisfies the recurrence relation

\psi (x+1)=\psi (x)+{\frac {1}{x}}.

Thus, it can be said to "telescope" $1 / x$ , for one has

\Delta [\psi ](x)={\frac {1}{x}}

where $Δ$ is the forward difference operator. This satisfies the recurrence relation of a partial sum of the harmonic series, thus implying the formula

\psi (n)=H_{n-1}-\gamma

where $γ$ is the Euler–Mascheroni constant.

More generally, one has

\psi (1+z)=-\gamma +\sum _{k=1}^{\infty }\left({\frac {1}{k}}-{\frac {1}{z+k}}\right).

for $Re(z)>0$ . Another series expansion is:

\psi (1+z)=\ln(z)+{\frac {1}{2z}}-\displaystyle \sum _{j=1}^{\infty }{\frac {B_{2j}}{2jz^{2j}}}

,

where $B_{2j}$ are the Bernoulli numbers. This series diverges for all $z$ and is known as the Stirling series.

Actually, $ψ$ is the only solution of the functional equation

F(x+1)=F(x)+{\frac {1}{x}}

that is monotonic on $ℝ +$ and satisfies $F (1) = - γ$ . This fact follows immediately from the uniqueness of the $Γ$ function given its recurrence equation and convexity restriction. This implies the useful difference equation:

\psi (x+N)-\psi (x)=\sum _{k=0}^{N-1}{\frac {1}{x+k}}

Some finite sums involving the digamma function

There are numerous finite summation formulas for the digamma function. Basic summation formulas, such as

\sum _{r=1}^{m}\psi \left({\frac {r}{m}}\right)=-m(\gamma +\ln m),

\sum _{r=1}^{m}\psi \left({\frac {r}{m}}\right)\cdot \exp {\dfrac {2\pi rki}{m}}=m\ln \left(1-\exp {\frac {2\pi ki}{m}}\right),\qquad k\in \mathbb {Z} ,\quad m\in \mathbb {N} ,\ k\neq m.

\sum _{r=1}^{m-1}\psi \left({\frac {r}{m}}\right)\cdot \cos {\dfrac {2\pi rk}{m}}=m\ln \left(2\sin {\frac {k\pi }{m}}\right)+\gamma ,\qquad k=1,2,\ldots ,m-1

\sum _{r=1}^{m-1}\psi \left({\frac {r}{m}}\right)\cdot \sin {\frac {2\pi rk}{m}}={\frac {\pi }{2}}(2k-m),\qquad k=1,2,\ldots ,m-1

are due to Gauss.^[8]^[9] More complicated formulas, such as

\sum _{r=0}^{m-1}\psi \left({\frac {2r+1}{2m}}\right)\cdot \cos {\frac {(2r+1)k\pi }{m}}=m\ln \left(\tan {\frac {\pi k}{2m}}\right),\qquad k=1,2,\ldots ,m-1

\sum _{r=0}^{m-1}\psi \left({\frac {2r+1}{2m}}\right)\cdot \sin {\dfrac {(2r+1)k\pi }{m}}=-{\frac {\pi m}{2}},\qquad k=1,2,\ldots ,m-1

\sum _{r=1}^{m-1}\psi \left({\frac {r}{m}}\right)\cdot \cot {\frac {\pi r}{m}}=-{\frac {\pi (m-1)(m-2)}{6}}

\sum _{r=1}^{m-1}\psi \left({\frac {r}{m}}\right)\cdot {\frac {r}{m}}=-{\frac {\gamma }{2}}(m-1)-{\frac {m}{2}}\ln m-{\frac {\pi }{2}}\sum _{r=1}^{m-1}{\frac {r}{m}}\cdot \cot {\frac {\pi r}{m}}

\sum _{r=1}^{m-1}\psi \left({\frac {r}{m}}\right)\cdot \cos {\dfrac {(2\ell +1)\pi r}{m}}=-{\frac {\pi }{m}}\sum _{r=1}^{m-1}{\frac {r\cdot \sin {\dfrac {2\pi r}{m}}}{\cos {\dfrac {2\pi r}{m}}-\cos {\dfrac {(2\ell +1)\pi }{m}}}},\qquad \ell \in \mathbb {Z}

\sum _{r=1}^{m-1}\psi \left({\frac {r}{m}}\right)\cdot \sin {\dfrac {(2\ell +1)\pi r}{m}}=-(\gamma +\ln 2m)\cot {\frac {(2\ell +1)\pi }{2m}}+\sin {\dfrac {(2\ell +1)\pi }{m}}\sum _{r=1}^{m-1}{\frac {\ln \sin {\dfrac {\pi r}{m}}}{\cos {\dfrac {2\pi r}{m}}-\cos {\dfrac {(2\ell +1)\pi }{m}}}},\qquad \ell \in \mathbb {Z}

\sum _{r=1}^{m-1}\psi ^{2}\left({\frac {r}{m}}\right)=(m-1)\gamma ^{2}+m(2\gamma +\ln 4m)\ln {m}-m(m-1)\ln ^{2}2+{\frac {\pi ^{2}(m^{2}-3m+2)}{12}}+m\sum _{\ell =1}^{m-1}\ln ^{2}\sin {\frac {\pi \ell }{m}}

are due to works of certain modern authors (see e.g. Appendix B in Blagouchine (2014)^[10]).

Gauss's digamma theorem

For positive integers $r$ and $m$ ( $r < m$ ), the digamma function may be expressed in terms of Euler's constant and a finite number of elementary functions

\psi \left({\frac {r}{m}}\right)=-\gamma -\ln(2m)-{\frac {\pi }{2}}\cot \left({\frac {r\pi }{m}}\right)+2\sum _{n=1}^{\left\lfloor {\frac {m-1}{2}}\right\rfloor }\cos \left({\frac {2\pi nr}{m}}\right)\ln \sin \left({\frac {\pi n}{m}}\right)

which holds, because of its recurrence equation, for all rational arguments.

Computation and approximation

According to the Euler–Maclaurin formula applied to^[11]

\sum_{n=1}^x \frac{1}{n}

the digamma function for $x$ , a real number, can be approximated by

\psi (x)\approx \ln(x)-{\frac {1}{2x}}-{\frac {1}{12x^{2}}}+{\frac {1}{120x^{4}}}-{\frac {1}{252x^{6}}}+{\frac {1}{240x^{8}}}-{\frac {5}{660x^{10}}}+{\frac {691}{32760x^{12}}}-{\frac {1}{12x^{14}}}

which is the beginning of the asymptotical expansion of $ψ (x)$ . The full asymptotic series of this expansions is

\psi (x)\sim \ln(x)-{\frac {1}{2x}}+\sum _{n=1}^{\infty }{\frac {\zeta (1-2n)}{x^{2n}}}=\ln(x)-{\frac {1}{2x}}-\sum _{n=1}^{\infty }{\frac {B_{2n}}{2nx^{2n}}}

where $B k$ is the $k$ th Bernoulli number and $ζ$ is the Riemann zeta function. Although the infinite sum does not converge for any $x$ , this expansion becomes more accurate for larger values of $x$ and any finite partial sum cut off from the full series. To compute $ψ (x)$ for small $x$ , the recurrence relation

\psi(x+1) = \frac{1}{x} + \psi(x)

can be used to shift the value of $x$ to a higher value. Beal^[12] suggests using the above recurrence to shift $x$ to a value greater than 6 and then applying the above expansion with terms above $x 14$ cut off, which yields "more than enough precision" (at least 12 digits except near the zeroes).

As $x$ goes to infinity, $ψ (x)$ gets arbitrarily close to both $ln(x - 1/2)$ and $ln x$ . Going down from $x + 1$ to $x$ , $ψ$ decreases by $1 / x$ , $ln(x - 1/2)$ decreases by $ln (x + 1/2) / (x - 1/2)$ , which is more than $1 / x$ , and $ln x$ decreases by $ln (1 + 1 / x)$ , which is less than $1 / x$ . From this we see that for any positive $x$ greater than $1/2$ ,

\psi (x)\in \left(\ln \left(x-{\tfrac {1}{2}}\right),\ln x\right)

or, for any positive $x$ ,

\exp \psi (x)\in \left(x-{\tfrac {1}{2}},x\right).

The exponential $exp ψ (x)$ is approximately $x - 1/2$ for large $x$ , but gets closer to $x$ at small $x$ , approaching 0 at $x = 0$ .

For $x < 1$ , we can calculate limits based on the fact that between 1 and 2, $ψ (x) \in [- γ, 1 - γ]$ , so

\psi (x)\in \left(-{\frac {1}{x}}-\gamma ,1-{\frac {1}{x}}-\gamma \right),\quad x\in (0,1)

or

\exp \psi (x)\in \left(\exp \left(-{\frac {1}{x}}-\gamma \right),e\exp \left(-{\frac {1}{x}}-\gamma \right)\right).

From the above asymptotic series for $ψ$ , one can derive an asymptotic series for $exp(- ψ (x))$ . The series matches the overall behaviour well, that is, it behaves asymptotically as it should for large arguments, and has a zero of unbounded multiplicity at the origin too.

{\frac {1}{\exp \psi (x)}}\sim {\frac {1}{x}}+{\frac {1}{2\cdot x^{2}}}+{\frac {5}{4\cdot 3!\cdot x^{3}}}+{\frac {3}{2\cdot 4!\cdot x^{4}}}+{\frac {47}{48\cdot 5!\cdot x^{5}}}-{\frac {5}{16\cdot 6!\cdot x^{6}}}+\cdots

This is similar to a Taylor expansion of $exp(- ψ (1 / y))$ at $y = 0$ , but it does not converge.^[13] (The function is not analytic at infinity.) A similar series exists for $exp(ψ (x))$ which starts with $\exp \psi (x)\sim x-{\frac {1}{2}}.$

If one calculates the asymptotic series for $ψ (x +1/2)$ it turns out that there are no odd powers of $x$ (there is no $x$ ⁻¹ term). This leads to the following asymptotic expansion, which saves computing terms of even order.

\exp \psi \left(x+{\tfrac {1}{2}}\right)\sim x+{\frac {1}{4!\cdot x}}-{\frac {37}{8\cdot 6!\cdot x^{3}}}+{\frac {10313}{72\cdot 8!\cdot x^{5}}}-{\frac {5509121}{384\cdot 10!\cdot x^{7}}}+\cdots

Special values

The digamma function has values in closed form for rational numbers, as a result of Gauss's digamma theorem. Some are listed below:

{\begin{aligned}\psi (1)&=-\gamma \\\psi \left({\tfrac {1}{2}}\right)&=-2\ln {2}-\gamma \\\psi \left({\tfrac {1}{3}}\right)&=-{\frac {\pi }{2{\sqrt {3}}}}-{\frac {3\ln {3}}{2}}-\gamma \\\psi \left({\tfrac {1}{4}}\right)&=-{\frac {\pi }{2}}-3\ln {2}-\gamma \\\psi \left({\tfrac {1}{6}}\right)&=-{\frac {\pi {\sqrt {3}}}{2}}-2\ln {2}-{\frac {3\ln {3}}{2}}-\gamma \\\psi \left({\tfrac {1}{8}}\right)&=-{\frac {\pi }{2}}-4\ln {2}-{\frac {\pi +\ln \left(2+{\sqrt {2}}\right)-\ln \left(2-{\sqrt {2}}\right)}{\sqrt {2}}}-\gamma .\end{aligned}}

Moreover, by the series representation, it can easily be deduced that at the imaginary unit,

{\begin{aligned}\operatorname {Re} \psi (i)&=-\gamma -\sum _{n=0}^{\infty }{\frac {n-1}{n^{3}+n^{2}+n+1}},\\[8px]\operatorname {Im} \psi (i)&=\sum _{n=0}^{\infty }{\frac {1}{n^{2}+1}}\\&={\frac {1}{2}}+{\frac {\pi }{2}}\coth \pi .\end{aligned}}

Roots of the digamma function

The roots of the digamma function are the saddle points of the complex-valued gamma function. Thus they lie all on the real axis. The only one on the positive real axis is the unique minimum of the real-valued gamma function on $ℝ +$ at $x 0 = 7000146163214496800♠ 1.461 632144968 ...$ . All others occur single between the poles on the negative axis:

{\begin{aligned}x_{1}&=-0.504\,083\,008\ldots ,\\x_{2}&=-1.573\,498\,473\ldots ,\\x_{3}&=-2.610\,720\,868\ldots ,\\x_{4}&=-3.635\,293\,366\ldots ,\\&\qquad \vdots \end{aligned}}

Already in 1881, Charles Hermite observed^[14] that

x_{n}=-n+{\frac {1}{\ln n}}+O\left({\frac {1}{(\ln n)^{2}}}\right)

holds asymptotically. A better approximation of the location of the roots is given by

x_n \approx -n + \frac{1}{\pi}\arctan\left(\frac{\pi}{\ln n}\right)\qquad n \ge 2

and using a further term it becomes still better

x_n \approx -n + \frac{1}{\pi}\arctan\left(\frac{\pi}{\ln n + \frac{1}{8n}}\right)\qquad n \ge 1

which both spring off the reflection formula via

0=\psi (1-x_{n})=\psi (x_{n})+{\frac {\pi }{\tan \pi x_{n}}}

and substituting $ψ (x n)$ by its not convergent asymptotic expansion. The correct second term of this expansion is $1 / 2 n$ , where the given one works good to approximate roots with small $n$ .

Another improvement of Hermite's formula can be given:^[4]

x_{n}=-n+{\frac {1}{\log n}}-{\frac {1}{2n(\log n)^{2}}}+O\left({\frac {1}{n^{2}(\log n)^{2}}}\right).

Regarding the zeros, the following infinite sum identities were recently proved by István Mező and Michael Hoffman^[4]

{\begin{aligned}\sum _{n=0}^{\infty }{\frac {1}{x_{n}^{2}}}&=\gamma ^{2}+{\frac {\pi ^{2}}{2}},\\\sum _{n=0}^{\infty }{\frac {1}{x_{n}^{3}}}&=-4\zeta (3)-\gamma ^{3}-{\frac {\gamma \pi ^{2}}{2}},\\\sum _{n=0}^{\infty }{\frac {1}{x_{n}^{4}}}&=\gamma ^{4}+{\frac {\pi ^{4}}{9}}+{\frac {2}{3}}\gamma ^{2}\pi ^{2}+4\gamma \zeta (3).\end{aligned}}

In general, the function

Z(k)=\sum _{n=0}^{\infty }{\frac {1}{x_{n}^{k}}}

can be determined and it is studied in detail by the cited authors.

The following results^[4]

{\begin{aligned}\sum _{n=0}^{\infty }{\frac {1}{x_{n}^{2}+x_{n}}}&=-2,\\\sum _{n=0}^{\infty }{\frac {1}{x_{n}^{2}-x_{n}}}&=\gamma +{\frac {\pi ^{2}}{6\gamma }}\end{aligned}}

also hold true.

Here $γ$ is the Euler–Mascheroni constant.

Regularization

The digamma function appears in the regularization of divergent integrals

\int _{0}^{\infty }{\frac {dx}{x+a}},

this integral can be approximated by a divergent general Harmonic series, but the following value can be attached to the series

\sum _{n=0}^{\infty }{\frac {1}{n+a}}=-\psi (a).

References

1 2 Abramowitz, M.; Stegun, I. A., eds. (1972). "6.3 psi (Digamma) Function.". Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (10th ed.). New York: Dover. pp. 258–259.
↑ Weisstein, Eric W. "Digamma function". MathWorld.
↑ Alzer, Horst; Jameson, Graham (2017). "A harmonic mean inequality for the digamma function and related results" (PDF). Rendiconti del Seminario Matematico della Università di Padova. 70 (201): 203–209. doi:10.4171/RSMUP/137-10. ISSN 0041-8994. LCCN 50046633. OCLC 01761704.
1 2 3 4 Mező, István; Hoffman, Michael E. (2017). "Zeros of the digamma function and its Barnes G-function analogue". Integral Transforms and Special Functions. 28 (28): 846–858. doi:10.1080/10652469.2017.1376193.
↑ Nörlund, N. E. (1924). Vorlesungen über Differenzenrechnung. Berlin: Springer.
1 2 3 4 5 6 7 Blagouchine, Ia. V. (2018). "Three Notes on Ser's and Hasse's Representations for the Zeta-functions" (PDF). Integers (Electronic Journal of Combinatorial Number Theory). 18A: 1–45. arXiv:1606.02044. Bibcode:2016arXiv160602044B.
1 2 Blagouchine, Ia. V. (2016). "Two series expansions for the logarithm of the gamma function involving Stirling numbers and containing only rational coefficients for certain arguments related to $π -1$ ". Journal of Mathematical Analysis and Applications. 442: 404–434. arXiv:1408.3902. Bibcode:2014arXiv1408.3902B.
↑ R. Campbell. Les intégrales eulériennes et leurs applications, Dunod, Paris, 1966.
↑ H.M. Srivastava and J. Choi. Series Associated with the Zeta and Related Functions, Kluwer Academic Publishers, the Netherlands, 2001.
↑ Blagouchine, Iaroslav V. (2014). "A theorem for the closed-form evaluation of the first generalized Stieltjes constant at rational arguments and some related summations". Journal of Number Theory. 148: 537–592. arXiv:1401.3724. doi:10.1016/j.jnt.2014.08.009.
↑ Bernardo, José M. (1976). "Algorithm AS 103 psi(digamma function) computation" (PDF). Applied Statistics. 25: 315–317.
↑ Beal, Matthew J. (2003). Variational Algorithms for Approximate Bayesian Inference (PDF) (PhD thesis). The Gatsby Computational Neuroscience Unit, University College London. pp. 265–266.
↑ If it converged to a function $f (y)$ then $ln(f (y) / y)$ would have the same Maclaurin series as $ln(1 / y) - φ (1 / y)$ . But this does not converge because the series given earlier for $φ (x)$ does not converge.
↑ Hermite, Charles (1881). "Sur l'intégrale Eulérienne de seconde espéce,". Journal für die reine und angewandte Mathematik (90): 332–338.

External links

OEIS sequence A020759 (Decimal expansion of (-1)*Gamma'(1/2)/Gamma(1/2) where Gamma(x) denotes the Gamma function)—psi(1/2)

A047787 psi(1/3),

A200064 psi(2/3),

A020777 psi(1/4),

A200134 psi(3/4),

A200135 to

A200138 psi(1/5) to psi(4/5).

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.

[AbramowitzStegun-1] 1 2 Abramowitz, M.; Stegun, I. A., eds. (1972). "6.3 psi (Digamma) Function.". Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables (10th ed.). New York: Dover. pp. 258–259.

[Weissstein-2] Weisstein, Eric W. "Digamma function". MathWorld.

[3] Alzer, Horst; Jameson, Graham (2017). "A harmonic mean inequality for the digamma function and related results" (PDF). Rendiconti del Seminario Matematico della Università di Padova. 70 (201): 203–209. doi:10.4171/RSMUP/137-10. ISSN 0041-8994. LCCN 50046633. OCLC 01761704.

[MezoHoffman-4] 1 2 3 4 Mező, István; Hoffman, Michael E. (2017). "Zeros of the digamma function and its Barnes G-function analogue". Integral Transforms and Special Functions. 28 (28): 846–858. doi:10.1080/10652469.2017.1376193.

[5] Nörlund, N. E. (1924). Vorlesungen über Differenzenrechnung. Berlin: Springer.

[blag2018-6] 1 2 3 4 5 6 7 Blagouchine, Ia. V. (2018). "Three Notes on Ser's and Hasse's Representations for the Zeta-functions" (PDF). Integers (Electronic Journal of Combinatorial Number Theory). 18A: 1–45. arXiv:1606.02044. Bibcode:2016arXiv160602044B.

[blag2016-7] 1 2 Blagouchine, Ia. V. (2016). "Two series expansions for the logarithm of the gamma function involving Stirling numbers and containing only rational coefficients for certain arguments related to $π -1$ ". Journal of Mathematical Analysis and Applications. 442: 404–434. arXiv:1408.3902. Bibcode:2014arXiv1408.3902B.

[8] R. Campbell. Les intégrales eulériennes et leurs applications, Dunod, Paris, 1966.

[9] H.M. Srivastava and J. Choi. Series Associated with the Zeta and Related Functions, Kluwer Academic Publishers, the Netherlands, 2001.

[iaroslav_07-10] Blagouchine, Iaroslav V. (2014). "A theorem for the closed-form evaluation of the first generalized Stieltjes constant at rational arguments and some related summations". Journal of Number Theory. 148: 537–592. arXiv:1401.3724. doi:10.1016/j.jnt.2014.08.009.

[11] Bernardo, José M. (1976). "Algorithm AS 103 psi(digamma function) computation" (PDF). Applied Statistics. 25: 315–317.

[12] Beal, Matthew J. (2003). Variational Algorithms for Approximate Bayesian Inference (PDF) (PhD thesis). The Gatsby Computational Neuroscience Unit, University College London. pp. 265–266.

[13] If it converged to a function $f (y)$ then $ln(f (y) / y)$ would have the same Maclaurin series as $ln(1 / y) - φ (1 / y)$ . But this does not converge because the series given earlier for $φ (x)$ does not converge.

[Hermite-14] Hermite, Charles (1881). "Sur l'intégrale Eulérienne de seconde espéce,". Journal für die reine und angewandte Mathematik (90): 332–338.