The Eulerian Functions

Total Page:16

File Type:pdf, Size:1020Kb

The Eulerian Functions Appendix A The Eulerian Functions Here we consider the so-called Eulerian functions, namely the well-known Gamma function and Beta function together with some special functions that turn out to be related to them, such as the Psi function and the incomplete gamma functions. We recall not only the main properties and representations of these functions, but we also briefly consider their applications in the evaluation of certain expressions relevant for the fractional calculus. A.1 The Gamma Function The Gamma function .z/ is the most widely used of all the special functions: it is usually discussed first because it appears in almost every integral or series representation of other advanced mathematical functions. We take as its definition the integral formula Z 1 .z/ WD uz1 eu du ; Re z >0: (A.1.1) 0 This integral representation is the most common for , even if it is valid only in the right half-plane of C. The analytic continuation to the left half-plane can be done in different ways. As will be shown later, the domain of analyticity D of is D D C nf0; 1; 2;:::;g : (A.1.2) Using integration by parts, (A.1.1) shows that, at least for Re z >0, satisfies the simple difference equation .z C 1/ D z .z/; (A.1.3) © Springer-Verlag Berlin Heidelberg 2014 269 R. Gorenflo et al., Mittag-Leffler Functions, Related Topics and Applications, Springer Monographs in Mathematics, DOI 10.1007/978-3-662-43930-2 270 A The Eulerian Functions which can be iterated to yield .z C n/ D z .z C 1/::: .z C n 1/ .z/; n2 N : (A.1.4) The recurrence formulas (A.1.3 and A.1.4) can be extended to any z 2 D .In particular, since .1/ D 1,wegetfornon-negative integer values .n C 1/ D nŠ ; n D 0; 1; 2; : : : : (A.1.5) As a consequence can be used to define the Complex Factorial Function zŠ WD .z C 1/ : (A.1.6) By the substitution u D v2 in (A.1.1)wegettheGaussian Integral Representation Z 1 2 .z/ D 2 ev v2z1dv ; Re .z/>0; (A.1.7) 0 which can be used to obtain when z assumes positive semi-integer values. Starting from  à Z C1 1 2 p D ev dv D 1:77245 ; (A.1.8) 2 1 we obtain for n 2 N,  à Z  à C1 1 2 1 .2n 1/ŠŠ p .2n/Š n C D v vn vD D :1 e d n 2n (A.1.9) 2 1 2 2 2 nŠ A.1.1 Analytic Continuation A common way to derive the domain of analyticity (A.1.2) is to carry out the analytic continuation by the mixed representation due to Mittag-Leffler: Z X1 .1/n 1 .z/ D C eu uz1 du ; z 2 D : (A.1.10) nŠ.z C n/ nD0 1 This representation can be obtained from the so-called Prym’s decomposition, namely by splitting the integral in (A.1.1) into two integrals, one over the interval 1The double factorial mŠŠ means mŠŠ D 1 3 :::mif m is odd, and mŠŠ D 2 4 :::mif m is even. A.1 The Gamma Function 271 0 Ä u Ä 1 which is then developed as a series, the other over the interval 1 Ä u Ä 1, which, being uniformly convergent inside C, provides an entire function. The terms of the series (uniformly convergent inside D ) provide the principal parts of at the corresponding poles zn Dn. So we recognize that is analytic in the entire complex plane except at the points zn Dn (n D 0; 1; : : : ), which turn n out to be simple poles with residues Rn D .1/ =nŠ. The point at infinity, being an accumulation point of poles, is an essential non-isolated singularity. Thus is a transcendental meromorphic function. A formal way to obtain the domain of analyticity D is to carry out the required analytical continuation via the Recurrence Formula .z C n/ .z/ D ; (A.1.11) .z C n 1/ .z C n 2/:::.z C 1/ z that is obtained by iterating (A.1.3) written as .z/ D .z C 1/=z.Inthiswaywe can enter the left half-plane step by step. The numerator on the R.H.S. of (A.1.11) is analytic for Re z > n; hence, the L.H.S. is analytic for Re z > n except for simple poles at z D 0; 1;:::;.n C 2/; .n C 1/.Sincen can be arbitrarily large, we deduce the properties discussed above. Another way to interpret the analytic continuation of the Gamma function is provided by the Cauchy–Saalschütz representation, which is obtained by iterated integration by parts in the basic representation (A.1.1). If n 0 denotes any non- negative integer, we have Z Ä 1 1 1 .z/D uz1 eu 1 C u u2 CC.1/nC1 un du (A.1.12) 0 2 n in the strip .n C 1/ < Re z < n. To prove this representation the starting point is provided by the integral Z 1 uz1 Œeu 1 du ; 1<Rez <0: 0 Integration by parts gives (the integrated terms vanish at both limits) Z Z 1 1 1 1 uz1 Œeu 1 du D uzeu du D .z C 1/ D .z/: 0 z 0 z So, by iteration, we get (A.1.12). A.1.2 The Graph of the Gamma Function on the Real Axis Plots of .x/ (continuous line) and 1=.x/ (dashed line) for 4<xÄ 4 are shown in Fig. A.1 and for 0<xÄ 3 in Fig. A.2. 272 A The Eulerian Functions 6 Γ(x) 4 2 1/Γ(x) 0 −2 −4 −6 −4 −3 −2 −1 0 1 2 3 4 x Fig. A.1 Plots of .x/ (continuous line)and1=.x/ (dashed line) 3 2.5 Γ(x) 2 1/Γ(x) 1.5 1 X: 1.462 Y: 0.8856 0.5 0 0 0.5 1 1.5 2 2.5 3 3.5 4 x Fig. A.2 Plots of .x/ (continuous line)and1=.x/ (dashed line) Hereafter we provide some analytical arguments that support the plots on the real axis. In fact, one can get an idea of the graph of the Gamma function on the real axis using the formulas .x/ .x C 1/ D x.x/; .x 1/ D ; x 1 to be iterated starting from the interval 0<xÄ 1,where.x/ !C1as x ! 0C and .1/ D 1. A.1 The Gamma Function 273 For x>0the integral representation (A.1.1) yields .x/ > 0 and 00.x/ > 0 since Z Z 1 1 .x/ D eu ux1 du ;00.x/ D eu ux1 .log u/2 du : 0 0 As a consequence, on the positive real axis .x/ turns out to be positive and convex so that it first decreases and then increases, exhibiting a minimum value. Since .1/ D .2/ D 1, we must have a minimum at some x0, 1<x0 <2. It turns out that x0 D 1:4616 : : : and .x0/ Dp0:8856 : : :; hence x0 is quite close to the point x D 1:5 where attains the value =2 D 0:8862 : : :. On the negative real axis .x/ exhibits vertical asymptotes at x Dn.nD 0; 1; 2; : : :/; it turns out to be positive for 2<x<1, 4<x<3, :::,and negative for 1<x<0, 3<x<2, :::. A.1.3 The Reflection or Complementary Formula .z/.1 z/ D : (A.1.13) sin z This formula, which shows the relationship between the function and the trigonometric sin function, is of great importance together with the recurrence formula (A.1.3). It can be proved in several ways; the simplest proof consists in proving (A.1.13)for0<Rez <1and extending the result by analytic continuation to C except at the points 0; ˙1; ˙2;::: The reflection formula shows that has no zeros. In fact, the zeros cannot be in z D 0; ˙1; ˙2;:::and, if vanished at a non-integer z, then because of (A.1.13), this zero would be a pole of .1 z/, which cannot be true. This fact implies that 1 1= is an entire function. Loosely speaking, .1z/ collects the positive zeros of 1 sin .z/, while .z/ collects the non-positive zeros. A.1.4 The Multiplication Formulas Gauss proved the following Multiplication Formula nY1 k .nz/ D .2/.1n/=2 nnz1=2 .z C /; nD 2;3;::: ; (A.1.14) n kD0 which reduces, for n D 2,toLegendre’s Duplication Formula 1 1 .2z/ D p 22z1=2 .z/.z C /; (A.1.15) 2 2 274 A The Eulerian Functions and, for n D 3,totheTriplication Formula 1 1 2 .3z/ D 33z1=2 .z/.z C /.z C /: (A.1.16) 2 3 3 A.1.5 Pochhammer’s Symbols Pochhammer’s symbols .z/n are defined for any non-negative integer n as .z C n/ .z/ WD z .z C 1/ .z C 2/:::.z C n 1/ D ;n2 N ; (A.1.17) n .z/ with .z/0 D 1. In particular, for z D 1=2 ; we obtain from (A.1.9)  à 1 .n C 1=2/ .2n 1/ŠŠ WD D : n 2 n .1=2/ 2 We extend the above notation to negative integers, defining .z C 1/ .z/ WD z .z 1/ .z 2/:::.z n C 1/ D ;n2 N : (A.1.18) n .z n C 1/ A.1.6 Hankel Integral Representations In 1864 Hankel provided a complex integral representation of the function 1=.z/ valid for unrestricted z;itreads: Z 1 1 et .z/ D t; z 2 C ; z d (A.1.19a) 2i Hat where Ha denotes the Hankel path defined as a contour that begins at t D1ia (a>0), encircles the branch cut that lies along the negative real axis, and ends up at t D1Cib (b>0).
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]
  • Conjugate Trigonometric Integrals! by R
    SOME THEOREMS ON FOURIER TRANSFORMS AND CONJUGATE TRIGONOMETRIC INTEGRALS! BY R. P. BOAS, JR. 1. Introduction. Functions whose Fourier transforms vanish for large val- ues of the argument have been studied by Paley and Wiener.f It is natural, then, to ask what properties are possessed by functions whose Fourier trans- forms vanish over a finite interval. The two classes of functions are evidently closely related, since a function whose Fourier transform vanishes on ( —A, A) is the difference of two suitably chosen functions, one of which has its Fourier transform vanishing outside (—A, A). In this paper, however, we obtain, for a function having a Fourier transform vanishing over a finite interval, a criterion which does not seem to be trivially related to known results. It is stated in terms of a transform which is closely related to the Hilbert trans- form, and which may be written in the simplest (and probably most interest- ing) case, 1 C" /(* + *) - fix - 0 (1) /*(*) = —- sin t dt. J IT J o t2 With this notation, we shall show that a necessary and sufficient condition for a function/(x) belonging to L2(—<», oo) to have a Fourier transform vanishing almost everywhere on ( — 1, 1) is that/**(x) = —/(*) almost every- where;/**^) is the result of applying the transform (1) to/*(x). The result holds also if (sin t)/t is replaced in the integrand in (1) by a certain more general function Hit). As a preliminary result, we discuss the effect of the general transform/*(x) on a trigonometric Stieltjes integral, (2) /(*) = f e^'dait), J -R where a(/) is a normalized! complex-valued function of bounded variation.
    [Show full text]
  • Probabilistic Proofs of Some Beta-Function Identities
    1 2 Journal of Integer Sequences, Vol. 22 (2019), 3 Article 19.6.6 47 6 23 11 Probabilistic Proofs of Some Beta-Function Identities Palaniappan Vellaisamy Department of Mathematics Indian Institute of Technology Bombay Powai, Mumbai-400076 India [email protected] Aklilu Zeleke Lyman Briggs College & Department of Statistics and Probability Michigan State University East Lansing, MI 48825 USA [email protected] Abstract Using a probabilistic approach, we derive some interesting identities involving beta functions. These results generalize certain well-known combinatorial identities involv- ing binomial coefficients and gamma functions. 1 Introduction There are several interesting combinatorial identities involving binomial coefficients, gamma functions, and hypergeometric functions (see, for example, Riordan [9], Bagdasaryan [1], 1 Vellaisamy [15], and the references therein). One of these is the following famous identity that involves the convolution of the central binomial coefficients: n 2k 2n − 2k =4n. (1) k n − k k X=0 In recent years, researchers have provided several proofs of (1). A proof that uses generating functions can be found in Stanley [12]. Combinatorial proofs can also be found, for example, in Sved [13], De Angelis [3], and Miki´c[6]. A related and an interesting identity for the alternating convolution of central binomial coefficients is n n n 2k 2n − 2k 2 n , if n is even; (−1)k = 2 (2) k n − k k=0 (0, if n is odd. X Nagy [7], Spivey [11], and Miki´c[6] discussed the combinatorial proofs of the above identity. Recently, there has been considerable interest in finding simple probabilistic proofs for com- binatorial identities (see Vellaisamy and Zeleke [14] and the references therein).
    [Show full text]
  • List of Reviews by Gianni Pagnini in Mathematical Reviews MR2145045
    List of reviews by Gianni Pagnini in Mathematical Reviews MR2145045 Saxena R.K., Ram J., Suthar D.L., On two-dimensional Saigo–Maeda frac- tional calculus involving two-dimensional H-transforms. Acta Cienc. Indica Math., 30(4), pp. 813–822 (2004) MR2220224 Nishimoto K., N-fractional calculus of a logarithmic function and generali- zed hypergeometric functions. J. Fract. Calc., 29(1), pp. 1–8 (2006) MR2224671 Saxena R.K., Ram J., Chandak S., On two-dimensional generalized Saigo fractional calculus associated with two-dimensional generalized H-transfroms. J. Indian Acad. Math., 27(1), pp. 167–180 (2005) MR2266353 Lin S.-D., Tu S.-T., Srivastava H.M., Some families of multiple infinite sums and associated fractional differintegral formulas for power and composite functions. J. Fract. Calc., 30, pp. 45–58 (2006) MR2286840 Biacino L., Derivatives of fractional orders of continuos functions. Ricerche di Matematica, LIII(2), pp. 231–254 (2004) MR2330471 Chaurasia V.B.L., Srivastava A., A unified approach to fractional calculus pertaining to H-functions. Soochow J. of Math., 33(2), pp. 211–221 (2007) MR2332241 Chaurasia V.B.L., Patni R., Shekhawat A.S., Applications of fractional deri- vatives of certain special functions. Soochow J. of Math., 33(2), pp. 325–334 (2007) MR2355703 Agrawal R., Bansal S.K., A study of unified integral operators involving a general multivariable polynomial and a product of two H¯ −functions. J. Rajasthan Acad. Phy. Sci., 6(3), pp. 289–300 (2007) MR2390179 Benchohra M., Hamani S., Ntouyas S.K., Boundary value problems for dif- ferential equations with fractional order. Surv.
    [Show full text]
  • A Novel Approach to Fractional Calculus -.:: Natural Sciences
    Progr. Fract. Differ. Appl. 4, No. 4, 463-478 (2018) 463 Progress in Fractional Differentiation and Applications An International Journal http://dx.doi.org/10.18576/pfda/040402 A Novel Approach to Fractional Calculus: Utilizing Fractional Integrals and Derivatives of the Dirac Delta Function Evan Camrud1,2 1 Department of Mathematics, Concordia College, Moorhead, MN, USA 2 Department of Mathematics, Iowa State University, Ames, IA, USA Received: 8 Jan. 2018, Revised: 28 Feb. 2018, Accepted: 2 Mar. 2018 Published online: 1 Oct. 2018 Abstract: While the definition of a fractional integral may be codified by Riemann and Liouville, an agreed-upon fractional derivative has eluded discovery for many years. This is likely a result of integral definitions including numerous constants of integration in their results. An elimination of constants of integration opens the door to an operator that reconciles all known fractional derivatives and shows surprising results in areas unobserved before, including the appearance of the Riemann Zeta function and fractional Laplace and Fourier transforms. A new class of functions, known as Zero Functions and closely related to the Dirac delta function, are necessary for one to perform elementary operations of functions without using constants. The operator also allows for a generalization of the Volterra integral equation, and provides a method of solving for Riemann’s complimentary function introduced during his research on fractional derivatives. Keywords: Fractional calculus, fractional differential equations, integral transforms, operations with distributions, special functions. 1 Introduction The concept of derivatives of non-integer order, commonly known as fractional derivatives, first appeared in a letter between L’Hopital and Leibniz in which the question of a half-order derivative was posed [1].
    [Show full text]
  • ME 285A: Foundations of the Theory of Continuous Media
    University of California, Berkeley Department of Mechanical Engineering ME 285A: Foundations of the Theory of Continuous Media Syllabus 1. Introduction. The nonlinear field theories. Mathematical preliminaries. Tensors. Fréchet derivative and differential. 2. Kinematics of Deformable Continua. Bodies. Deformations. Polar decomposition theorem. Derivatives of stretch and rotation tensors. Velocity and vorticity fields. Pull-back and push- forward operations. Lie derivative. Material transport of fields. Internal constraints. Relative deformation gradient. Rivlin-Ericksen tensors. 3. Balance Laws. Conservation of mass. Balance of linear and angular momenta for continua. Stress tensors. Balance of energy. A remark on the second law of thermodynamics. 4. General Constitutive Relations (Purely Mechanical Theory). Relationship between stress and history of deformation. Invariance requirements. Simple materials. Symmetry considerations. Cauchy-elastic and Green-elastic (or hyperelastic) materials. Internally constrained elastic materials. Viscous fluids. Materials with memory. Viscoelastic fluids and solids. Reiner-Rivlin fluids. Rivlin-Ericksen fluids. Materials with fading memory. Linear viscoelasticity. Pseudo-Rigid materials. 5. Thermoelastic Materials. Ideal gas. Entropy and Helmholtz free energy functions. Clausius- Duhem inequality. Thermoelasticity. Approach of Coleman and Noll. Approach of Rivlin: construction of entropy function. Internally constrained thermoelatic materials. Linear thermoelasticity. 6. Thermoelastic Materials with Viscosity. Construction of the entropy function for thermoelastic materials with viscosity, and subject to internal constraints. 7. Thermodynamics of Materials with Memory. Materials with fading memory. Instantaneous elasticity. Open questions. 8. Selected Topics. Mixture theory, nonlocal continua. ______________________________________________________________________________ Grading scheme: Homework (~40%), Midterm (~20%), Final Project (~40%). .
    [Show full text]
  • On the Continuability of Multivalued Analytic Functions to an Analytic Subset*
    Functional Analysis and Its Applications, Vol. 35, No. 1, pp. 51–59, 2001 On the Continuability of Multivalued Analytic Functions to an Analytic Subset* A. G. Khovanskii UDC 517.9 In the paper it is shown that a germ of a many-valued analytic function can be continued analytically along the branching set at least until the topology of this set is changed. This result is needed to construct the many-dimensional topological version of the Galois theory. The proof heavily uses the Whitney stratification. Introduction In the topological version of Galois theory for functions of one variable (see [1–4]) it is proved that the character of location of the Riemann surface of a function over the complex line can prevent the representability of this function by quadratures. This not only explains why many differential equations cannot be solved by quadratures but also gives the sharpest known results on their nonsolvability. I always thought that there is no many-dimensional topological version of Galois theory of full value. The point is that, to construct such a version in the many-variable case, it would be necessary to have information not only on the continuability of germs of functions outside their branching sets but also along these sets, and it seemed that there is nowhere one can obtain this information from. Only in spring of 1999 did I suddenly understand that germs of functions are sometimes automatically continued along the branching set. Therefore, a many-dimensional topological version of the Galois theory does exist. I am going to publish it in forthcoming papers.
    [Show full text]
  • Mathematics 136 – Calculus 2 Summary of Trigonometric Substitutions October 17 and 18, 2016
    Mathematics 136 – Calculus 2 Summary of Trigonometric Substitutions October 17 and 18, 2016 The trigonometric substitution method handles many integrals containing expressions like a2 x2, x2 + a2, x2 a2 p − p p − (possibly including expressions without the square roots!) The basis for this approach is the trigonometric identities 1 = sin2 θ + cos2 θ sec2 θ = tan2 θ +1. ⇒ from which we derive other related identities: a2 (a sin θ)2 = a2(1 sin2 θ) = √a2 cos2 θ = a cos θ p − q − (a tan θ)2 + a2 = a2(tan2 θ +1) = √a2 sec2 θ = a sec θ p q (a sec θ)2 a2 = a2(sec2 θ 1) = a2 tan2 θ = a tan θ p − p − p (Technical note: We usually assume that a > 0 and 0 <θ<π/2 here so that all the trig functions take positive values.) Hence, 1. If our integral contains √a2 x2, the substitution x = a sin θ will convert this radical to the simpler form a cos θ. − 2. If our integral contains √x2 + a2, the substitution x = a tan θ will convert this radical to the simpler form a sec θ. 3. If our integral contains √x2 a2, the substitution x = a sec θ will convert this radical to the simpler form a tan θ. − We substitute for the rest of the integral including the dx. For instance if x = a sin θ, then the dx = a cos θ dθ. If x = a tan θ, then dx = a sec2 θ dθ. If x = a sec θ, then dx = a sec θ tan θ dθ.
    [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]
  • Draftfebruary 16, 2021-- 02:14
    Exactification of Stirling’s Approximation for the Logarithm of the Gamma Function Victor Kowalenko School of Mathematics and Statistics The University of Melbourne Victoria 3010, Australia. February 16, 2021 Abstract Exactification is the process of obtaining exact values of a function from its complete asymptotic expansion. This work studies the complete form of Stirling’s approximation for the logarithm of the gamma function, which consists of standard leading terms plus a remainder term involving an infinite asymptotic series. To obtain values of the function, the divergent remainder must be regularized. Two regularization techniques are introduced: Borel summation and Mellin-Barnes (MB) regularization. The Borel-summed remainder is found to be composed of an infinite convergent sum of exponential integrals and discontinuous logarithmic terms from crossing Stokes sectors and lines, while the MB-regularized remainders possess one MB integral, with similar logarithmic terms. Because MB integrals are valid over overlapping domains of convergence, two MB-regularized asymptotic forms can often be used to evaluate the logarithm of the gamma function. Although the Borel- summed remainder is truncated, albeit at very large values of the sum, it is found that all the remainders when combined with (1) the truncated asymptotic series, (2) the leading terms of Stirling’s approximation and (3) their logarithmic terms yield identical valuesDRAFT that agree with the very high precision results obtained from mathematical software packages.February 16,
    [Show full text]
  • Standard Distributions
    Appendix A Standard Distributions A.1 Standard Univariate Discrete Distributions I. Binomial Distribution B(n, p) has the probability mass function n f(x)= px (x =0, 1,...,n). x The mean is np and variance np(1 − p). II. The Negative Binomial Distribution NB(r, p) arises as the distribution of X ≡{number of failures until the r-th success} in independent trials each with probability of success p. Thus its probability mass function is r + x − 1 r x fr(x)= p (1 − p) (x =0, 1, ··· , ···). x Let Xi denote the number of failures between the (i − 1)-th and i-th successes (i =2, 3,...,r), and let X1 be the number of failures before the first success. Then X1,X2,...,Xr and r independent random variables each having the distribution NB(1,p) with probability mass function x f1(x)=p(1 − p) (x =0, 1, 2,...). Also, X = X1 + X2 + ···+ Xr. Hence ∞ ∞ x x−1 E(X)=rEX1 = r p x(1 − p) = rp(1 − p) x(1 − p) x=0 x=1 ∞ ∞ d d = rp(1 − p) − (1 − p)x = −rp(1 − p) (1 − p)x x=1 dp dp x=1 © Springer-Verlag New York 2016 343 R. Bhattacharya et al., A Course in Mathematical Statistics and Large Sample Theory, Springer Texts in Statistics, DOI 10.1007/978-1-4939-4032-5 344 A Standard Distributions ∞ d d 1 = −rp(1 − p) (1 − p)x = −rp(1 − p) dp dp 1 − (1 − p) x=0 =p r(1 − p) = . (A.1) p Also, one may calculate var(X) using (Exercise A.1) var(X)=rvar(X1).
    [Show full text]
  • 3 Elementary Functions
    3 Elementary Functions We already know a great deal about polynomials and rational functions: these are analytic on their entire domains. We have thought a little about the square-root function and seen some difficulties. The remaining elementary functions are the exponential, logarithmic and trigonometric functions. 3.1 The Exponential and Logarithmic Functions (§30–32, 34) We have already defined the exponential function exp : C ! C : z 7! ez using Euler’s formula ez := ex cos y + iex sin y (∗) and seen that its real and imaginary parts satisfy the Cauchy–Riemann equations on C, whence exp C d z = z is entire (analytic on ). Indeed recall that dz e e . We have also seen several of the basic properties of the exponential function, we state these and several others for reference. Lemma 3.1. Throughout let z, w 2 C. 1. ez 6= 0. ez 2. ez+w = ezew and ez−w = ew 3. For all n 2 Z, (ez)n = enz. 4. ez is periodic with period 2pi. Indeed more is true: ez = ew () z − w = 2pin for some n 2 Z Proof. Part 1 follows trivially from (∗). To prove 2, recall the multiple-angle formulae for cosine and sine. Part 3 requires an induction using part 2 with z = w. Part 4 is more interesting: certainly ew+2pin = ew by the periodicity of sine and cosine. Now suppose ez = ew where z = x + iy and w = u + iv. Then, by considering the modulus and argument, ( ex = eu exeiy = eueiv =) y = v + 2pin for some n 2 Z We conclude that x = u and so z − w = i(y − v) = 2pin.
    [Show full text]