Bounds for the Logarithm of the Euler Gamma Function and Its Derivatives

Total Page:16

File Type:pdf, Size:1020Kb

Bounds for the Logarithm of the Euler Gamma Function and Its Derivatives Bounds for the logarithm of the Euler gamma function and its derivatives Harold G. Diamond∗ and Armin Strauby Department of Mathematics University of Illinois at Urbana-Champaign 1409 W. Green St, Urbana, IL 61801, United States August 12, 2015 Abstract We consider differences between log Γ(x) and truncations of certain clas- sical asymptotic expansions in inverse powers of x − λ whose coefficients are expressed in terms of Bernoulli polynomials Bn(λ), and we obtain conditions under which these differences are strictly completely monotonic. In the sym- metric cases λ = 0 and λ = 1=2, we recover results of Sonin, N¨orlundand Alzer. Also we show how to derive these asymptotic expansions using the functional equation of the logarithmic derivative of the Euler gamma function, the representation of 1=x as a difference F (x + 1) − F (x), and a backward induction. 1 Introduction The Euler gamma function Γ(x) is widely regarded as the most important special function. Good accounts and many formulas for Γ(x) can be found in [4], [5, x6.3], 2010 Mathematics subject classification. Primary 33B15, 26D07; Secondary 11B68 Key words and phrases. Gamma function, psi function, inequalities, asymptotic expansions, Bernoulli polynomials, complete monotonicity. ∗[email protected] [email protected] 1 [6, Chapter 1], and [12, Chapter 12]. Various bounds are known for Γ(x) and (x), the logarithmic derivative of Γ(x), also called the digamma function, e.g. [1], [2], [5, 6.3.21], [7], [9, Appendix C], [10]. The logarithm of the gamma function has the classical asymptotic expansion p 1 log Γ(x) ∼ log 2π + x − 2 log(x − λ) − (x − λ)(1) 1 X Bn(λ) + (−1)n (x − λ)1−n; n(n − 1) n=2 which we consider for real x > λ and λ 2 [0; 1]. As usual, Bn(λ) denotes the nth Bernoulli function. Our focus here is on differences (and their derivatives) between log Γ(x) and the truncations p 1 LN (λ; x) := log 2π + x − 2 log(x − λ) − (x − λ) N X Bn(λ) + (−1)n (x − λ)1−n n(n − 1) n=2 of this asymptotic series. In the most basic case, we are interested in conditions under which these truncations provide strict upper or lower bounds. These inequalities together with the functional equations provide an effective method for calculating log Γ(x) and its derivatives with arbitrary accuracy. Example 1.1. A pleasant property of (1) is that we can differentiate both sides to obtain asymptotic expansions for the derivatives of log Γ(x), that is, the derivatives of the truncations LN (λ; x) provide asymptotic expansions for the derivatives of log Γ(x). For instance, in the particularly symmetric case λ = 1=2, we have the asymptotic expansion 1 X B2n(1=2) (2) (x) ∼ log(x − 1=2) − (x − 1=2)−2n: 2n n=1 Our main theorem implies the following result of N¨orlund[11, x56, p. 106]: if the series in (2) is truncated at n = N with N ≥ 0 and even, the resulting sum is a lower bound for (x) valid for all x > 1=2; if N is odd, it is an upper bound. In other words, successive truncations of (2) provide lower and upper bounds valid for all x > 1=2. The first examples of these bounds are (3) (x) > log(x − 1=2) and (4) (x) < log(x − 1=2) + (x − 1=2)−2=24 : 2 A function f is said to be completely monotonic on an interval or half-line I if, for all x 2 I and all integers n ≥ 0, the derivative f (n)(x) exists and satisfies (−1)nf (n)(x) ≥ 0. If this inequality is strict, that is, if (−1)nf (n)(x) > 0, then f is called strictly completely monotonic on I. For an introduction to completely monotonic functions and their properties we refer to [13, Chapter IV]. Our main result, which is proved in Section4, is the following. As indicated in Example 1.4, this is an extension of a result of Alzer [1, Theorem 8], which concerns the case λ = 0. Theorem 1.2. Let N ≥ 1. Suppose that either N is even and λ 2 [0; λ0] where λ0 is the unique root of BN (λ) in [0; 1=2], or N is odd and λ 2 [λ0; 1=2] where λ0 is the unique root of BN+1(λ) in [0; 1=2]. Then, dN=2e (−1) (log Γ(x) − LN (λ; x)) is strictly completely monotonic on (λ, 1). We note that λ0 < 1=4 and that λ0 ! 1=4 very rapidly as N increases. This is explained in Example 3.3, which also illustrates that Theorem 1.2 would be false without the conditions on N and λ. On the other hand, the restriction to the case λ 2 [0; 1=2] is made for the sake of exposition, and we invite the interested reader to adjust the details for the case λ 2 [1=2; 1]. Example 1.3. The particular case λ = 1=2 and N = 3 shows that p 1 1 1=24 log Γ(x) − log 2π − x − 2 log x − 2 − 1 + 1 x − 2 is strictly completely monotonic on (1=2; 1). This generalizes the inequality (4). Example 1.4. In the case λ = 0, we find that the functions (x) − L2M (0; x), explicitly given by M p X B2n (x) − log 2π − x − 1 log(x) + x − ; 2 2n(2n − 1)x2n−1 n=1 are strictly completely monotonic on (0; 1) if M is even. Likewise, the functions L2M (0; x) − (x) are strictly completely monotonic on (0; 1) if M is odd. This is a result of Alzer [1, Theorem 8]. This paper is organized as follows. In Section2, we give a first illustration of our approach by proving the claims made in Example 1.1. These bounds for the digamma function are then generalized in Section3, before we prove the main result in Section4. A very brief indication of how to use these bounds for numerical calculations is then given in Section5. 3 2 Proving bounds for the digamma function In this section, we illustrate our approach by proving the claims made in Example 1.1. In particular, we show, by a method we believe new, how to derive asymptotic expansions for (x), such as (2), whose successive truncations provide lower and upper bounds. Assuming there is an approximation of the form (x) ≈ F (x), then the functional equation (x + 1) = (x) + 1=x would give (5) F (x + 1) − F (x) ≈ 1=x: Our method is to find functions F giving upper (resp. lower) bounds in (5) and then show by a \backward induction" that these functions provide bounds for (x). We use one additional fact about the digamma function (x), namely the asymptotic formula (6) (x) = log x + o(1); x ! 1: For completeness, we prove this at the end of the article. Example 2.1. To illustrate our method, we first obtain the lower bound (3). By convexity, Z x+1=2 du 1 log(x + 1=2) − log(x − 1=2) = > : x−1=2 u x This inequality and the functional equation give (x) > (x + 1) − (log(x + 1=2) − log(x − 1=2)) for all x > 1=2, or, equivalently, (7) (x) − log(x − 1=2) > (x + 1) − log(x + 1=2): Now we prove (3) by an induction: repeating (7) with x replaced by x + 1; x + 2;::: and using (6), we see that (x) − log(x − 1=2) > lim inf( (y) − log(y − 1=2)) = 0: y!1 The key to our treatment of the general case is to develop a suitable expression for 1=x with the aid of the Euler{Maclaurin sum formula. The following version of Euler{Maclaurin, given in [9, Appendix B] is convenient for our application. 4 Theorem 2.2 (Euler{Maclaurin). Let K be a positive integer, and suppose that f is K times continuously differentiable on [a; b]. Let Bk(·) denote the k-th Bernoulli function. Then, X Z b f(n) = f(x)dx a<n≤b a K X (−1)k + (B (fbg)f (k−1)(b) − B (fag)f (k−1)(a)) k! k k k=1 K Z b (−1) (K) − BK (fxg)f (x)dx; K! a where ftg := t − btc. We now show the following result of N¨orlund[11, x56, p. 106] stated in Exam- ple 1.1. Theorem 2.3. The asymptotic expansion (2) for (x) holds as x ! 1. Moreover, if the series in (2) is truncated at n = N with N ≥ 0 and even, the resulting sum is a lower bound for (x) valid for all x > 1=2; if N is odd, it is an upper bound. Proof. Take f(t) := 1=(x+t), a = −1=2, b = 1=2 in Theorem 2.2. Because Bk(1=2) = 0 for all odd k, we restrict the series to even indices. Starting with the one-term (!) sum 1 X 1 = ; x x + n −1=2<n≤1=2 we have Z 1=2 N 1 dt X B2n(1=2) 1 1 = − − + R ; x x + t 2n (x + 1=2)2n (x − 1=2)2n N −1=2 n=1 where Z 1=2 dt RN := − B2N (ftg) 2N+1 : −1=2 (x + t) For the claimed inequalities, we establish the sign of RN with the aid of several properties of Bernoulli functions, all of which are given in [9]. For −1=2 < t < 0, B2N (ftg) = B2N (1 + t) = B2N (−t): 5 Thus, we have Z 1=2 1 1 RN = − B2N (t) 2N+1 + 2N+1 dt : 0 (x − t) (x + t) Next, Z y B2N+1(y) = (2N + 1) B2N (t) dt: 0 Integrating by parts and noting that the integrated term vanishes because B2N+1(1=2) = B2N+1(0) = 0, we find Z 1=2 1 1 RN = B2N+1(t) 2N+2 − 2N+2 dt : 0 (x − t) (x + t) N+1 For 0 < t < 1=2; sgn(B2N+1(t)) = (−1) .
Recommended publications
  • 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]
  • 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]
  • 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]
  • Euler-Maclaurin and Euler-Boole Formulas
    Appendix A Euler-MacLaurin and Euler-Boole Formulas A.1 A Taylor Formula The classical Taylor formula Z Xm xk x .x t/m f .x/ D @kf .0/ C @mC1f .t/dt kŠ 0 mŠ kD0 xk can be generalized if we replace the polynomial kŠ by other polynomials (Viskov 1988; Bourbaki 1959). Definition If is a linear form on C0.R/ such that .1/ D 1,wedefinethe polynomials .Pn/ by: P0 D 1 @Pn D Pn1 , .Pn/ D 0 for n 1 P . / k The Generating Function k0 Pk x z We have formally X X X k k k @x. Pk.x/z / D Pk1.x/z D z. Pk.x/z / k0 k1 k0 © Springer International Publishing AG 2017 175 B. Candelpergher, Ramanujan Summation of Divergent Series, Lecture Notes in Mathematics 2185, DOI 10.1007/978-3-319-63630-6 176 A Euler-MacLaurin and Euler-Boole Formulas thus X k xz Pk.x/z D C.z/e k0 To evaluate C.z/ we use the notation x for and by definition of .Pn/ we can write X X k k x. Pk.x/z / D x.Pk.x//z D 1 k0 k0 X k xz xz x. Pk.x/z / D x.C.z/e / D C.z/x.e / k0 1 this gives C.z/ D xz . Thus the generating function of the sequence .Pn/ is x.e / X n xz Pn.x/z D e =M.z/ n xz where the function M is defined by M.z/ D x.e /: Examples P xn n xz .
    [Show full text]
  • On Some Series Representations of the Hurwitz Zeta Function Mark W
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Elsevier - Publisher Connector Journal of Computational and Applied Mathematics 216 (2008) 297–305 www.elsevier.com/locate/cam On some series representations of the Hurwitz zeta function Mark W. Coffey Department of Physics, Colorado School of Mines, Golden, CO 80401, USA Received 21 November 2006; received in revised form 3 May 2007 Abstract A variety of infinite series representations for the Hurwitz zeta function are obtained. Particular cases recover known results, while others are new. Specialization of the series representations apply to the Riemann zeta function, leading to additional results. The method is briefly extended to the Lerch zeta function. Most of the series representations exhibit fast convergence, making them attractive for the computation of special functions and fundamental constants. © 2007 Elsevier B.V. All rights reserved. MSC: 11M06; 11M35; 33B15 Keywords: Hurwitz zeta function; Riemann zeta function; Polygamma function; Lerch zeta function; Series representation; Integral representation; Generalized harmonic numbers 1. Introduction (s, a)= ∞ (n+a)−s s> a> The Hurwitz zeta function, defined by n=0 for Re 1 and Re 0, extends to a meromorphic function in the entire complex s-plane. This analytic continuation to C has a simple pole of residue one. This is reflected in the Laurent expansion ∞ n 1 (−1) n (s, a) = + n(a)(s − 1) , (1) s − 1 n! n=0 (a) (a)=−(a) =/ wherein k are designated the Stieltjes constants [3,4,9,13,18,20] and 0 , where is the digamma a a= 1 function.
    [Show full text]
  • A New Entire Factorial Function
    A new entire factorial function Matthew D. Klimek Laboratory for Elementary Particle Physics, Cornell University, Ithaca, NY, 14853, USA Department of Physics, Korea University, Seoul 02841, Republic of Korea Abstract We introduce a new factorial function which agrees with the usual Euler gamma function at both the positive integers and at all half-integers, but which is also entire. We describe the basic features of this function. The canonical extension of the factorial to the complex plane is given by Euler's gamma function Γ(z). It is the only such extension that satisfies the recurrence relation zΓ(z) = Γ(z + 1) and is logarithmically convex for all positive real z [1, 2]. (For an alternative function theoretic characterization of Γ(z), see [3].) As a consequence of the recurrence relation, Γ(z) is a meromorphic function when continued onto − the complex plane. It has poles at all non-positive integer values of z 2 Z0 . However, dropping the two conditions above, there are infinitely many other extensions of the factorial that could be constructed. Perhaps the second best known factorial function after Euler's is that of Hadamard [4] which can be expressed as sin πz z z + 1 H(z) = Γ(z) 1 + − ≡ Γ(z)[1 + Q(z)] (1) 2π 2 2 where (z) is the digamma function d Γ0(z) (z) = log Γ(z) = : (2) dz Γ(z) The second term in brackets is given the name Q(z) for later convenience. We see that the Hadamard gamma is a certain multiplicative modification to the Euler gamma.
    [Show full text]
  • Transcendence of Various Infinite Series Applications of Baker's
    Transcendence of Various Infinite Series and Applications of Baker's Theorem by Chester Jay Weatherby A thesis submitted to the Department of Mathematics and Statistics in conformity with the requirements for the degree of Doctor of Philosophy Queen's University Kingston, Ontario, Canada April 2009 Copyright c Chester Jay Weatherby, 2009 Library and Archives Bibliothèque et Canada Archives Canada Published Heritage Direction du Branch Patrimoine de l’édition 395 Wellington Street 395, rue Wellington Ottawa ON K1A 0N4 Ottawa ON K1A 0N4 Canada Canada Your file Votre référence ISBN: 978-0-494-65340-1 Our file Notre référence ISBN: 978-0-494-65340-1 NOTICE: AVIS: The author has granted a non- L’auteur a accordé une licence non exclusive exclusive license allowing Library and permettant à la Bibliothèque et Archives Archives Canada to reproduce, Canada de reproduire, publier, archiver, publish, archive, preserve, conserve, sauvegarder, conserver, transmettre au public communicate to the public by par télécommunication ou par l’Internet, prêter, telecommunication or on the Internet, distribuer et vendre des thèses partout dans le loan, distribute and sell theses monde, à des fins commerciales ou autres, sur worldwide, for commercial or non- support microforme, papier, électronique et/ou commercial purposes, in microform, autres formats. paper, electronic and/or any other formats. The author retains copyright L’auteur conserve la propriété du droit d’auteur ownership and moral rights in this et des droits moraux qui protège cette thèse. Ni thesis. Neither the thesis nor la thèse ni des extraits substantiels de celle-ci substantial extracts from it may be ne doivent être imprimés ou autrement printed or otherwise reproduced reproduits sans son autorisation.
    [Show full text]
  • POCKETBOOKOF MATHEMATICAL FUNCTIONS Abridged Edition of Handbook of Mathematical Functions Milton Abramowitz and Irene A
    POCKETBOOKOF MATHEMATICAL FUNCTIONS Abridged edition of Handbook of Mathematical Functions Milton Abramowitz and Irene A. Stegun (eds.) Material selected by Michael Danos and Johann Rafelski 1984 Verlag Harri Deutsch - Thun - Frankfurt/Main CONTENTS Forewordtothe Original NBS Handbook 5 Pref ace 6 2. PHYSICAL CONSTANTS AND CONVERSION FACTORS 17 A.G. McNish, revised by the editors Table 2.1. Common Units and Conversion Factors 17 Table 2.2. Names and Conversion Factors for Electric and Magnetic Units 17 Table 2.3. AdjustedValuesof Constants 18 Table 2.4. Miscellaneous Conversion Factors 19 Table 2.5. FactorsforConvertingCustomaryU.S. UnitstoSIUnits 19 Table 2.6. Geodetic Constants 19 Table 2.7. Physical andNumericalConstants 20 Table 2.8. Periodic Table of the Elements 21 Table 2.9. Electromagnetic Relations 22 Table 2.10. Radioactivity and Radiation Protection 22 3. ELEMENTARY ANALYTICAL METHODS 23 Milton Abramowitz 3.1. Binomial Theorem and Binomial Coefficients; Arithmetic and Geometrie Progressions; Arithmetic, Geometrie, Harmonie and Generalized Means 23 3.2. Inequalities 23 3.3. Rules for Differentiation and Integration 24 3.4. Limits, Maxima and Minima 26 3.5. Absolute and Relative Errors 27 3.6. Infinite Series 27 3.7. Complex Numbers and Functions 29 3.8. Algebraic Equations 30 3.9. Successive Approximation Methods 31 3.10. TheoremsonContinuedFractions 32 4. ELEMENTARY TRANSCENDENTAL FUNCTIONS 33 Logarithmic, Exponential, Circular and Hyperbolic Functions Ruth Zucker 4.1. Logarithmic Function 33 4.2. Exponential Function 35 4.3. Circular Functions 37 4.4. Inverse Circular Functions 45 4.5. Hyperbolic Functions 49 4.6. InverseHyberbolicFunctions 52 5. EXPONENTIAL INTEGRAL AND RELATED FUNCTIONS .... 56 Walter Gautschi and William F.
    [Show full text]
  • A Note on the Zeros and Local Extrema of Digamma Related Functions
    A note on the zeros and local extrema of Digamma related functions Istv´an Mez˝o 1 School of Mathematics and Statistics, Nanjing University of Information Science and Technology, No.219 Ningliu Rd, Pukou, Nanjing, Jiangsu, P. R. China Abstract Little is known about the zeros of the Digamma function. Establishing some Weier- strassian infinite product representations for a given regularization of the Digamma function we find interesting sums of its zeros. In addition, we study the same ques- tions for the zeros of the logarithmic derivative of the Barnes G function. At the end of the paper we provide rather accurate approximations of the hyperfactorial where, rather interestingly, the Lambert function appears. Key words: Euler Gamma function; Digamma function; Barnes G function; zeros; Weierstrass Product Theorem 1991 MSC: 33B15, 30C15, 30D99 arXiv:1409.2971v2 [math.CV] 9 Feb 2016 1 Introduction The Euler gamma function is defined by the improper integral ∞ Γ(z)= e−ztz−1dt (ℜ(z) > 0). Z0 The logarithmic derivative of Γ is the Digamma function Γ′(z) ψ(z)= (z ∈ C \{0, −1, −2,... }). Γ(z) Email address: [email protected] (Istv´an Mez˝o). 1 The research of Istv´an Mez˝owas supported by the Scientific Research Foundation of Nanjing University of Information Science & Technology, the Startup Foundation for Introducing Talent of NUIST. Project no.: S8113062001, and the National Natural Science Foundation for China. Grant no. 11501299. Preprint submitted to 10 February 2016 This function is analytic everywhere except the non positive integers, and it has first order poles in these points. It is known that the Digamma function has only real and simple zeros, and all of them except only one are negative.
    [Show full text]
  • Clustering, Features Selection, and Outliers Rejection
    Statistical Modeling for Simultaneous Data Clustering, Features Selection, and Outliers Rejection Khaled Almakadmeh A Thesis in The Concordia Institute for Information Systems Engineering Presented in Partial Fulfillment of the Requirements for the Degree of Master of Applied Science in Information Systems Security at Concordia University Montréal, Québec, Canada May 2010 © Khaled Almakadmeh, 2010 Library and Archives Bibliothèque et ?F? Canada Archives Canada Published Heritage Direction du Branch Patrimoine de l'édition 395 Wellington Street 395, rue Wellington Ottawa ON K1A 0N4 OttawaONK1A0N4 Canada Canada Your We Votre référence ISBN: 978-0-494-71008-1 Our file Notre référence ISBN: 978-0-494-71008-1 NOTICE: AVIS: The author has granted a non- L'auteur a accordé une licence non exclusive exclusive license allowing Library and permettant à la Bibliothèque et Archives Archives Canada to reproduce, Canada de reproduire, publier, archiver, publish, archive, preserve, conserve, sauvegarder, conserver, transmettre au public communicate to the public by par télécommunication ou par l'Internet, prêter, telecommunication or on the Internet, distribuer et vendre des thèses partout dans le loan, distribute and sell theses monde, à des fins commerciales ou autres, sur worldwide, for commercial or non- support microforme, papier, électronique et/ou commercial purposes, in microform, autres formats. paper, electronic and/or any other formats. The author retains copyright L'auteur conserve la propriété du droit d'auteur ownership and moral rights in this et des droits moraux qui protège cette thèse. Ni thesis. Neither the thesis nor la thèse ni des extraits substantiels de celle-ci substantial extracts from it may be ne doivent être imprimés ou autrement printed or otherwise reproduced reproduits sans son autorisation.
    [Show full text]
  • Some Inequalities for the Trigamma Function in Terms of the Digamma
    SOME INEQUALITIES FOR THE TRIGAMMA FUNCTION IN TERMS OF THE DIGAMMA FUNCTION FENG QI AND CRISTINEL MORTICI Abstract. In the paper, the authors establish three kinds of double inequal- ities for the trigamma function in terms of the exponential function to powers of the digamma function. These newly established inequalities extend some known results. The method in the paper utilizes some facts from the asymp- totic theory and is a natural way to solve problems for approximating some quantities for large values of the variable. 1. Motivations and main results The classical Euler gamma function may be defined for (z) > 0 by ∞ ℜ − − Γ(z)= tz 1e t d t. Z0 The logarithmic derivative of the gamma function Γ(z) is denoted by d Γ′(z) ψ(z)= [ln Γ(z)] = d z Γ(z) and called the digamma function. The derivatives ψ′(z) and ψ′′(z) are called the trigamma and tetragamma functions respectively. As a whole, the functions ψ(k)(z) for k 0 N are called the polygamma functions. These functions are widely used in theoretical∈{ }∪ and practical problems in all branches of mathematical science. Con- sequently, many mathematicians were preoccupied to establish new results about the gamma function, polygamma functions, and other related functions. 1.1. The first main result. In 2007, Alzer and Batir [3, Corollary] discovered that the double inequality 1 1 √2π xx exp x ψ(x + α) < Γ(x) < √2π xx exp x ψ(x + β) (1.1) arXiv:1503.03020v1 [math.CA] 7 Mar 2015 − − 2 − − 2 1 holds for x> 0 if and only if α 3 and β 0.
    [Show full text]
  • Arxiv:1705.06547V1 [Math.CA] 18 May 2017 Hysatisfy They Am Ucin Oyam Ucin,Inequalities
    INEQUALITIES FOR THE INVERSES OF THE POLYGAMMA FUNCTIONS NECDET BATIR Abstract. We provide an elementary proof of the left side in- equality and improve the right inequality in 1 n+1 n! − − < ((−1)n 1ψ(n)) 1(x) x − (x−1/n + α)−n 1 n! n+1 < , x − (x−1/n + β)−n where α = [(n − 1)!]−1/n and β = [n!ζ(n + 1)]−1/n, which was proved in [6], and we prove the following inequalities for the inverse of the digamma function ψ. 1 − 1 < ψ 1(x) <ex + , x ∈ R. log(1 + e−x) 2 The proofs are based on nice applications of the mean value theo- rem for differentiation and elementary properties of the polygamma functions. 1. introduction As it is well known for a positive real number x the gamma function Γ is defined to be ∞ Γ(x)= ux−1e−udu. Z0 It is a common knowledge that the gamma function plays a special role arXiv:1705.06547v1 [math.CA] 18 May 2017 in the theory of special functions. The most important function related to the gamma function is the digamma or psi function ψ(x), which is defined by logarithmic derivative of the gamma function Γ(x), that is, Γ′(x) ψ(x)= . Γ(x) The digamma function ψ is closely related with the Euler-Mascheroni 1 ··· 1 constant γ (= 0.57721...) and harmonic numbers Hn =1+ 2 + + n . They satisfy ψ(n +1) = −γ + Hn (n ∈ N). In [14] it is proved that 2000 Mathematics Subject Classification. Primary: 33B15; Secondary: 26D07.
    [Show full text]