DOCSLIB.ORG

Arxiv:1803.07841V2 [Math.CA] 5 Jul 2018 Singularity at This Point

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Arxiv:1803.07841V2 [Math.CA] 5 Jul 2018 Singularity at This Point ASYMPTOTIC EXPANSIONS FOR THE INCOMPLETE GAMMA FUNCTION IN THE TRANSITION REGIONS GERGO} NEMES AND ADRI B. OLDE DAALHUIS Abstract. We construct asymptotic expansions for the normalised incom- plete gamma function Q(a; z) = Γ(a; z)=Γ(a) that are valid in the transition regions, including the case z ≈ a, and have simple polynomial coefficients. For Bessel functions, these type of expansions are well known, but for the nor- malised incomplete gamma function they were missing from the literature. A detailed historical overview is included. We also derive an asymptotic expan- sion for the corresponding inverse problem, which has importance in proba- bility theory and mathematical statistics. The coefficients in this expansion are again simple polynomials, and therefore its implementation is straightfor- ward. As a byproduct, we give the first complete asymptotic expansion as a ! −∞ of the unique negative zero of the regularised incomplete gamma function γ∗(a; x). 1. Introduction and main results The normalised incomplete gamma function Q(a; z) = Γ(a; z)=Γ(a) is one of the most widely used special functions of two variables. It is used in constructing gamma distributions [9, Ch. 17], which appear naturally in the theory associated with normally distributed random variables. In many applications, Poisson random variables are used in which the Poisson rate is not fixed, and the function Q(a; z) plays an important role in these cases (see, for example, [6]). Consequently, efficient and accurate approximations for the incomplete gamma function are essential. In many applications the inverse problem also arises naturally, that is, the problem of determining the solution x, real and positive, of the equation Q(a; x) = q for given 0 < q < 1 and a > 0. In the important papers [27, 31], Temme established asymptotic expansions for the normalised incomplete gamma function as a ! 1, that are valid uniformly for jzj 2 [0; 1). The asymptotic expansions of the corresponding inverse functions are discussed in [29]. The large region of validity is important in applications, and many papers in physics, applied statistics, network and control theory, and so on, refer to these uniform asymptotic expansions (see, e.g., [11, 14, 21, 22, 26, 33]). However, applying Temme's expansions in the neighbourhood of the (most in- teresting) point z = a can be difficult because its coefficients possess a removable arXiv:1803.07841v2 [math.CA] 5 Jul 2018 singularity at this point. To overcome this difficulty, Temme [27] gave power series expansions for these coefficients. Having to use local power series expansions for the coefficients in uniform asymptotic expansions, is not a very elegant solution. In this paper, we propose an alternative approach to the problem. We will construct what we call transition region expansions, and provide full details of the inversion of these new expansions. These are expansions that are valid in the regions in which Q(a; z) changes dramatically, and their coefficients are 2010 Mathematics Subject Classification. 33B20, 41A60. Key words and phrases. incomplete gamma functions, asymptotic expansions, transition regions. The authors' research was supported by a research grant (GRANT11863412/70NANB15H221) from the National Institute of Standards and Technology. 1 2 G. NEMES AND A. B. OLDE DAALHUIS polynomials satisfying simple recurrence relations. The region of validity overlaps with those of the non-uniform \outer" expansions. Furthermore, the coefficients of their inversions are simple polynomials, whose computation and implementation are straightforward. It is surprising to us that these transition region expansions for the normalised incomplete gamma function have not yet been discussed in the literature, given the fact that expansions of similar type for Bessel functions are well known (see [18, x10.19(iii)]). What our new expansions and the transition region expansions for the Bessel functions have in common is that both mimic the corresponding uniform expansions; compare the similarities between (1.1) and (2.6), and between [18, Eq. 10.19.8] and [18, Eq. 10.20.4]. The ideas in this paper should also be applicable to other cumulative distribution functions using the results of [28]. The proofs of our new results rely heavily on the uniform asymptotic expansions by Temme. Thus, we shall provide full details of his expansions, and even include new recurrence relations for the coefficients in the local power series expansions mentioned above (see AppendixB). For more information about existing results in the literature and their connection with those we prove in this paper, the reader is referred to the historical overview in Section2. The following theorem is the main result of the paper. Theorem 1.1. The normalised incomplete gamma function admits the asymptotic expansions 2 1 1 1 − 1 1 τ X Cn(τ) (1.1) Q a; a + τa 2 ∼ erfc 2 2 τ + p exp − 2 2 an=2 2πa n=0 and ±πia e 1 ±πi 1 − 1 Q −a; (a + τa 2 )e ∼ ± erfc ±2 2 iτ 2i sin(πa) 2 (1.2) 2 1 n i τ X (−i) Cn(iτ) − p exp ; 2 an=2 2πa n=0 as a ! 1 in the sector j arg aj ≤ π − δ < π, uniformly with respect to bounded complex values of τ. Here, erfc denotes the complementary error function [18, Eq. 7.2.2]. The coefficients Cn(τ) are polynomials in τ of degree 3n + 2 and satisfy 1 1 C (τ) = τ 2 − ; 0 3 3 0 00 2 2 0 00 Cn(τ) + τCn(τ) − Cn(τ) = τ(τ − 2)Cn−1(τ) − (2τ − 1)Cn−1(τ) + τCn−1(τ) for n ≥ 1. In addition, the even- and odd-order polynomials are even and odd functions, respectively. The explicit forms of the next two coefficients in (1.1) are as follows: 1 11 1 1 29 133 23 1 C (τ) = τ 5 − τ 3 + τ; C (τ) = τ 8 − τ 6 + τ 4 − τ 2 − : 1 18 36 12 2 162 324 540 540 540 For higher coefficients and for an algorithm to generate them, see the end of Section 3. The most common case in applications is when both variables of Q(a; z) are real and positive. It is possible to combine our transition region expansion (1.1), with the outer expansions that exist in the literature and cover the complete range τ > −a1=2 (or equivalently, the range z > 0). First we note that in Theorem 1.1 we assume that τ is bounded. We can weaken this restriction as follows. Since µ Cn(τ) is a polynomial of degree 3n + 2 it follows that we can allow τ = O (jaj ) as 1 a ! 1, as long as µ < 6 . See the proposition below. In Section2, we will give ASYMPTOTICS FOR THE INCOMPLETE GAMMA FUNCTION 3 the expansions (2.1) and (2.2). These are also asymptotic expansions with simple coefficients, and with the notation of this section, these two outer expansions are valid for jτj ≥ aµ as a ! +1, with any µ > 0 (see AppendixA). Therefore, by combining these three expansions, we can cover the whole range τ > −a1=2. Proposition 1.2. The expansions (1.1) and (1.2) are still valid in the sense of gen- eralised asymptotic expansions (see, e.g., [18, x2.1(v)]) provided that τ = O (jajµ) 1 as a ! 1 with some constant µ < 6 . In the inverse problem we are trying to find x = x(a; q) such that (1.3) Q(a; x) = q holds, where 0 < q < 1 and a > 0 are given. The inversion of the normalised incom- plete gamma function is an important topic in mathematical statistics and probabil- ity theory, in particular for computing percentage points of the χ2-distribution (see, for example, [8]). A suitable inversion of the expansion (1.1) yields the following result. Theorem 1.3. For any 0 < q < 1, define τ0 to be the unique real root of the equation 1 1 − 2 (1.4) q = 2 erfc 2 τ0 : Then the inverse function x = x(a; q) that satisfies the equation (1.3) has the as- ymptotic expansion 1 1 X dn(τ0) (1.5) x(a; q) ∼ a + τ a 2 + ; 0 an=2 n=0 as a ! +1, provided that q is bounded away from 0 and 1. The coefficients dn(τ0) are polynomials in τ0 of degree n + 2. In addition, the even- and odd-order polynomials are even and odd functions, respectively. The first few coefficients in (1.5) are as follows: 1 1 1 7 1 7 8 d (τ ) = τ 2 − ; d (τ ) = τ 3 − τ ; d (τ ) = − τ 4 − τ 2 + : 0 0 3 0 3 1 0 36 0 36 0 2 0 270 0 810 0 405 For higher coefficients and for an algorithm to generate them, see Section5. In Theorem 1.3, we assume that q is bounded away from 0 and 1. In practice, these bounds can be very small. It follows from (1.4) that τ0 ! +1 as q ! 0 and τ0 ! −∞ as q ! 1. Since 2 1 1 τ 1 − 2 0 2 erfc 2 τ0 ∼ p exp − 2πτ0 2 as τ0 ! +1, and 2 1 1 τ 1 − 2 0 2 erfc 2 τ0 ∼ 1 + p exp − 2πτ0 2 p as τ0 ! −∞, we can infer that jτ0j ∼ −2 log(q(1 − q)) as q ! 0 or 1. Since the dn(τ0)'s are polynomials in τ0 of degree n + 2, for (1.5) to be an asymptotic 1=2 expansion we need τ0 = o jaj as a ! 1. Hence, we need log(q(1−q)) = o(jaj), i.e., q can be sub-exponentially close in jaj to 0 or 1.
Load more

Recommended publications
  • Lecture 5: Complex Logarithm and Trigonometric Functions
    LECTURE 5: COMPLEX LOGARITHM AND TRIGONOMETRIC FUNCTIONS Let C∗ = C \{0}. Recall that exp : C → C∗ is surjective (onto), that is, given w ∈ C∗ with w = ρ(cos φ + i sin φ), ρ = |w|, φ = Arg w we have ez = w where z = ln ρ + iφ (ln stands for the real log) Since exponential is not injective (one one) it does not make sense to talk about the inverse of this function. However, we also know that exp : H → C∗ is bijective. So, what is the inverse of this function? Well, that is the logarithm. We start with a general definition Definition 1. For z ∈ C∗ we define log z = ln |z| + i argz. Here ln |z| stands for the real logarithm of |z|. Since argz = Argz + 2kπ, k ∈ Z it follows that log z is not well defined as a function (it is multivalued), which is something we find difficult to handle. It is time for another definition. Definition 2. For z ∈ C∗ the principal value of the logarithm is defined as Log z = ln |z| + i Argz. Thus the connection between the two definitions is Log z + 2kπ = log z for some k ∈ Z. Also note that Log : C∗ → H is well defined (now it is single valued). Remark: We have the following observations to make, (1) If z 6= 0 then eLog z = eln |z|+i Argz = z (What about Log (ez)?). (2) Suppose x is a positive real number then Log x = ln x + i Argx = ln x (for positive real numbers we do not get anything new).
    [Show full text]
  • Asymptotic Estimates of Stirling Numbers
    Asymptotic Estimates of Stirling Numbers By N. M. Temme New asymptotic estimates are given of the Stirling numbers s~mJ and (S)~ml, of first and second kind, respectively, as n tends to infinity. The approxima­ tions are uniformly valid with respect to the second parameter m. 1. Introduction The Stirling numbers of the first and second kind, denoted by s~ml and (S~~m>, respectively, are defined through the generating functions n x(x-l)···(x-n+l) = [. S~m)Xm, ( 1.1) m= 0 n L @~mlx(x -1) ··· (x - m + 1), ( 1.2) m= 0 where the left-hand side of (1.1) has the value 1 if n = 0; similarly, for the factors in the right-hand side of (1.2) if m = 0. This gives the 'boundary values' Furthermore, it is convenient to agree on s~ml = @~ml= 0 if m > n. The Stirling numbers are integers; apart from the above mentioned zero values, the numbers of the second kind are positive; those of the first kind have the sign of ( - l)n + m. Address for correspondence: Professor N. M. Temme, CWI, P.O. Box 4079, 1009 AB Amsterdam, The Netherlands. e-mail: [email protected]. STUDIES IN APPLIED MATHEMATICS 89:233-243 (1993) 233 Copyright © 1993 by the Massachusetts Institute of Technology Published by Elsevier Science Publishing Co., Inc. 0022-2526 /93 /$6.00 234 N. M. Temme Alternative generating functions are [ln(x+1)r :x n " s(ln)~ ( 1.3) m! £..,, n n!' n=m ( 1.4) The Stirling numbers play an important role in difference calculus, combina­ torics, and probability theory.
    [Show full text]
  • Asymptotics of Singularly Perturbed Volterra Type Integro-Differential Equation
    Konuralp Journal of Mathematics, 8 (2) (2020) 365-369 Konuralp Journal of Mathematics Journal Homepage: www.dergipark.gov.tr/konuralpjournalmath e-ISSN: 2147-625X Asymptotics of Singularly Perturbed Volterra Type Integro-Differential Equation Fatih Say1 1Department of Mathematics, Faculty of Arts and Sciences, Ordu University, Ordu, Turkey Abstract This paper addresses the asymptotic behaviors of a linear Volterra type integro-differential equation. We study a singular Volterra integro equation in the limiting case of a small parameter with proper choices of the unknown functions in the equation. We show the effectiveness of the asymptotic perturbation expansions with an instructive model equation by the methods in superasymptotics. The methods used in this study are also valid to solve some other Volterra type integral equations including linear Volterra integro-differential equations, fractional integro-differential equations, and system of singular Volterra integral equations involving small (or large) parameters. Keywords: Singular perturbation, Volterra integro-differential equations, asymptotic analysis, singularity 2010 Mathematics Subject Classification: 41A60; 45M05; 34E15 1. Introduction A systematic approach to approximation theory can find in the subject of asymptotic analysis which deals with the study of problems in the appropriate limiting regimes. Approximations of some solutions of differential equations, usually containing a small parameter e, are essential in the analysis. The subject has made tremendous growth in recent years and has a vast literature. Developed techniques of asymptotics are successfully applied to problems including the classical long-standing problems in mathematics, physics, fluid me- chanics, astrodynamics, engineering and many diverse fields, for instance, see [1, 2, 3, 4, 5, 6, 7]. There exists a plethora of examples of asymptotics which includes a wide variety of problems in rich results from the very practical to the highly theoretical analysis.
    [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]
  • Complex Analysis
    Complex Analysis Andrew Kobin Fall 2010 Contents Contents Contents 0 Introduction 1 1 The Complex Plane 2 1.1 A Formal View of Complex Numbers . .2 1.2 Properties of Complex Numbers . .4 1.3 Subsets of the Complex Plane . .5 2 Complex-Valued Functions 7 2.1 Functions and Limits . .7 2.2 Infinite Series . 10 2.3 Exponential and Logarithmic Functions . 11 2.4 Trigonometric Functions . 14 3 Calculus in the Complex Plane 16 3.1 Line Integrals . 16 3.2 Differentiability . 19 3.3 Power Series . 23 3.4 Cauchy's Theorem . 25 3.5 Cauchy's Integral Formula . 27 3.6 Analytic Functions . 30 3.7 Harmonic Functions . 33 3.8 The Maximum Principle . 36 4 Meromorphic Functions and Singularities 37 4.1 Laurent Series . 37 4.2 Isolated Singularities . 40 4.3 The Residue Theorem . 42 4.4 Some Fourier Analysis . 45 4.5 The Argument Principle . 46 5 Complex Mappings 47 5.1 M¨obiusTransformations . 47 5.2 Conformal Mappings . 47 5.3 The Riemann Mapping Theorem . 47 6 Riemann Surfaces 48 6.1 Holomorphic and Meromorphic Maps . 48 6.2 Covering Spaces . 52 7 Elliptic Functions 55 7.1 Elliptic Functions . 55 7.2 Elliptic Curves . 61 7.3 The Classical Jacobian . 67 7.4 Jacobians of Higher Genus Curves . 72 i 0 Introduction 0 Introduction These notes come from a semester course on complex analysis taught by Dr. Richard Carmichael at Wake Forest University during the fall of 2010. The main topics covered include Complex numbers and their properties Complex-valued functions Line integrals Derivatives and power series Cauchy's Integral Formula Singularities and the Residue Theorem The primary reference for the course and throughout these notes is Fisher's Complex Vari- ables, 2nd edition.
    [Show full text]
  • An Introduction to Asymptotic Analysis Simon JA Malham
    An introduction to asymptotic analysis Simon J.A. Malham Department of Mathematics, Heriot-Watt University Contents Chapter 1. Order notation 5 Chapter 2. Perturbation methods 9 2.1. Regular perturbation problems 9 2.2. Singular perturbation problems 15 Chapter 3. Asymptotic series 21 3.1. Asymptotic vs convergent series 21 3.2. Asymptotic expansions 25 3.3. Properties of asymptotic expansions 26 3.4. Asymptotic expansions of integrals 29 Chapter 4. Laplace integrals 31 4.1. Laplace's method 32 4.2. Watson's lemma 36 Chapter 5. Method of stationary phase 39 Chapter 6. Method of steepest descents 43 Bibliography 49 Appendix A. Notes 51 A.1. Remainder theorem 51 A.2. Taylor series for functions of more than one variable 51 A.3. How to determine the expansion sequence 52 A.4. How to find a suitable rescaling 52 Appendix B. Exam formula sheet 55 3 CHAPTER 1 Order notation The symbols , o and , were first used by E. Landau and P. Du Bois- Reymond and areOdefined as∼ follows. Suppose f(z) and g(z) are functions of the continuous complex variable z defined on some domain C and possess D ⊂ limits as z z0 in . Then we define the following shorthand notation for the relative!propertiesD of these functions in the limit z z . ! 0 Asymptotically bounded: f(z) = (g(z)) as z z ; O ! 0 means that: there exists constants K 0 and δ > 0 such that, for 0 < z z < δ, ≥ j − 0j f(z) K g(z) : j j ≤ j j We say that f(z) is asymptotically bounded by g(z) in magnitude as z z0, or more colloquially, and we say that f(z) is of `order big O' of g(z).
    [Show full text]
  • A Note on the Asymptotic Expansion of the Lerch's Transcendent Arxiv
    A note on the asymptotic expansion of the Lerch’s transcendent Xing Shi Cai∗1 and José L. Lópezy2 1Department of Mathematics, Uppsala University, Uppsala, Sweden 2Departamento de Estadística, Matemáticas e Informática and INAMAT, Universidad Pública de Navarra, Pamplona, Spain June 5, 2018 Abstract In [7], the authors derived an asymptotic expansion of the Lerch’s transcendent Φ(z; s; a) for large jaj, valid for <a > 0, <s > 0 and z 2 C n [1; 1). In this paper we study the special case z ≥ 1 not covered in [7], deriving a complete asymptotic expansion of the Lerch’s transcendent Φ(z; s; a) for z > 1 and <s > 0 as <a goes to infinity. We also show that when a is a positive integer, this expansion is convergent for Pm n s <z ≥ 1. As a corollary, we get a full asymptotic expansion for the sum n=1 z =n for fixed z > 1 as m ! 1. Some numerical results show the accuracy of the approximation. 1 Introduction The Lerch’s transcendent (Hurwitz-Lerch zeta function) [2, §25.14(i)] is defined by means of the power series 1 X zn Φ(z; s; a) = ; a 6= 0; −1; −2;:::; (a + n)s n=0 on the domain jzj < 1 for any s 2 C or jzj ≤ 1 for <s > 1. For other values of the variables arXiv:1806.01122v1 [math.CV] 4 Jun 2018 z; s; a, the function Φ(z; s; a) is defined by analytic continuation. In particular [7], 1 Z 1 xs−1e−ax Φ(z; s; a) = −x dx; <a > 0; z 2 C n [1; 1) and <s > 0: (1) Γ(s) 0 1 − ze ∗This work is supported by the Knut and Alice Wallenberg Foundation and the Ministerio de Economía y Competitividad of the spanish government (MTM2017-83490-P).
    [Show full text]
  • Complex Analysis
    COMPLEX ANALYSIS MARCO M. PELOSO Contents 1. Holomorphic functions 1 1.1. The complex numbers and power series 1 1.2. Holomorphic functions 3 1.3. Exercises 7 2. Complex integration and Cauchy's theorem 9 2.1. Cauchy's theorem for a rectangle 12 2.2. Cauchy's theorem in a disk 13 2.3. Cauchy's formula 15 2.4. Exercises 17 3. Examples of holomorphic functions 18 3.1. Power series 18 3.2. The complex logarithm 20 3.3. The binomial series 23 3.4. Exercises 23 4. Consequences of Cauchy's integral formula 25 4.1. Expansion of a holomorphic function in Taylor series 25 4.2. Further consequences of Cauchy's integral formula 26 4.3. The identity principle 29 4.4. The open mapping theorem and the principle of maximum modulus 30 4.5. The general form of Cauchy's theorem 31 4.6. Exercises 36 5. Isolated singularities of holomorphic functions 37 5.1. The residue theorem 40 5.2. The Riemann sphere 42 5.3. Evalutation of definite integrals 43 5.4. The argument principle and Rouch´e's theorem 48 5.5. Consequences of Rouch´e'stheorem 49 5.6. Exercises 51 6. Conformal mappings 53 6.1. Fractional linear transformations 56 6.2. The Riemann mapping theorem 58 6.3. Exercises 63 7. Harmonic functions 65 7.1. Maximum principle 65 7.2. The Dirichlet problem 68 7.3. Exercises 71 8. Entire functions 72 8.1. Infinite products 72 Appunti per il corso Analisi Complessa per i Corsi di Laurea in Matematica dell'Universit`adi Milano.
    [Show full text]
  • Complex Numbers and Functions
    Complex Numbers and Functions Richard Crew January 20, 2018 This is a brief review of the basic facts of complex numbers, intended for students in my section of MAP 4305/5304. I will discuss basic facts of com- plex arithmetic, limits and derivatives of complex functions, power series and functions like the complex exponential, sine and cosine which can be defined by convergent power series. This is a preliminary version and will be added to later. 1 Complex Numbers 1.1 Arithmetic. A complex number is an expression a + bi where i2 = −1. Here the real number a is the real part of the complex number and bi is the imaginary part. If z is a complex number we write <(z) and =(z) for the real and imaginary parts respectively. Two complex numbers are equal if and only if their real and imaginary parts are equal. In particular a + bi = 0 only when a = b = 0. The set of complex numbers is denoted by C. Complex numbers are added, subtracted and multiplied according to the usual rules of algebra: (a + bi) + (c + di) = (a + c) + (b + di) (1.1) (a + bi) − (c + di) = (a − c) + (b − di) (1.2) (a + bi)(c + di) = (ac − bd) + (ad + bc)i (1.3) (note how i2 = −1 has been used in the last equation). Division performed by rationalizing the denominator: a + bi (a + bi)(c − di) (ac − bd) + (bc − ad)i = = (1.4) c + di (c + di)(c − di) c2 + d2 Note that denominator only vanishes if c + di = 0, so that a complex number can be divided by any nonzero complex number.
    [Show full text]
  • Universit`A Degli Studi Di Perugia the Lambert W Function on Matrices
    Universita` degli Studi di Perugia Facolta` di Scienze Matematiche, Fisiche e Naturali Corso di Laurea Triennale in Informatica The Lambert W function on matrices Candidato Relatore MassimilianoFasi BrunoIannazzo Contents Preface iii 1 The Lambert W function 1 1.1 Definitions............................. 1 1.2 Branches.............................. 2 1.3 Seriesexpansions ......................... 10 1.3.1 Taylor series and the Lagrange Inversion Theorem. 10 1.3.2 Asymptoticexpansions. 13 2 Lambert W function for scalar values 15 2.1 Iterativeroot-findingmethods. 16 2.1.1 Newton’smethod. 17 2.1.2 Halley’smethod . 18 2.1.3 K¨onig’s family of iterative methods . 20 2.2 Computing W ........................... 22 2.2.1 Choiceoftheinitialvalue . 23 2.2.2 Iteration.......................... 26 3 Lambert W function for matrices 29 3.1 Iterativeroot-findingmethods. 29 3.1.1 Newton’smethod. 31 3.2 Computing W ........................... 34 3.2.1 Computing W (A)trougheigenvectors . 34 3.2.2 Computing W (A) trough an iterative method . 36 A Complex numbers 45 A.1 Definitionandrepresentations. 45 B Functions of matrices 47 B.1 Definitions............................. 47 i ii CONTENTS C Source code 51 C.1 mixW(<branch>, <argument>) ................. 51 C.2 blockW(<branch>, <argument>, <guess>) .......... 52 C.3 matW(<branch>, <argument>) ................. 53 Preface Main aim of the present work was learning something about a not- so-widely known special function, that we will formally call Lambert W function. This function has many useful applications, although its presence goes sometimes unrecognised, in mathematics and in physics as well, and we found some of them very curious and amusing. One of the strangest situation in which it comes out is in writing in a simpler form the function .
    [Show full text]
  • Chapter 2 Complex Analysis
    Chapter 2 Complex Analysis In this part of the course we will study some basic complex analysis. This is an extremely useful and beautiful part of mathematics and forms the basis of many techniques employed in many branches of mathematics and physics. We will extend the notions of derivatives and integrals, familiar from calculus, to the case of complex functions of a complex variable. In so doing we will come across analytic functions, which form the centerpiece of this part of the course. In fact, to a large extent complex analysis is the study of analytic functions. After a brief review of complex numbers as points in the complex plane, we will ¯rst discuss analyticity and give plenty of examples of analytic functions. We will then discuss complex integration, culminating with the generalised Cauchy Integral Formula, and some of its applications. We then go on to discuss the power series representations of analytic functions and the residue calculus, which will allow us to compute many real integrals and in¯nite sums very easily via complex integration. 2.1 Analytic functions In this section we will study complex functions of a complex variable. We will see that di®erentiability of such a function is a non-trivial property, giving rise to the concept of an analytic function. We will then study many examples of analytic functions. In fact, the construction of analytic functions will form a basic leitmotif for this part of the course. 2.1.1 The complex plane We already discussed complex numbers briefly in Section 1.3.5.
    [Show full text]
  • NM Temme 1. Introduction the Incomplete Gamma Functions Are Defined by the Integrals 7(A,*)
    Methods and Applications of Analysis © 1996 International Press 3 (3) 1996, pp. 335-344 ISSN 1073-2772 UNIFORM ASYMPTOTICS FOR THE INCOMPLETE GAMMA FUNCTIONS STARTING FROM NEGATIVE VALUES OF THE PARAMETERS N. M. Temme ABSTRACT. We consider the asymptotic behavior of the incomplete gamma func- tions 7(—a, —z) and r(—a, —z) as a —► oo. Uniform expansions are needed to describe the transition area z ~ a, in which case error functions are used as main approximants. We use integral representations of the incomplete gamma functions and derive a uniform expansion by applying techniques used for the existing uniform expansions for 7(0, z) and V(a,z). The result is compared with Olver's uniform expansion for the generalized exponential integral. A numerical verification of the expansion is given. 1. Introduction The incomplete gamma functions are defined by the integrals 7(a,*)= / T-Vcft, r(a,s)= / t^e^dt, (1.1) where a and z are complex parameters and ta takes its principal value. For 7(0,, z), we need the condition ^Ra > 0; for r(a, z), we assume that |arg2:| < TT. Analytic continuation can be based on these integrals or on series representations of 7(0,2). We have 7(0, z) + r(a, z) = T(a). Another important function is defined by 7>>*) = S7(a,s) = =^ fu^e—du. (1.2) This function is a single-valued entire function of both a and z and is real for pos- itive and negative values of a and z. For r(a,z), we have the additional integral representation e-z poo -zt j.-a r(a, z) = — r / dt, ^a < 1, ^z > 0, (1.3) L [i — a) JQ t ■+■ 1 which can be verified by differentiating the right-hand side with respect to z.
    [Show full text]