The Gamma Function

Total Page:16

File Type:pdf, Size:1020Kb

The Gamma Function The Gamma Function EMILY GULLERUD, michael vaughan University of Wisconsin - Eau Claire May 10, 2017 EMILY GULLERUD, michael vaughan The Gamma Function May 10, 2017 1 / 15 What is the Gamma Function? Motivation The gamma function (Γ) is an extension of the factorial function to − C n Z [ 0. It is a solution to the interpolation problem of connecting the discrete points of the factorial function. EMILY GULLERUD, michael vaughan The Gamma Function May 10, 2017 2 / 15 What is the Gamma Function? Definitions Factorial function definition + If n 2 Z , then Γ(n) = (n − 1)! Improper integral definition R 1 z−1 −x If z 2 fx + iy j x > 0g, then Γ(z) = 0 x e dx. Analytic continuation is required to extend the integral definition to − z 2 fx + iy j y 6= 0 when x 2 Z [ 0g. This function has simple poles at all the non-positive integers. EMILY GULLERUD, michael vaughan The Gamma Function May 10, 2017 3 / 15 What is the Gamma Function? More Definitions Definition as an infinite product 1 1 z 1 Y 1 + Γ(z) = n z 1 + z n=1 n Weierstrass's definition 1 e−γz Y z −1 Γ(z) = 1 + ez=n z n n=1 EMILY GULLERUD, michael vaughan The Gamma Function May 10, 2017 4 / 15 Properties Identities and Formulas Requirements for Γ to be an extension of the factorial function: Γ(1) = 1 zΓ(z) = Γ(z + 1) for Re(z) > 0 Property of conjugation: Γ(z) = Γ(z) Complement Formula: π Γ(z)Γ(1 − z) = for z 2= sin(πz) Z Duplication Formula: 1 1−2z p Γ(z)Γ(z + 2 ) = Γ(2z) · 2 π EMILY GULLERUD, michael vaughan The Gamma Function May 10, 2017 5 / 15 Properties Proof of zΓ(z) = Γ(z + 1) Let α 2 C such that Re(α) > 0 Z 1 Γ(α) = xα−1e−x dx 0 1 Z 1 α−1 −x −x α−2 = x (−e ) − −e (α − 1)x dx 0 0 α−1 Z 1 x −x α−2 = − lim x + e (α − 1)x dx x!1 e 0 Z 1 = 0 + (α − 1) e−x xα−1−1dx 0 Γ(α) = (α − 1)Γ(α − 1) Take z = α − 1. Γ(z + 1) = zΓ(z) EMILY GULLERUD, michael vaughan The Gamma Function May 10, 2017 6 / 15 Properties Values of the Gamma Function 3 4p Γ − = π 2 3 1 p Γ − = −2 π 2 Γ(1) = 1 1 p Γ(2) = 1 Γ = π 2 Γ(3) = 2 3 1p Γ = π Γ(4) = 6 2 2 Γ(5) = 24 5 3p Γ = π 2 4 7 15p Γ = π 2 8 EMILY GULLERUD, michael vaughan The Gamma Function May 10, 2017 7 / 15 Properties p 1 π Example: Proving 2 ! = 2 1 3 ! = Γ 2 2 1 1 = Γ 2 2 1=2 Z 1 Let u = t 1 −1=2 −t 1 −1=2 = t e dt Then du = 2 t dt 2 0 1 Z 2 = e−u p0 π = 2 EMILY GULLERUD, michael vaughan The Gamma Function May 10, 2017 8 / 15 Applications and Connections to Other Functions The log-gamma function The Gamma function grows rapidly, so taking the natural logarithm yields a function which grows much more slowly: ln Γ(z) = ln Γ(z + 1) − ln z This function is used in many computing environments and in the context of wave propogation. EMILY GULLERUD, michael vaughan The Gamma Function May 10, 2017 9 / 15 Applications and Connections to Other Functions Derivative of the log-gamma function: Digamma function The Digamma function is defined to be the logarithmic derivative of the Gamma function: d Γ0 (z) (z) = ln (Γ(z)) = dz Γ(z) A general form of the Digamma function: the Polygamma function defined to be the (m + 1)th logarithmic derivative of the Gamma function: dm (m)(z) = ln (Γ(z)) dzm Notice that (0)(z) = (z) The Polygamma function is meromorphic on C (holomorphic on − C n Z [ 0 The nonpositive integers have poles of order m + 1 EMILY GULLERUD, michael vaughan The Gamma Function May 10, 2017 10 / 15 Applications and Connections to Other Functions Digamma and Polygamma functions EMILY GULLERUD, michael vaughan The Gamma Function May 10, 2017 11 / 15 Applications and Connections to Other Functions Digamma and Polygamma functions Finding a property of the Digamma function: Recall that Γ(z + 1) = zΓ(z) Take the derivative: 0 0 Γ (z + 1) = zΓ (z) + Γ(z) Divide by Γ(z + 1) = zΓ(z): Γ0 (z + 1) Γ0 (z) 1 = + Γ(z + 1) Γ(z) z Substitute in function: 1 (z + 1) = (z) + z EMILY GULLERUD, michael vaughan The Gamma Function May 10, 2017 12 / 15 Applications and Connections to Other Functions Incomplete Gamma Functions upper incomplete gamma function R 1 s−1 −t Γ(s; x) = x x e dt lower incomplete gamma function R x s−1 −t γ(s; x) = 0 x e dt EMILY GULLERUD, michael vaughan The Gamma Function May 10, 2017 13 / 15 Applications and Connections to Other Functions Pi and Beta Functions Definition of the Beta Function R 1 x−1 y−1 B(x; y) = 0 t (1 − t) dt for Re(x); Re(y) > 0 Beta function in terms of the Gamma function: Γ(x)Γ(y) B(x; y) = Γ(x + y) Definition of the Pi Function Π(z) = Γ(z + 1) = zΓ(z) EMILY GULLERUD, michael vaughan The Gamma Function May 10, 2017 14 / 15 Applications and Connections to Other Functions Riemann Zeta Function The functional equation: πs ζ(s) = 2s πs−1 sin Γ(1 − s)ζ(1 − s) 2 The functional equation in another form: − s s − 1−s 1 − s π 2 Γ ζ(s) = π 2 Γ ζ(1 − s) 2 2 Another relation: Z 1 uz−1 ζ(z)Γ(z) = u du for Re(z) > 1 0 e − 1 EMILY GULLERUD, michael vaughan The Gamma Function May 10, 2017 15 / 15.
Recommended publications
  • Introduction to Analytic Number Theory the Riemann Zeta Function and Its Functional Equation (And a Review of the Gamma Function and Poisson Summation)
    Math 229: Introduction to Analytic Number Theory The Riemann zeta function and its functional equation (and a review of the Gamma function and Poisson summation) Recall Euler’s identity: ∞ ∞ X Y X Y 1 [ζ(s) :=] n−s = p−cps = . (1) 1 − p−s n=1 p prime cp=0 p prime We showed that this holds as an identity between absolutely convergent sums and products for real s > 1. Riemann’s insight was to consider (1) as an identity between functions of a complex variable s. We follow the curious but nearly universal convention of writing the real and imaginary parts of s as σ and t, so s = σ + it. We already observed that for all real n > 0 we have |n−s| = n−σ, because n−s = exp(−s log n) = n−σe−it log n and e−it log n has absolute value 1; and that both sides of (1) converge absolutely in the half-plane σ > 1, and are equal there either by analytic continuation from the real ray t = 0 or by the same proof we used for the real case. Riemann showed that the function ζ(s) extends from that half-plane to a meromorphic function on all of C (the “Riemann zeta function”), analytic except for a simple pole at s = 1. The continuation to σ > 0 is readily obtained from our formula ∞ ∞ 1 X Z n+1 X Z n+1 ζ(s) − = n−s − x−s dx = (n−s − x−s) dx, s − 1 n=1 n n=1 n since for x ∈ [n, n + 1] (n ≥ 1) and σ > 0 we have Z x −s −s −1−s −1−σ |n − x | = s y dy ≤ |s|n n so the formula for ζ(s) − (1/(s − 1)) is a sum of analytic functions converging absolutely in compact subsets of {σ + it : σ > 0} and thus gives an analytic function there.
    [Show full text]
  • The Digamma Function and Explicit Permutations of the Alternating Harmonic Series
    The digamma function and explicit permutations of the alternating harmonic series. Maxim Gilula February 20, 2015 Abstract The main goal is to present a countable family of permutations of the natural numbers that provide explicit rearrangements of the alternating harmonic series and that we can easily define by some closed expression. The digamma function presents its ubiquity in mathematics once more by being the key tool in computing explicitly the simple rearrangements presented in this paper. The permutations are simple in the sense that composing one with itself will give the identity. We show that the count- able set of rearrangements presented are dense in the reals. Then, slight generalizations are presented. Finally, we reprove a result given originally by J.H. Smith in 1975 that for any conditionally convergent real series guarantees permutations of infinite cycle type give all rearrangements of the series [4]. This result provides a refinement of the well known theorem by Riemann (see e.g. Rudin [3] Theorem 3.54). 1 Introduction A permutation of order n of a conditionally convergent series is a bijection φ of the positive integers N with the property that φn = φ ◦ · · · ◦ φ is the identity on N and n is the least such. Given a conditionally convergent series, a nat- ural question to ask is whether for any real number L there is a permutation of order 2 (or n > 1) such that the rearrangement induced by the permutation equals L. This turns out to be an easy corollary of [4], and is reproved below with elementary methods. Other \simple"rearrangements have been considered elsewhere, such as in Stout [5] and the comprehensive references therein.
    [Show full text]
  • Analytic Continuation of Massless Two-Loop Four-Point Functions
    CERN-TH/2002-145 hep-ph/0207020 July 2002 Analytic Continuation of Massless Two-Loop Four-Point Functions T. Gehrmanna and E. Remiddib a Theory Division, CERN, CH-1211 Geneva 23, Switzerland b Dipartimento di Fisica, Universit`a di Bologna and INFN, Sezione di Bologna, I-40126 Bologna, Italy Abstract We describe the analytic continuation of two-loop four-point functions with one off-shell external leg and internal massless propagators from the Euclidean region of space-like 1 3decaytoMinkowskian regions relevant to all 1 3and2 2 reactions with one space-like or time-like! off-shell external leg. Our results can be used! to derive two-loop! master integrals and unrenormalized matrix elements for hadronic vector-boson-plus-jet production and deep inelastic two-plus-one-jet production, from results previously obtained for three-jet production in electron{positron annihilation. 1 Introduction In recent years, considerable progress has been made towards the extension of QCD calculations of jet observables towards the next-to-next-to-leading order (NNLO) in perturbation theory. One of the main ingredients in such calculations are the two-loop virtual corrections to the multi leg matrix elements relevant to jet physics, which describe either 1 3 decay or 2 2 scattering reactions: two-loop four-point functions with massless internal propagators and→ up to one off-shell→ external leg. Using dimensional regularization [1, 2] with d = 4 dimensions as regulator for ultraviolet and infrared divergences, the large number of different integrals6 appearing in the two-loop Feynman amplitudes for 2 2 scattering or 1 3 decay processes can be reduced to a small number of master integrals.
    [Show full text]
  • An Explicit Formula for Dirichlet's L-Function
    University of Tennessee at Chattanooga UTC Scholar Student Research, Creative Works, and Honors Theses Publications 5-2018 An explicit formula for Dirichlet's L-Function Shannon Michele Hyder University of Tennessee at Chattanooga, [email protected] Follow this and additional works at: https://scholar.utc.edu/honors-theses Part of the Mathematics Commons Recommended Citation Hyder, Shannon Michele, "An explicit formula for Dirichlet's L-Function" (2018). Honors Theses. This Theses is brought to you for free and open access by the Student Research, Creative Works, and Publications at UTC Scholar. It has been accepted for inclusion in Honors Theses by an authorized administrator of UTC Scholar. For more information, please contact [email protected]. An Explicit Formula for Dirichlet's L-Function Shannon M. Hyder Departmental Honors Thesis The University of Tennessee at Chattanooga Department of Mathematics Thesis Director: Dr. Andrew Ledoan Examination Date: April 9, 2018 Members of Examination Committee Dr. Andrew Ledoan Dr. Cuilan Gao Dr. Roger Nichols c 2018 Shannon M. Hyder ALL RIGHTS RESERVED i Abstract An Explicit Formula for Dirichlet's L-Function by Shannon M. Hyder The Riemann zeta function has a deep connection to the distribution of primes. In 1911 Landau proved that, for every fixed x > 1, X T xρ = − Λ(x) + O(log T ) 2π 0<γ≤T as T ! 1. Here ρ = β + iγ denotes a complex zero of the zeta function and Λ(x) is an extension of the usual von Mangoldt function, so that Λ(x) = log p if x is a positive integral power of a prime p and Λ(x) = 0 for all other real values of x.
    [Show full text]
  • The Riemann and Hurwitz Zeta Functions, Apery's Constant and New
    The Riemann and Hurwitz zeta functions, Apery’s constant and new rational series representations involving ζ(2k) Cezar Lupu1 1Department of Mathematics University of Pittsburgh Pittsburgh, PA, USA Algebra, Combinatorics and Geometry Graduate Student Research Seminar, February 2, 2017, Pittsburgh, PA A quick overview of the Riemann zeta function. The Riemann zeta function is defined by 1 X 1 ζ(s) = ; Re s > 1: ns n=1 Originally, Riemann zeta function was defined for real arguments. Also, Euler found another formula which relates the Riemann zeta function with prime numbrs, namely Y 1 ζ(s) = ; 1 p 1 − ps where p runs through all primes p = 2; 3; 5;:::. A quick overview of the Riemann zeta function. Moreover, Riemann proved that the following ζ(s) satisfies the following integral representation formula: 1 Z 1 us−1 ζ(s) = u du; Re s > 1; Γ(s) 0 e − 1 Z 1 where Γ(s) = ts−1e−t dt, Re s > 0 is the Euler gamma 0 function. Also, another important fact is that one can extend ζ(s) from Re s > 1 to Re s > 0. By an easy computation one has 1 X 1 (1 − 21−s )ζ(s) = (−1)n−1 ; ns n=1 and therefore we have A quick overview of the Riemann function. 1 1 X 1 ζ(s) = (−1)n−1 ; Re s > 0; s 6= 1: 1 − 21−s ns n=1 It is well-known that ζ is analytic and it has an analytic continuation at s = 1. At s = 1 it has a simple pole with residue 1.
    [Show full text]
  • INTEGRALS of POWERS of LOGGAMMA 1. Introduction The
    PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 139, Number 2, February 2011, Pages 535–545 S 0002-9939(2010)10589-0 Article electronically published on August 18, 2010 INTEGRALS OF POWERS OF LOGGAMMA TEWODROS AMDEBERHAN, MARK W. COFFEY, OLIVIER ESPINOSA, CHRISTOPH KOUTSCHAN, DANTE V. MANNA, AND VICTOR H. MOLL (Communicated by Ken Ono) Abstract. Properties of the integral of powers of log Γ(x) from 0 to 1 are con- sidered. Analytic evaluations for the first two powers are presented. Empirical evidence for the cubic case is discussed. 1. Introduction The evaluation of definite integrals is a subject full of interconnections of many parts of mathematics. Since the beginning of Integral Calculus, scientists have developed a large variety of techniques to produce magnificent formulae. A partic- ularly beautiful formula due to J. L. Raabe [12] is 1 Γ(x + t) (1.1) log √ dx = t log t − t, for t ≥ 0, 0 2π which includes the special case 1 √ (1.2) L1 := log Γ(x) dx =log 2π. 0 Here Γ(x)isthegamma function defined by the integral representation ∞ (1.3) Γ(x)= ux−1e−udu, 0 for Re x>0. Raabe’s formula can be obtained from the Hurwitz zeta function ∞ 1 (1.4) ζ(s, q)= (n + q)s n=0 via the integral formula 1 t1−s (1.5) ζ(s, q + t) dq = − 0 s 1 coupled with Lerch’s formula ∂ Γ(q) (1.6) ζ(s, q) =log √ . ∂s s=0 2π An interesting extension of these formulas to the p-adic gamma function has recently appeared in [3].
    [Show full text]
  • Bernstein's Analytic Continuation of Complex Powers of Polynomials
    Bernsteins analytic continuation of complex p owers c Paul Garrett garrettmathumnedu version January Analytic continuation of distributions Statement of the theorems on analytic continuation Bernsteins pro ofs Pro of of the Lemma the Bernstein p olynomial Pro of of the Prop osition estimates on zeros Garrett Bernsteins analytic continuation of complex p owers Let f b e a p olynomial in x x with real co ecients For complex s let n s f b e the function dened by s s f x f x if f x s f x if f x s Certainly for s the function f is lo cally integrable For s in this range s we can dened a distribution denoted by the same symb ol f by Z s s f x x dx f n R n R the space of compactlysupp orted smo oth realvalued where is in C c n functions on R s The ob ject is to analytically continue the distribution f as a meromorphic distributionvalued function of s This typ e of question was considered in several provo cative examples in IM Gelfand and GE Shilovs Generalized Functions volume I One should also ask ab out analytic continuation as a temp ered distribution In a lecture at the Amsterdam Congress IM Gelfand rened this question to require further that one show that the p oles lie in a nite numb er of arithmetic progressions Bernstein proved the result in under a certain regularity hyp othesis on the zeroset of the p olynomial f Published in Journal of Functional Analysis and Its Applications translated from Russian The present discussion includes some background material from complex function theory and
    [Show full text]
  • [Math.AG] 2 Jul 2015 Ic Oetime: Some and Since Time, of People
    MONODROMY AND NORMAL FORMS FABRIZIO CATANESE Abstract. We discuss the history of the monodromy theorem, starting from Weierstraß, and the concept of monodromy group. From this viewpoint we compare then the Weierstraß, the Le- gendre and other normal forms for elliptic curves, explaining their geometric meaning and distinguishing them by their stabilizer in PSL(2, Z) and their monodromy. Then we focus on the birth of the concept of the Jacobian variety, and the geometrization of the theory of Abelian functions and integrals. We end illustrating the methods of complex analysis in the simplest issue, the difference equation f(z)= g(z + 1) g(z) on C. − Contents Introduction 1 1. The monodromy theorem 3 1.1. Riemann domain and sheaves 6 1.2. Monodromy or polydromy? 6 2. Normalformsandmonodromy 8 3. Periodic functions and Abelian varieties 16 3.1. Cohomology as difference equations 21 References 23 Introduction In Jules Verne’s novel of 1874, ‘Le Tour du monde en quatre-vingts jours’ , Phileas Fogg is led to his remarkable adventure by a bet made arXiv:1507.00711v1 [math.AG] 2 Jul 2015 in his Club: is it possible to make a tour of the world in 80 days? Idle questions and bets can be very stimulating, but very difficult to answer when they deal with the history of mathematics, and one asks how certain ideas, which have been a common knowledge for long time, did indeed evolve and mature through a long period of time, and through the contributions of many people. In short, there are three idle questions which occupy my attention since some time: Date: July 3, 2015.
    [Show full text]
  • Chapter 4 the Riemann Zeta Function and L-Functions
    Chapter 4 The Riemann zeta function and L-functions 4.1 Basic facts We prove some results that will be used in the proof of the Prime Number Theorem (for arithmetic progressions). The L-function of a Dirichlet character χ modulo q is defined by 1 X L(s; χ) = χ(n)n−s: n=1 P1 −s We view ζ(s) = n=1 n as the L-function of the principal character modulo 1, (1) (1) more precisely, ζ(s) = L(s; χ0 ), where χ0 (n) = 1 for all n 2 Z. We first prove that ζ(s) has an analytic continuation to fs 2 C : Re s > 0gnf1g. We use an important summation formula, due to Euler. Lemma 4.1 (Euler's summation formula). Let a; b be integers with a < b and f :[a; b] ! C a continuously differentiable function. Then b X Z b Z b f(n) = f(x)dx + f(a) + (x − [x])f 0(x)dx: n=a a a 105 Remark. This result often occurs in the more symmetric form b Z b Z b X 1 1 0 f(n) = f(x)dx + 2 (f(a) + f(b)) + (x − [x] − 2 )f (x)dx: n=a a a Proof. Let n 2 fa; a + 1; : : : ; b − 1g. Then Z n+1 Z n+1 x − [x]f 0(x)dx = (x − n)f 0(x)dx n n h in+1 Z n+1 Z n+1 = (x − n)f(x) − f(x)dx = f(n + 1) − f(x)dx: n n n By summing over n we get b Z b X Z b (x − [x])f 0(x)dx = f(n) − f(x)dx; a n=a+1 a which implies at once Lemma 4.1.
    [Show full text]
  • COMPLETE MONOTONICITY for a NEW RATIO of FINITE MANY GAMMA FUNCTIONS Feng Qi
    COMPLETE MONOTONICITY FOR A NEW RATIO OF FINITE MANY GAMMA FUNCTIONS Feng Qi To cite this version: Feng Qi. COMPLETE MONOTONICITY FOR A NEW RATIO OF FINITE MANY GAMMA FUNCTIONS. 2020. hal-02511909 HAL Id: hal-02511909 https://hal.archives-ouvertes.fr/hal-02511909 Preprint submitted on 19 Mar 2020 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. COMPLETE MONOTONICITY FOR A NEW RATIO OF FINITE MANY GAMMA FUNCTIONS FENG QI Dedicated to people facing and fighting COVID-19 Abstract. In the paper, by deriving an inequality involving the generating function of the Bernoulli numbers, the author introduces a new ratio of finite many gamma functions, finds complete monotonicity of the second logarithmic derivative of the ratio, and simply reviews complete monotonicity of several linear combinations of finite many digamma or trigamma functions. Contents 1. Preliminaries and motivations 1 2. A lemma 3 3. Complete monotonicity 4 4. A simple review 5 References 7 1. Preliminaries and motivations Let f(x) be an infinite differentiable function on (0; 1). If (−1)kf (k)(x) ≥ 0 for all k ≥ 0 and x 2 (0; 1), then we call f(x) a completely monotonic function on (0; 1).
    [Show full text]
  • Two Series Expansions for the Logarithm of the Gamma Function Involving Stirling Numbers and Containing Only Rational −1 Coefficients for Certain Arguments Related to Π
    J. Math. Anal. Appl. 442 (2016) 404–434 Contents lists available at ScienceDirect Journal of Mathematical Analysis and Applications www.elsevier.com/locate/jmaa Two series expansions for the logarithm of the gamma function involving Stirling numbers and containing only rational −1 coefficients for certain arguments related to π Iaroslav V. Blagouchine ∗ University of Toulon, France a r t i c l e i n f o a b s t r a c t Article history: In this paper, two new series for the logarithm of the Γ-function are presented and Received 10 September 2015 studied. Their polygamma analogs are also obtained and discussed. These series Available online 21 April 2016 involve the Stirling numbers of the first kind and have the property to contain only Submitted by S. Tikhonov − rational coefficients for certain arguments related to π 1. In particular, for any value of the form ln Γ( 1 n ± απ−1)andΨ ( 1 n ± απ−1), where Ψ stands for the Keywords: 2 k 2 k 1 Gamma function kth polygamma function, α is positive rational greater than 6 π, n is integer and k Polygamma functions is non-negative integer, these series have rational terms only. In the specified zones m −2 Stirling numbers of convergence, derived series converge uniformly at the same rate as (n ln n) , Factorial coefficients where m =1, 2, 3, ..., depending on the order of the polygamma function. Explicit Gregory’s coefficients expansions into the series with rational coefficients are given for the most attracting Cauchy numbers −1 −1 1 −1 −1 1 −1 values, such as ln Γ(π ), ln Γ(2π ), ln Γ( 2 + π ), Ψ(π ), Ψ( 2 + π )and −1 Ψk(π ).
    [Show full text]
  • Introduction to Analytic Number Theory More About the Gamma Function We Collect Some More Facts About Γ(S)
    Math 259: Introduction to Analytic Number Theory More about the Gamma function We collect some more facts about Γ(s) as a function of a complex variable that will figure in our treatment of ζ(s) and L(s, χ). All of these, and most of the Exercises, are standard textbook fare; one basic reference is Ch. XII (pp. 235–264) of [WW 1940]. One reason for not just citing Whittaker & Watson is that some of the results concerning Euler’s integrals B and Γ have close analogues in the Gauss and Jacobi sums associated to Dirichlet characters, and we shall need these analogues before long. The product formula for Γ(s). Recall that Γ(s) has simple poles at s = 0, −1, −2,... and no zeros. We readily concoct a product that has the same behavior: let ∞ 1 Y . s g(s) := es/k 1 + , s k k=1 the product converging uniformly in compact subsets of C − {0, −1, −2,...} because ex/(1 + x) = 1 + O(x2) for small x. Then Γ/g is an entire function with neither poles nor zeros, so it can be written as exp α(s) for some entire function α. We show that α(s) = −γs, where γ = 0.57721566490 ... is Euler’s constant: N X 1 γ := lim − log N + . N→∞ k k=1 That is, we show: Lemma. The Gamma function has the product formulas ∞ N ! e−γs Y . s 1 Y k Γ(s) = e−γsg(s) = es/k 1 + = lim N s . (1) s k s N→∞ s + k k=1 k=1 Proof : For s 6= 0, −1, −2,..., the quotient g(s + 1)/g(s) is the limit as N→∞ of N N ! N s Y 1 + s s X 1 Y k + s e1/k k = exp s + 1 1 + s+1 s + 1 k k + s + 1 k=1 k k=1 k=1 N ! N X 1 = s · · exp − log N + .
    [Show full text]