1 Introduction

Total Page:16

File Type:pdf, Size:1020Kb

1 Introduction A Method to Construct Continued-Fraction Approximations and Its Applications Chao-Ping Chen1;∗, H. M. Srivastava2;3 and Qin Wang4 1School of Mathematics and Informatics, Henan Polytechnic University, Jiaozuo 454000, Henan, People’s Republic of China E-Mail: [email protected] ∗Corresponding Author 2Department of Mathematics and Statistics, University of Victoria, Victoria, British Columbia V8W 3R4, Canada 3Department of Medical Research, China Medical University Hospital, China Medical University, Taichung 40402, Taiwan, Republic of China E-Mail: [email protected] 4School of Mathematics and Informatics, Henan Polytechnic University, Jiaozuo 454000, Henan, People’s Republic of China E-Mail: [email protected] Abstract In this paper, we provide a method to construct a continued-fraction approx- imation based upon a given asymptotic expansion. As applications of the method developed here, we establish several continued-fraction approximations for the gam- ma and the digamma (or psi) functions. Finally, some closely-related open problems are also presented. 2010 Mathematics Subject Classification. Primary 11A55, 33B15; Secondary 11Y60. Key Words and Phrases. Gamma function; Digamma (or psi) function; Euler-Mascheroni constant; Asymptotic expansions; Continued-fraction approximations. 1 Introduction The gamma function Γ(x) given by Z 1 Γ(x) = tx−1 e−t dt (x > 0) 0 A Method to Construct Continued-Fraction Approximations and Its Applications 2 is one of the most important functions in mathematical analysis and has applications in many diverse areas. The logarithmic derivative (x) of the gamma function Γ(x) given by Γ0(x) Z x (x) = or ln Γ(x) = (x) dt Γ(x) 1 is known as the psi (or digamma) function. The psi function (x) is connected to the Euler- Mascheroni constant γ and the harmonic numbers Hn by means of the following well-known relation (see [1, p. 258, Eq. (6.3.2)]): (n + 1) = −γ + Hn (n 2 N := f1; 2; 3; · · · g); (1.1) where n X 1 H := (n 2 ) n k N k=1 is the nth harmonic number and γ is the Euler-Mascheroni constant defined by γ = lim Dn = 0:577215664 ··· ; n!1 where n X 1 D = − ln n: (1.2) n k k=1 Various approximations of the psi function (x) are used in the relation (1.1) and interpreted as approximation for the harmonic number Hn or as approximations of the Euler-Mascheroni constant γ. There has been significant interest and research on γ as evidenced by survey papers (see, for details, [14]) and expository books (see, for example, [19]), which reveal its essential properties and surprising connections with other areas of the mathematical sciences. The following two-sided inequality for the difference Dn − γ was established in [28, 33]: 1 1 < D − γ < (n 2 ): 2(n + 1) n 2n N The convergence of the sequence Dn to γ is very slow. By changing the logarithmic term in (1.2), DeTemple [15, 16] presented the following inequality: 1 1 < R − γ < (n 2 ); (1.3) 24(n + 1)2 n 24n2 N where 1 R = H − ln n + : (1.4) n n 2 A Method to Construct Continued-Fraction Approximations and Its Applications 3 On the other hand, Negoi [26] proved that the sequence fTngn2N given by n X 1 1 1 T = − ln n + + (1.5) n k 2 24n k=1 is strictly increasing and convergent to γ. Moreover, Negoi [26] proved that 1 1 < γ − T < (n 2 ): (1.6) 48(n + 1)3 n 48n3 N The following faster approximation formulas can be found in [10, 12]: 1 1 1 X := H − ln n + + − = γ + O n−4 (n ! 1) (1.7) n n 2 24n 48n2 and 1 1 1 23 Y := H − ln n + + − + = γ + O n−5 (n ! 1): (1.8) n n 2 24n 48n2 5760n3 In view of (1.3), (1.6), (1.7) and (1.8), Chen and Mortici [10] posed the following open prob- lem: Open Problem 1. For a given m 2 N; find the constants pj (j = 1; 2; 3; ··· ; m) such that m ! X pj H − ln n 1 + (1.9) n nj j=1 is the fastest sequence which would converge to γ. Yang [31] first presented the solution of Open Problem 1 by using the Bell polynomials of a logarithmic type. Subsequently, other proofs of Open Problem 1 (1.9) were published by Gavrea and Ivan [17, 18], Lin [20], Chen et al. [8], and Chen and Choi [6]. The following familiar Stirling’s formula: p nn n! ∼ 2πn (n ! 1) (1.10) e has many applications in statistical physics, probability theory and number theory. Actually, it was first discovered in 1733 by the French mathematician, Abraham de Moivre (1667–1754), in the form given by p nn n! ∼ constant · n (n ! 1) e A Method to Construct Continued-Fraction Approximations and Its Applications 4 when he was studying the Gaussian distribution and the central limit theorem. Afterwards,p the Scottish mathematician, James Stirling (1692–1770), found the missing constant 2π when he was attempting to give the normal approximation of the binomial distribution. Recently, Sandor´ and Debnath [29, Theorem 5] proved the following inequality for all posi- tive integers n = 2: 1 p nn p nn n 2 2πn < Γ(n + 1) < 2πn (1.11) e e n − 1 and proposed the approximation formula given below: 1 p nn n 2 Γ(n + 1) ∼ 2πn (n ! 1): (1.12) e n − 1 Motivated by the above-mentioned work by Sandor´ and Debnath (1.12), Mortici and Srivas- tava [25] introduced the following class of approximations for all real numbers a and b: ! p nn 1 −b Γ(n + 1) ∼ µ (a; b) n ! 1; µ (a; b) := 2πn 1 + : (1.13) n n e n + a We note that Stirling’s formula (1.10) is obtained from the Mortici-Srivastava result (1.13) in its special case when b = 0. Furthermore, the approximation formula (1.12) can be written as follows: 1 Γ(n + 1) ∼ µ −1; − (n ! 1): n 2 Indeed, in the aforecited work, Mortici and Srivastava [25] proved that 1 1 Γ(n + 1) ∼ µ − ; − (n ! 1) n 2 12 is the best approximation among all approximations given by (1.13). We choose to write this best approximation as follows: ! 1 p nn 1 12 Γ(n + 1) ∼ 2πn 1 + 1 (n ! 1): (1.14) e n − 2 Mortici and Srivastava [25] also developed the approximation formula (1.14) in order to produce a complete asymptotic expansion (see [25, Theorem 2]). In this paper, we provide a method to construct a continued-fraction approximation which is based upon a given asymptotic expansion (see Remarks 1 and 2). As applications of our A Method to Construct Continued-Fraction Approximations and Its Applications 5 continued-fraction approximation, we develop the approximation formula (1.14) to produce sev- eral further continued-fraction approximations (see Theorems 3, 4 and 5). We also establish continued-fraction approximation for the psi function (see Theorem 6). Finally, we present the higher-order estimate for the Euler-Mascheroni constant γ (see Theorem 7). The following lemmas will be useful in our present investigation. Lemma 1 (see [8]). Let 1 X −n g(x) ∼ bnx (x ! 1) n=1 be a given asymptotic expansion. Then the composition exp g(x) has asymptotic expansion of the following form: 1 X −n exp g(x) ∼ anx (x ! 1); n=0 where n 1 X a = 1 and a = kb a (n 2 ): (1.15) 0 n n k n−k N k=1 Lemma 2 (see [5, Theorem 9]). Let k = 1 and n = 0 be integers. Then; for all real numbers x > 0: k+1 (k) Sk(2n; x) < (−1) (x) < Sk(2n + 1; x); (1.16) where p " k−1 # (k − 1)! k! X Y 1 S (p; x) = + + B (2i + j) ; k xk 2xk+1 2i x2i+k i=1 j=1 fBngn2N0 (N0 = N [ f0g) are the Bernoulli numbers defined by 1 t X tn = B : et − 1 n n! n=0 It follows from (1.16) that 1 1 1 1 1 1 + + − + − < 0(x) x 2x2 6x3 30x5 42x7 30x9 1 1 1 1 1 1 5 < + + − + − + (x > 0): (1.17) x 2x2 6x3 30x5 42x7 30x9 66x11 By using the recurrence formula: 1 0(x + 1) = 0(x) − ; x2 A Method to Construct Continued-Fraction Approximations and Its Applications 6 we deduce from (1.17) that, for x > 0, 1 1 1 1 1 1 − + − + − < 0(x + 1) x 2x2 6x3 30x5 42x7 30x9 1 1 1 1 1 1 5 < − + − + − + : (1.18) x 2x2 6x3 30x5 42x7 30x9 66x11 The numerical calculations presented in this work were performed by using the Maple soft- ware for symbolic computations. 2 A Method for the Construction of Continued-Fraction Approxi- mations In this section, we present a method to construct a continued-fraction approximation based upon a given asymptotic expansion (see Remark 1 and Remark 2 below). Theorem 1 generalizes an earlier result [13, Lemma 1.1]. Theorem 1. Let a` 6= 0 (` 2 N) and 1 X aj A(x) ∼ (x ! 1) (2.1) xj j=` be a given asymptotic expansion. Define the function B(x) by a A(x) = ` : B(x) a` Then the function B(x) = A(x) has the following asymptotic expansion: 1 X bj B(x) ∼ x` + b x`−1 + b x`−2 + ··· + b x + b + (x ! 1); −(`−1) −(`−2) −1 0 xj j=1 A Method to Construct Continued-Fraction Approximations and Its Applications 7 where 8 a`+1 >b−(`−1) = − > a` > > a + a b > `+2 `+1 −(`−1) >b−(`−2) = − > a` > > > > .
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]