Series and integrals for the digamma function
TheoremSeries for the derivative of \(\Gamma'(z+1)/\Gamma(z+1)\)#✓fully proved
Let \(\gamma\) be the Euler–Mascheroni constant and let \(\varphi(w) = \gamma + \Gamma'(w+1)/\Gamma(w+1)\) for complex \(w\). For every complex number \(z\) other than \(-1, -2, -3, \dots\), the function \(\varphi\) is complex differentiable at \(z\), and the series \(\sum_{k=1}^{\infty} \frac{1}{(z+k)^2}\) converges absolutely with sum \(\varphi'(z)\).
theorem ramanujan_part1_ch6_entry2_stirlingseries (z : ℂ)
(hz : ∀ k : ℕ, z ≠ -((k + 1 : ℕ) : ℂ)) :
DifferentiableAt ℂ chapter6HarmonicInterpolation z ∧
HasSum (fun k : ℕ => 1 / (z + (k + 1 : ℕ)) ^ 2)
(deriv chapter6HarmonicInterpolation z) ∧
Summable (fun k : ℕ => ‖(1 / (z + (k + 1 : ℕ)) ^ 2 : ℂ)‖)
def chapter6HarmonicInterpolation (z : ℂ) : ℂ :=
(Real.eulerMascheroniConstant : ℂ) +
deriv Complex.Gamma (z + 1) / Complex.Gamma (z + 1)
MathlibExt/Analysis/Ramanujan/Part1Ch6Entry2Stirlingseries.lean:49
Sources
- Bruce C. Berndt, Ramanujan's Notebooks, Part I, Chapter 6.
TheoremRamanujan's integral for \(\psi(1 - 1/x)\) (Chapter 8, Entry 10)#✓fully proved
Let \(x\) be a complex number with \(\operatorname{Re} x < 0\). For real \(u\) put \(u^x = e^{x \log |u|}\) when \(u \ne 0\), and define \(K_x(0) = 0\), \(K_x(1) = -1/x\) and \(K_x(u) = \frac{1 - u}{u\,(u^x - 1)}\) for every other real \(u\), the quotient being taken to be \(0\) at \(u = -1\), where its denominator vanishes. Then \(x \ne 0\) and \(\Gamma(1 - 1/x) \ne 0\); for every \(u \in (0, 1)\) the denominator \(u\,(u^x - 1)\) is nonzero and \(K_x(u) = \frac{1 - u}{u\,(u^x - 1)}\); \(K_x\) is integrable on \([0, 1]\); and
where \(\psi\) is the digamma function and \(\gamma\) is the Euler–Mascheroni constant.
theorem ramanujan_part1_ch8_entry10_qgammalimit (x : ℂ) (hx : x.re < 0) :
x ≠ 0 ∧
Complex.Gamma (1 - 1 / x) ≠ 0 ∧
chapter8Entry10Kernel x 0 = 0 ∧
chapter8Entry10Kernel x 1 = -1 / x ∧
(∀ u : ℝ, u ∈ Set.Ioo 0 1 →
(u : ℂ) * (Entry9Qgamma.chapter8UnitPow x u - 1) ≠ 0 ∧
chapter8Entry10Kernel x u =
((1 - u : ℝ) : ℂ) /
((u : ℂ) * (Entry9Qgamma.chapter8UnitPow x u - 1))) ∧
IntervalIntegrable (chapter8Entry10Kernel x) volume 0 1 ∧
Complex.digamma (1 - 1 / x) +
(Real.eulerMascheroniConstant : ℂ) =
-x * ∫ u in (0 : ℝ)..1, chapter8Entry10Kernel x u
def chapter8Entry10Kernel (x : ℂ) (u : ℝ) : ℂ :=
if u = 0 then 0
else if u = 1 then -1 / x
else
((1 - u : ℝ) : ℂ) /
((u : ℂ) * (Entry9Qgamma.chapter8UnitPow x u - 1))
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry10Qgammalimit.lean:76
def chapter8UnitPow (x : ℂ) (u : ℝ) : ℂ :=
if u = 0 then 0 else Complex.exp (x * Real.log u)
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry9Qgamma.lean:81
Sources
- Bruce C. Berndt, Ramanujan's Notebooks, Part I, Chapter 8.
TheoremDigamma as the logarithmic derivative of Gamma (used for Chapter 8, Entry 11)#✓fully proved
For every complex number \(z\), \(\Gamma'(z)/\Gamma(z) = \psi(z)\), where \(\psi\) is the digamma function.
At the poles \(z = 0, -1, -2, \dots\) the conventions \(\Gamma(z) = 0\) and division by zero giving \(0\) make both sides \(0\), so the identity holds for every \(z\).
theorem chapter8Digamma_eq_digamma (z : ℂ) :
chapter8Digamma z = Complex.digamma z
def chapter8Digamma (z : ℂ) : ℂ :=
deriv Complex.Gamma z / Complex.Gamma z
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry11Qgammafunctional.lean:79
TheoremThe Entry 11 and Entry 12 series terms agree#✓fully proved
For every complex number \(x\) and every natural number \(j\), with \(k = j + 1\), the term \(\frac{2}{(kx)^3 - kx}\) as defined for Entry 11 equals the term given by the same formula for Entry 12.
In both definitions a term with zero denominator is \(0\).
theorem chapter8PhiTerm_eq_entry12Qgammaintegers (x : ℂ) (j : ℕ) :
chapter8PhiTerm x j =
MathlibExt.Analysis.Ramanujan.Part1Ch8.Entry12Qgammaintegers.chapter8PhiTerm x j
def chapter8PhiTerm (x : ℂ) (j : ℕ) : ℂ :=
let k := j + 1
2 / ((((k : ℂ) * x) ^ 3) - (k : ℂ) * x)
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry11Qgammafunctional.lean:82
def chapter8PhiTerm (x : ℂ) (j : ℕ) : ℂ :=
let k := j + 1
2 / ((((k : ℂ) * x) ^ 3) - (k : ℂ) * x)
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry12Qgammaintegers.lean:56
Used by 1
TheoremThe Entry 11 and Entry 12 definitions of \(\varphi\) agree#✓fully proved
For every complex number \(x\), the function \(\varphi(x) = 1 + \sum_{k \ge 1} \frac{2}{(kx)^3 - kx}\) as defined for Entry 11 equals the function given by the same formula for Entry 12.
In both definitions the sum is the unconditional sum, taken to be \(0\) when it does not converge, and a term with zero denominator is \(0\).
theorem chapter8Phi_eq_entry12Qgammaintegers (x : ℂ) :
chapter8Phi x =
MathlibExt.Analysis.Ramanujan.Part1Ch8.Entry12Qgammaintegers.chapter8Phi x
def chapter8Phi (x : ℂ) : ℂ :=
1 + ∑' j : ℕ, chapter8PhiTerm x j
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry11Qgammafunctional.lean:86
def chapter8Phi (x : ℂ) : ℂ :=
1 + ∑' j : ℕ, chapter8PhiTerm x j
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry12Qgammaintegers.lean:60
Uses 1
TheoremChapter 8, Entry 11 in terms of the digamma function#✓fully proved
Let \(x\) be a nonzero complex number with \(x \ne 1/k\) and \(x \ne -1/k\) for every positive integer \(k\), and put \(\varphi(x) = 1 + \sum_{k \ge 1} \frac{2}{(kx)^3 - kx}\). Then the series \(\sum_{k \ge 1} \frac{2}{(kx)^3 - kx}\) converges unconditionally, \(\Gamma(1/x) \ne 0\), \(\Gamma(1 - 1/x) \ne 0\), and
where \(\psi\) is the digamma function and \(\gamma\) is the Euler–Mascheroni constant.
theorem ramanujan_part1_ch8_entry11_qgammafunctional_digamma
(x : ℂ) (hx0 : x ≠ 0)
(hpoles : ∀ k : ℕ,
x ≠ 1 / ((k + 1 : ℕ) : ℂ) ∧
x ≠ -(1 / ((k + 1 : ℕ) : ℂ))) :
Summable (chapter8PhiTerm x) ∧
Complex.Gamma (1 / x) ≠ 0 ∧
Complex.Gamma (1 - 1 / x) ≠ 0 ∧
Complex.digamma (1 / x) + Complex.digamma (1 - 1 / x) =
-2 * (Real.eulerMascheroniConstant : ℂ) - x * chapter8Phi x
def chapter8PhiTerm (x : ℂ) (j : ℕ) : ℂ :=
let k := j + 1
2 / ((((k : ℂ) * x) ^ 3) - (k : ℂ) * x)
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry11Qgammafunctional.lean:82
def chapter8Phi (x : ℂ) : ℂ :=
1 + ∑' j : ℕ, chapter8PhiTerm x j
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry11Qgammafunctional.lean:86
TheoremRamanujan's identity for \(\Gamma'/\Gamma\) at \(1/x\) and \(1 - 1/x\) (Chapter 8, Entry 11)#✓fully proved
Let \(x\) be a nonzero complex number with \(x \ne 1/k\) and \(x \ne -1/k\) for every positive integer \(k\), and put \(\varphi(x) = 1 + \sum_{k \ge 1} \frac{2}{(kx)^3 - kx}\). Then the series \(\sum_{k \ge 1} \frac{2}{(kx)^3 - kx}\) converges unconditionally, \(\Gamma(1/x) \ne 0\), \(\Gamma(1 - 1/x) \ne 0\), and
where \(\gamma\) is the Euler–Mascheroni constant.
theorem ramanujan_part1_ch8_entry11_qgammafunctional
(x : ℂ) (hx0 : x ≠ 0)
(hpoles : ∀ k : ℕ,
x ≠ 1 / ((k + 1 : ℕ) : ℂ) ∧
x ≠ -(1 / ((k + 1 : ℕ) : ℂ))) :
Summable (chapter8PhiTerm x) ∧
Complex.Gamma (1 / x) ≠ 0 ∧
Complex.Gamma (1 - 1 / x) ≠ 0 ∧
chapter8Digamma (1 / x) + chapter8Digamma (1 - 1 / x) =
-2 * (Real.eulerMascheroniConstant : ℂ) - x * chapter8Phi x
def chapter8PhiTerm (x : ℂ) (j : ℕ) : ℂ :=
let k := j + 1
2 / ((((k : ℂ) * x) ^ 3) - (k : ℂ) * x)
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry11Qgammafunctional.lean:82
def chapter8Digamma (z : ℂ) : ℂ :=
deriv Complex.Gamma z / Complex.Gamma z
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry11Qgammafunctional.lean:79
def chapter8Phi (x : ℂ) : ℂ :=
1 + ∑' j : ℕ, chapter8PhiTerm x j
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry11Qgammafunctional.lean:86
Sources
- Bruce C. Berndt, Ramanujan's Notebooks, Part I, Chapter 8.
TheoremExponential and principal powers agree for positive bases (Chapter 8, Entry 12)#✓fully proved
For a complex number \(x\) and a real number \(u\), let \(P_x(u) = 0\) if \(u = 0\) and \(P_x(u) = e^{x \log |u|}\) otherwise. Then for every complex \(x\) and every real \(u > 0\), \(P_x(u) = u^x\), where \(u^x\) is the principal complex power.
theorem chapter8UnitPow_eq_cpow (x : ℂ) (u : ℝ) (hu : 0 < u) :
chapter8UnitPow x u = (u : ℂ) ^ x
def chapter8UnitPow (x : ℂ) (u : ℝ) : ℂ :=
if u = 0 then 0 else Complex.exp (x * Real.log u)
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry12Qgammaintegers.lean:46
TheoremRamanujan's integral for \(\varphi(x) - 1\) with principal powers (Chapter 8, Entry 12)#✓fully proved
Let \(x\) be a complex number with \(\operatorname{Re} x > 1\). For real \(u\) and complex \(s\) put \(P_s(u) = 0\) if \(u = 0\) and \(P_s(u) = e^{s \log |u|}\) otherwise, and define \(K_x(0) = 0\), \(K_x(1) = 0\) and \(K_x(u) = \frac{P_{x-2}(u)\,(1 - u)^2}{1 - P_x(u)}\) for every other real \(u\), the quotient being taken to be \(0\) at \(u = -1\), where its denominator vanishes. Let \(\varphi(x) = 1 + \sum_{k \ge 1} \frac{2}{(kx)^3 - kx}\). Then \((kx)^3 - kx \ne 0\) for every positive integer \(k\); the series \(\sum_{k \ge 1} \frac{2}{(kx)^3 - kx}\) converges unconditionally; \(K_x(0) = 0\) and \(K_x(1) = 0\); for every \(u \in (0, 1)\), with \(u^x\) and \(u^{x-2}\) the principal complex powers, \(1 - u^x \ne 0\) and \(K_x(u) = \frac{u^{x-2}\,(1 - u)^2}{1 - u^x}\); \(K_x\) is integrable on \([0, 1]\); and
For \(u > 0\) the powers \(P_s(u)\) and \(u^s\) agree, so this is the same statement as the version written with \(P_s\).
MathlibExt.Analysis.Ramanujan.Part1Ch8.Entry12Qgammaintegers.ramanujan_part1_ch8_entry12_qgammaintegers_cpow · 15 lines
theorem ramanujan_part1_ch8_entry12_qgammaintegers_cpow (x : ℂ) (hx : 1 < x.re) :
(∀ j : ℕ,
let k : ℂ := j + 1
(k * x) ^ 3 - k * x ≠ 0) ∧
Summable (chapter8PhiTerm x) ∧
chapter8Entry12Kernel x 0 = 0 ∧
chapter8Entry12Kernel x 1 = 0 ∧
(∀ u : ℝ, u ∈ Set.Ioo 0 1 →
1 - (u : ℂ) ^ x ≠ 0 ∧
chapter8Entry12Kernel x u =
(u : ℂ) ^ (x - 2) * ((1 - u : ℝ) : ℂ) ^ 2 /
(1 - (u : ℂ) ^ x)) ∧
IntervalIntegrable (chapter8Entry12Kernel x) volume 0 1 ∧
(∫ u in (0 : ℝ)..1, chapter8Entry12Kernel x u) =
chapter8Phi x - 1
def chapter8PhiTerm (x : ℂ) (j : ℕ) : ℂ :=
let k := j + 1
2 / ((((k : ℂ) * x) ^ 3) - (k : ℂ) * x)
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry12Qgammaintegers.lean:56
def chapter8Entry12Kernel (x : ℂ) (u : ℝ) : ℂ :=
if u = 0 then 0
else if u = 1 then 0
else
chapter8UnitPow (x - 2) u * ((1 - u : ℝ) : ℂ) ^ 2 /
(1 - chapter8UnitPow x u)
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry12Qgammaintegers.lean:49
def chapter8Phi (x : ℂ) : ℂ :=
1 + ∑' j : ℕ, chapter8PhiTerm x j
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry12Qgammaintegers.lean:60
Uses 1 · Used by 1
TheoremRamanujan's integral for \(\varphi(x) - 1\) (Chapter 8, Entry 12)#✓fully proved
Let \(x\) be a complex number with \(\operatorname{Re} x > 1\). For real \(u\) and complex \(s\) put \(P_s(u) = 0\) if \(u = 0\) and \(P_s(u) = e^{s \log |u|}\) otherwise, and define \(K_x(0) = 0\), \(K_x(1) = 0\) and \(K_x(u) = \frac{P_{x-2}(u)\,(1 - u)^2}{1 - P_x(u)}\) for every other real \(u\), the quotient being taken to be \(0\) at \(u = -1\), where its denominator vanishes. Let \(\varphi(x) = 1 + \sum_{k \ge 1} \frac{2}{(kx)^3 - kx}\). Then \((kx)^3 - kx \ne 0\) for every positive integer \(k\); the series \(\sum_{k \ge 1} \frac{2}{(kx)^3 - kx}\) converges unconditionally; \(K_x(0) = 0\) and \(K_x(1) = 0\); for every \(u \in (0, 1)\), \(1 - P_x(u) \ne 0\) and \(K_x(u) = \frac{P_{x-2}(u)\,(1 - u)^2}{1 - P_x(u)}\); \(K_x\) is integrable on \([0, 1]\); and
MathlibExt.Analysis.Ramanujan.Part1Ch8.Entry12Qgammaintegers.ramanujan_part1_ch8_entry12_qgammaintegers · 15 lines
theorem ramanujan_part1_ch8_entry12_qgammaintegers (x : ℂ) (hx : 1 < x.re) :
(∀ j : ℕ,
let k : ℂ := j + 1
(k * x) ^ 3 - k * x ≠ 0) ∧
Summable (chapter8PhiTerm x) ∧
chapter8Entry12Kernel x 0 = 0 ∧
chapter8Entry12Kernel x 1 = 0 ∧
(∀ u : ℝ, u ∈ Set.Ioo 0 1 →
1 - chapter8UnitPow x u ≠ 0 ∧
chapter8Entry12Kernel x u =
chapter8UnitPow (x - 2) u * ((1 - u : ℝ) : ℂ) ^ 2 /
(1 - chapter8UnitPow x u)) ∧
IntervalIntegrable (chapter8Entry12Kernel x) volume 0 1 ∧
(∫ u in (0 : ℝ)..1, chapter8Entry12Kernel x u) =
chapter8Phi x - 1
def chapter8PhiTerm (x : ℂ) (j : ℕ) : ℂ :=
let k := j + 1
2 / ((((k : ℂ) * x) ^ 3) - (k : ℂ) * x)
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry12Qgammaintegers.lean:56
def chapter8Entry12Kernel (x : ℂ) (u : ℝ) : ℂ :=
if u = 0 then 0
else if u = 1 then 0
else
chapter8UnitPow (x - 2) u * ((1 - u : ℝ) : ℂ) ^ 2 /
(1 - chapter8UnitPow x u)
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry12Qgammaintegers.lean:49
def chapter8UnitPow (x : ℂ) (u : ℝ) : ℂ :=
if u = 0 then 0 else Complex.exp (x * Real.log u)
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry12Qgammaintegers.lean:46
def chapter8Phi (x : ℂ) : ℂ :=
1 + ∑' j : ℕ, chapter8PhiTerm x j
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry12Qgammaintegers.lean:60
Sources
- Bruce C. Berndt, Ramanujan's Notebooks, Part I, Chapter 8.
TheoremDigamma as the logarithmic derivative of Gamma (used for Chapter 8, Entry 8)#✓fully proved
For every complex number \(z\), \(\Gamma'(z)/\Gamma(z) = \psi(z)\), where \(\psi\) is the digamma function.
At the poles \(z = 0, -1, -2, \dots\) the conventions \(\Gamma(z) = 0\) and division by zero giving \(0\) make both sides \(0\), so the identity holds for every \(z\).
theorem chapter8Digamma_eq_Complex_digamma (z : ℂ) :
chapter8Digamma z = Complex.digamma z
def chapter8Digamma (z : ℂ) : ℂ :=
deriv Complex.Gamma z / Complex.Gamma z
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry8Digammaseries.lean:68
Used by 1
TheoremAlternating series \(\sum_{k \ge 1} (-1)^{k-1}/(a + kb)\) via the digamma function#✓fully proved
Let \(a\) and \(b\) be complex numbers with \(b \ne 0\) and \(a + kb \ne 0\) for every integer \(k \ge 1\), and let \(\psi\) be the digamma function, the logarithmic derivative of \(\Gamma\). Then the partial sums of \(\sum_{k=1}^{\infty} \frac{(-1)^{k-1}}{a + kb}\) converge to
theorem ramanujan_part1_ch8_entry8_digammaseries_complex_digamma
(a b : ℂ) (hb : b ≠ 0)
(hdenom : ∀ j : ℕ, a + b * ((j + 1 : ℕ) : ℂ) ≠ 0) :
HasSum (chapter8Entry8Term a b)
((Complex.digamma ((a + 2 * b) / (2 * b)) -
Complex.digamma ((a + b) / (2 * b))) /
(2 * b))
(SummationFilter.conditional ℕ)
def chapter8Entry8Term (a b : ℂ) (j : ℕ) : ℂ :=
(-1 : ℂ) ^ j / (a + b * ((j + 1 : ℕ) : ℂ))
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry8Digammaseries.lean:71
Uses 1 · Used by 1
TheoremRamanujan's alternating digamma series (Chapter 8, Entry 8)#✓fully proved
Let \(a\) and \(b\) be complex numbers with \(b \ne 0\) and \(a + kb \ne 0\) for every integer \(k \ge 1\), suppose that \(\Gamma\left(\frac{a+2b}{2b}\right) \ne 0\) and \(\Gamma\left(\frac{a+b}{2b}\right) \ne 0\), and let \(\psi(z) = \Gamma'(z)/\Gamma(z)\). Then the partial sums of \(\sum_{k=1}^{\infty} \frac{(-1)^{k-1}}{a + kb}\) converge to
Here \(\Gamma\) takes the value \(0\) at its poles \(0, -1, -2, \dots\), so the two conditions on \(\Gamma\) say that these arguments are not poles; they already follow from the condition on \(a + kb\).
theorem ramanujan_part1_ch8_entry8_digammaseries
(a b : ℂ) (hb : b ≠ 0)
(hleft : Complex.Gamma ((a + 2 * b) / (2 * b)) ≠ 0)
(hright : Complex.Gamma ((a + b) / (2 * b)) ≠ 0)
(hdenom : ∀ j : ℕ, a + b * ((j + 1 : ℕ) : ℂ) ≠ 0) :
HasSum (chapter8Entry8Term a b)
((chapter8Digamma ((a + 2 * b) / (2 * b)) -
chapter8Digamma ((a + b) / (2 * b))) /
(2 * b))
(SummationFilter.conditional ℕ)
def chapter8Entry8Term (a b : ℂ) (j : ℕ) : ℂ :=
(-1 : ℂ) ^ j / (a + b * ((j + 1 : ℕ) : ℂ))
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry8Digammaseries.lean:71
def chapter8Digamma (z : ℂ) : ℂ :=
deriv Complex.Gamma z / Complex.Gamma z
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry8Digammaseries.lean:68
Sources
- Bruce C. Berndt, Ramanujan's Notebooks, Part I, Chapter 8.
TheoremDigamma as the logarithmic derivative of Gamma (Chapter 8, Entry 9)#✓fully proved
For every complex number \(z\), \(\Gamma'(z)/\Gamma(z) = \psi(z)\), where \(\Gamma'\) is the complex derivative of the Gamma function and \(\psi\) is the digamma function. At the poles of \(\Gamma\) both sides are \(0\) by convention.
theorem chapter8Digamma_eq_digamma (z : ℂ) :
chapter8Digamma z = Complex.digamma z
def chapter8Digamma (z : ℂ) : ℂ :=
deriv Complex.Gamma z / Complex.Gamma z
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry9Qgamma.lean:78
Used by 1
TheoremThe Entry 9 power function on positive reals (Chapter 8)#✓fully proved
For a complex number \(x\) and a real number \(u\), let \(p_x(u) = e^{x \log\lvert u\rvert}\) if \(u \ne 0\) and \(p_x(0) = 0\). If \(u > 0\), then \(p_x(u) = u^x\), the principal complex power of \(u\).
theorem chapter8UnitPow_eq_cpow (x : ℂ) (u : ℝ) (hu : 0 < u) :
chapter8UnitPow x u = (u : ℂ) ^ x
def chapter8UnitPow (x : ℂ) (u : ℝ) : ℂ :=
if u = 0 then 0 else Complex.exp (x * Real.log u)
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry9Qgamma.lean:81
TheoremRamanujan's integral for the digamma function (Chapter 8, Entry 9)#✓fully proved
Let \(x\) be a complex number with \(\operatorname{Re} x > 0\), and for real \(u\) write \(u^x = e^{x \log u}\) if \(u > 0\) and \(0^x = 0\). Let \(\psi\) be the digamma function. Then \(x \ne 0\); the function \(u \mapsto u^x/(1 + u^x)\) is integrable on \([0, 1]\); \(1 + u^x \ne 0\) for every \(u \in [0, 1]\); \(\Gamma(1/(2x) + 1) \ne 0\) and \(\Gamma(1/x + 1) \ne 0\); and
theorem ramanujan_part1_ch8_entry9_qgamma_digamma (x : ℂ) (hx : 0 < x.re) :
x ≠ 0 ∧
IntervalIntegrable (chapter8UnitPowRatio x) volume 0 1 ∧
(∀ u : ℝ, u ∈ Set.Icc 0 1 →
1 + chapter8UnitPow x u ≠ 0) ∧
Complex.Gamma (1 / (2 * x) + 1) ≠ 0 ∧
Complex.Gamma (1 / x + 1) ≠ 0 ∧
Complex.digamma (1 / (2 * x) + 1) =
Complex.digamma (1 / x + 1) - (Real.log 2 : ℂ) +
x * ∫ u in (0 : ℝ)..1, chapter8UnitPowRatio x u
def chapter8UnitPowRatio (x : ℂ) (u : ℝ) : ℂ :=
chapter8UnitPow x u / (1 + chapter8UnitPow x u)
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry9Qgamma.lean:84
def chapter8UnitPow (x : ℂ) (u : ℝ) : ℂ :=
if u = 0 then 0 else Complex.exp (x * Real.log u)
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry9Qgamma.lean:81
Sources
- Bruce C. Berndt, Ramanujan's Notebooks, Part I, Chapter 8.
Used by 1
TheoremRamanujan's integral for \(\Gamma'/\Gamma\) (Chapter 8, Entry 9)#✓fully proved
Let \(x\) be a complex number with \(\operatorname{Re} x > 0\), and for real \(u\) write \(u^x = e^{x \log u}\) if \(u > 0\) and \(0^x = 0\). Write \(\Gamma'\) for the complex derivative of the Gamma function. Then \(x \ne 0\); the function \(u \mapsto u^x/(1 + u^x)\) is integrable on \([0, 1]\); \(1 + u^x \ne 0\) for every \(u \in [0, 1]\); \(\Gamma(1/(2x) + 1) \ne 0\) and \(\Gamma(1/x + 1) \ne 0\); and
theorem ramanujan_part1_ch8_entry9_qgamma (x : ℂ) (hx : 0 < x.re) :
x ≠ 0 ∧
IntervalIntegrable (chapter8UnitPowRatio x) volume 0 1 ∧
(∀ u : ℝ, u ∈ Set.Icc 0 1 →
1 + chapter8UnitPow x u ≠ 0) ∧
Complex.Gamma (1 / (2 * x) + 1) ≠ 0 ∧
Complex.Gamma (1 / x + 1) ≠ 0 ∧
chapter8Digamma (1 / (2 * x) + 1) =
chapter8Digamma (1 / x + 1) - (Real.log 2 : ℂ) +
x * ∫ u in (0 : ℝ)..1, chapter8UnitPowRatio x u
def chapter8UnitPowRatio (x : ℂ) (u : ℝ) : ℂ :=
chapter8UnitPow x u / (1 + chapter8UnitPow x u)
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry9Qgamma.lean:84
def chapter8UnitPow (x : ℂ) (u : ℝ) : ℂ :=
if u = 0 then 0 else Complex.exp (x * Real.log u)
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry9Qgamma.lean:81
def chapter8Digamma (z : ℂ) : ℂ :=
deriv Complex.Gamma z / Complex.Gamma z
MathlibExt/Analysis/Ramanujan/Part1Ch8Entry9Qgamma.lean:78
Sources
- Bruce C. Berndt, Ramanujan's Notebooks, Part I, Chapter 8.
Loading dependency explorer…
The dependency explorer needs JavaScript and its graph data.