Integral Formula for the Bessel Function of the First Kind

Total Page:16

File Type:pdf, Size:1020Kb

Integral Formula for the Bessel Function of the First Kind INTEGRAL FORMULA FOR THE BESSEL FUNCTION OF THE FIRST KIND ENRICO DE MICHELI ABSTRACT. In this paper, we prove a new integral representation for the Bessel function of the first kind Jµ (z), which holds for any µ,z C. ∈ 1. INTRODUCTION In this note, we prove the following integral representation of the Bessel function of the first kind Jµ (z), which holds for unrestricted complex values of the order µ: z µ π 1 izcosθ 1 iθ C C (1.1) Jµ (z)= e γ∗ µ, 2 ize dθ (µ ;z ), 2π 2 π ∈ ∈ Z− where γ∗ denotes the Tricomi version of the (lower) incomplete gamma function. Lower and upper incomplete gamma functions arise from the decomposition of the Euler integral for the gamma function: w t µ 1 (1.2a) γ(µ,w)= e− t − dt (Re µ > 0), 0 Z ∞ t µ 1 (1.2b) Γ(µ,w)= e− t − dt ( argw < π). w | | Z Analytical continuation with respect to both parameters can be based either on the integrals in (1.2) or on series expansions of γ(µ,w), e.g. [4, Eq. 9.2(4)]: ∞ ( 1)n wµ+n (1.3) γ(µ,w)= ∑ − , n=0 n! µ + n which shows that γ(µ,w) has simple poles at µ = 0, 1, 2,... and, in general, is a multi- valued (when µ is not a positive integer) function with− a− branch point at w = 0. Similarly, even Γ(µ,w) is in general multi-valued but it is an entire function of µ. Inconveniences re- lated to poles and multi-valuedness of γ can be circumvented by introducing what Tricomi arXiv:2012.02887v1 [math.CA] 4 Dec 2020 called the fundamental function of the incomplete gamma theory [11, §2, p. 264]: µ w w w− t µ 1 e− (1.4) γ∗(µ,w)= e− t − dt = M(1,1 + µ;w), Γ(µ) Γ(1 + µ) Z0 where M(a,b;z) is the Kummer confluent hypergeometric function [8, Eq. 13.2.2]. The function γ (µ,w) is entire in both w and µ, real-valued for real w and µ and γ ( n,w)= wn ∗ ∗ − for n = 0,1,2,... B¨ohmer [1] and later Tricomi [10, 11] (see also [5]) assumed γ∗(µ,w) as base point for the theory, defining the incomplete gamma functions by the relations: µ (1.5) γ(µ,w)= Γ(µ)w γ∗(µ,w), µ (1.6) Γ(µ,w)= Γ(µ)[1 w γ∗(µ,w)]. − 2010 Mathematics Subject Classification. 33C10,40C10,33B20. Key words and phrases. Bessel functions, Incomplete gamma functions, Integral representations. 1 2 ENRICO DE MICHELI For a detailed account of the incomplete gamma functions theory, see [5], [9, Chapter 11] and [4, 8], where massive compilations of results are stated. For what concerns details on the Bessel function, the reader is referred to the classical treatise of Watson [12]. The strong link between Bessel functions and incomplete gamma functions has long been known. Examples of this connection are represented, e.g., by the integral representa- tions of γ(µ,w) in terms of Jµ (z) for Re µ > 0 [4, Eq. 9.3(4)] and of Γ(µ,w) in terms of modified Bessel function of the second kind Kµ (z) for Re µ < 1 [8, Eq. 8.6.6]. In this re- spect, it is worth mentioning also the expansions of incomplete gamma functions in series of Bessel functions: [5, Eqs. (2.8), (2.9)], [11, Eqs. (40)], and of modified Bessel func- tions: [8, Eqs. 8.7.4, 8.7.5], [10, Eqs. (48), (49)]. Formula (1.1), whose validity extends to (µ,z) C2, makes even more explicit the intimate link between these two classes of functions.∈ Classical formulae of the Bessel function of the first kind Jµ (z) given by integrals on the real line can be found in [8, Chapter 10.9], [4, Sect. 7.3] and [6, Sect. 8.41]. Many of them hold only for restricted values of µ: for instance, the classical Bessel integral [7, p. 70], µ π ( i) izcosθ (1.7) Jµ (z)= − e cos µθ dθ (µ Z), π 0 ∈ Z and Poisson’s integral [7, Eq. (14.5), p. 38], µ π 1 z izcosθ 2µ 1 (1.8) Jµ (z)= 1 e (sinθ) dθ (Re µ > 2 ), Γ(µ + )√π 2 0 − 2 Z 1 hold only for m Z and Re µ > 2 , respectively. Additional examples are the Gubler ∈ − 1 representations [4, Eq. 7.3(11)], which holds for Re µ < 2 , the Mehler-Sonine integral 1 formulae for Re µ < 2 [4, Eq. 7.12(12)] or for Re µ < 1 [4, Eq. 7.12(14)]. Only the well-known Schl¨afli’s| | and Heine’s representations| hold fo| r unrestricted complex values of the index µ. Schl¨afli’s representation [4, Eq. 7.3(9)] holds for µ C, with Rez > 0 (it holds also in case Rez = 0 provided that Re µ > 0), while the two Heine’s∈ expressions [4, Eqs. 7.3(31) & (32)] hold for µ C, but separately in the upper and lower half-plane of the complex z-plane, respectively.∈ In both cases, Schl¨afli’s and Heine’s, the representations are made up of the sum of two integrals. Standard generalizations of Poisson’s integral are obtained via its extension to complex contour integrals, examples of which are the classical Hankel’s representations [7, Eqs. (19.7), (19.8)]. As Poisson’s representation of Jµ (z) can be seen as the generalization of 1 Bessel’s integral to complex values of µ with Re µ > 2 , mere visual inspection shows that formula (1.1) can be viewed as the generalization− of Poisson’s integral to represent 2µ Jµ (z) at values of µ belonging to the entire complex plane, the role of (sinθ) in (1.8) 1 iθ being now taken by the incomplete gamma function γ∗(µ, 2 ize ). The proof of formula (1.1) is given in the next section and, essentially, makes use of elementary properties of the special functions involved. A preliminary result was given in [3]. 2. INTEGRAL REPRESENTATION FOR Jµ (z) We denote by N0 = 0,1,2,... the set of non-negative integers. { } Theorem 2.1. The following integral representation of the Bessel function of the first kind holds true for any complex order µ C and for z C (slit along the real negative axis ∈ ∈ INTEGRALFORMULAFORBESSELFUNCTION 3 when µ Z): 6∈ µ π 1 z izcosθ 1 iθ C C (2.1) Jµ (z)= e γ∗ µ, 2 ize dθ (µ ;z ). 2π 2 π ∈ ∈ Z− Proof. Our starting point is the Gegenbauer generalization of Poisson’s integral represen- tation of the Bessel functions of the first kind [12, §3.32, Eq. (1), p. 50], which holds for 1 ℓ N0 and Reν > : ∈ − 2 ℓ ν π ( i) Γ(2ν)ℓ! z izcosu 2ν ν (2.2) Jν+ℓ(z)= −1 1 e (sinu) Cℓ (cosu)du, Γ(ν + )Γ( )Γ(2ν + ℓ) 2 0 2 2 Z ν where Cℓ (t) is the Gegenbauer polynomial of order ν and degree ℓ. Now, we recall the ν following representation of the Gegenbauer polynomials Cℓ (cosu), holding for Reν > 0 and ℓ N0, that we proved in [2, Proposition 1]: (2.3) ∈ iπν 2π u ν Γ(ℓ + 2ν) e− − i(ℓ+ν)t ν 1 Cℓ (cosu)= ν 2 e (cosu cost) − dt (u [0,π]). 2 ℓ![Γ(ν)] sinu (2ν 1) u − ∈ ( ) − Z Plugging (2.3) into (2.2) and using the Legendre duplication formula for the gamma func- tion, we obtain: ℓ iπν ν π 2π u ( i) e− z izcosu − i(ℓ+ν)θ ν 1 (2.4) Jν+ℓ(z)= − du sinue e (cosu cosθ) − dθ. 2π Γ(ν) 0 u − Z Z Interchanging the order of integration, (2.4) can be written as follows: ℓ iπν ν π θ ( i) e− z i(ℓ+ν)θ izcosu ν 1 Jν+ℓ(z)= − dθ e e (cosu cosθ) − sinudu 2π Γ(ν) 0 0 − (2.5) Z Z 2π 2π θ i(ℓ+ν)θ − izcosu ν 1 + dθ e e (cosu cosθ) − sinudu . π 0 − Z Z Next, changing the integration variables: θ 2π θ and u u, the second integral on the r.h.s. of (2.5) becomes: − → →− 0 θ iν2π i(ℓ+ν)θ izcosu ν 1 e dθ e e (cosu cosθ) − sinudu, π 0 − Z− Z which, inserted in (2.5), yields: ℓ ν 0 θ ( i) z iνπ i(ℓ+ν)θ izcosu ν 1 Jν+ℓ(z)= − e dθ e e (cosu cosθ) − sinudu 2π Γ(ν) π 0 − (2.6) Z− Z π θ iνπ i(ℓ+ν)θ izcosu ν 1 +e− dθ e e (cosu cosθ) − sinudu . 0 0 − Z Z Formula (2.6) allows us to write Jν+ℓ(z) as the Fourier coefficient of a suitable 2π-periodic function, i.e.: π ℓ iℓθ (2.7) Jν+ℓ(z) = ( i) Aν,z(θ) e dθ (ℓ N0,Reν > 0), − π ∈ Z− where Aν,z(θ) denotes the Abel-type integral ν 1 z iν[θ π sgn(θ)] izt ν 1 (2.8) A (θ)= e − e (t cosθ) − dt (Reν > 0), ν,z 2πΓ(ν) Zcosθ − 4 ENRICO DE MICHELI and sgn( ) denotes the sign function. Changing in (2.8) the integration variable iz(t · − − cosθ) t and recalling (1.2a), Aν,z(θ) can be rewritten as follows: → ν i iν[θ π sgn(θ)] izcosθ (2.9) Aν,z(θ)= e − e P(ν, iz(1 cosθ)) (Reν > 0), 2π − − .
Recommended publications
  • Branch Points and Cuts in the Complex Plane
    BRANCH POINTS AND CUTS IN THE COMPLEX PLANE Link to: physicspages home page. To leave a comment or report an error, please use the auxiliary blog and include the title or URL of this post in your comment. Post date: 6 July 2021. We’ve looked at contour integration in the complex plane as a technique for evaluating integrals of complex functions and finding infinite integrals of real functions. In some cases, the complex functions that are to be integrated are multi- valued. As a preliminary to the contour integration of such functions, we’ll look at the concepts of branch points and branch cuts here. The stereotypical function that is used to introduce branch cuts in most books is the complex logarithm function logz which is defined so that elogz = z (1) If z is real and positive, this reduces to the familiar real logarithm func- tion. (Here I’m using natural logs, so the real natural log function is usually written as ln. In complex analysis, the term log is usually used, so be careful not to confuse it with base 10 logs.) To generalize it to complex numbers, we write z in modulus-argument form z = reiθ (2) and apply the usual rules for taking a log of products and exponentials: logz = logr + iθ (3) = logr + iargz (4) To see where problems arise, suppose we start with z on the positive real axis and increase θ. Everything is fine until θ approaches 2π. When θ passes 2π, the original complex number z returns to its starting value, as given by 2.
    [Show full text]
  • Riemann Surfaces
    RIEMANN SURFACES AARON LANDESMAN CONTENTS 1. Introduction 2 2. Maps of Riemann Surfaces 4 2.1. Defining the maps 4 2.2. The multiplicity of a map 4 2.3. Ramification Loci of maps 6 2.4. Applications 6 3. Properness 9 3.1. Definition of properness 9 3.2. Basic properties of proper morphisms 9 3.3. Constancy of degree of a map 10 4. Examples of Proper Maps of Riemann Surfaces 13 5. Riemann-Hurwitz 15 5.1. Statement of Riemann-Hurwitz 15 5.2. Applications 15 6. Automorphisms of Riemann Surfaces of genus ≥ 2 18 6.1. Statement of the bound 18 6.2. Proving the bound 18 6.3. We rule out g(Y) > 1 20 6.4. We rule out g(Y) = 1 20 6.5. We rule out g(Y) = 0, n ≥ 5 20 6.6. We rule out g(Y) = 0, n = 4 20 6.7. We rule out g(C0) = 0, n = 3 20 6.8. 21 7. Automorphisms in low genus 0 and 1 22 7.1. Genus 0 22 7.2. Genus 1 22 7.3. Example in Genus 3 23 Appendix A. Proof of Riemann Hurwitz 25 Appendix B. Quotients of Riemann surfaces by automorphisms 29 References 31 1 2 AARON LANDESMAN 1. INTRODUCTION In this course, we’ll discuss the theory of Riemann surfaces. Rie- mann surfaces are a beautiful breeding ground for ideas from many areas of math. In this way they connect seemingly disjoint fields, and also allow one to use tools from different areas of math to study them.
    [Show full text]
  • The Monodromy Groups of Schwarzian Equations on Closed
    Annals of Mathematics The Monodromy Groups of Schwarzian Equations on Closed Riemann Surfaces Author(s): Daniel Gallo, Michael Kapovich and Albert Marden Reviewed work(s): Source: Annals of Mathematics, Second Series, Vol. 151, No. 2 (Mar., 2000), pp. 625-704 Published by: Annals of Mathematics Stable URL: http://www.jstor.org/stable/121044 . Accessed: 15/02/2013 18:57 Your use of the JSTOR archive indicates your acceptance of the Terms & Conditions of Use, available at . http://www.jstor.org/page/info/about/policies/terms.jsp . JSTOR is a not-for-profit service that helps scholars, researchers, and students discover, use, and build upon a wide range of content in a trusted digital archive. We use information technology and tools to increase productivity and facilitate new forms of scholarship. For more information about JSTOR, please contact [email protected]. Annals of Mathematics is collaborating with JSTOR to digitize, preserve and extend access to Annals of Mathematics. http://www.jstor.org This content downloaded on Fri, 15 Feb 2013 18:57:11 PM All use subject to JSTOR Terms and Conditions Annals of Mathematics, 151 (2000), 625-704 The monodromy groups of Schwarzian equations on closed Riemann surfaces By DANIEL GALLO, MICHAEL KAPOVICH, and ALBERT MARDEN To the memory of Lars V. Ahlfors Abstract Let 0: 7 (R) -* PSL(2, C) be a homomorphism of the fundamental group of an oriented, closed surface R of genus exceeding one. We will establish the following theorem. THEOREM. Necessary and sufficient for 0 to be the monodromy represen- tation associated with a complex projective stucture on R, either unbranched or with a single branch point of order 2, is that 0(7ri(R)) be nonelementary.
    [Show full text]
  • Sequential Discontinuities of Feynman Integrals and the Monodromy Group
    Sequential Discontinuities of Feynman Integrals and the Monodromy Group Jacob L. Bourjaily1,2, Holmfridur Hannesdottir3, Andrew J. McLeod1, Matthew D. Schwartz3, and Cristian Vergu1 1Niels Bohr International Academy and Discovery Center, Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, DK-2100, Copenhagen Ø, Denmark 2Institute for Gravitation and the Cosmos, Department of Physics, Pennsylvania State University, University Park, PA 16892, USA 3Department of Physics, Harvard University, Cambridge, MA 02138, USA July 29, 2020 Abstract We generalize the relation between discontinuities of scattering amplitudes and cut diagrams to cover sequential discontinuities (discontinuities of discontinuities) in arbitrary momentum channels. The new relations are derived using time-ordered per- turbation theory, and hold at phase-space points where all cut momentum channels are simultaneously accessible. As part of this analysis, we explain how to compute se- quential discontinuities as monodromies and explore the use of the monodromy group in characterizing the analytic properties of Feynman integrals. We carry out a number of cross-checks of our new formulas in polylogarithmic examples, in some cases to all loop orders. arXiv:2007.13747v1 [hep-th] 27 Jul 2020 Contents 1 Introduction1 2 Cutting rules: a review4 2.1 Cutkosky, 't Hooft and Veltman . .5 2.2 Time-ordered perturbation theory . .7 3 Discontinuities 10 3.1 Covariant approach . 11 3.2 Discontinuities in TOPT . 14 4 Discontinuities as monodromies 17 4.1 Warm-up: the natural logarithm . 17 4.2 The monodromy group . 21 4.3 Monodromies of propagators . 31 5 Sequential discontinuities 33 5.1 Sequential discontinuities in the same channel . 34 5.2 Sequential discontinuities in different channels .
    [Show full text]
  • CHAPTER 4 Elementary Functions Dr. Pulak Sahoo
    CHAPTER 4 Elementary Functions BY Dr. Pulak Sahoo Assistant Professor Department of Mathematics University Of Kalyani West Bengal, India E-mail : sahoopulak1@gmail:com 1 Module-3: Multivalued Functions-I 1 Introduction We recall that a function w = f(z) is called multivalued function if for all or some z of the domain, we find more than one value of w. Thus, a function f is said to be single valued if f satisfies f(z) = f[z(r; θ)] = f[z(r; θ + 2π)]: Otherwise, f is called as a multivalued function. We know that a multivalued function can be considered as a collection of single valued functions. For analytical properties of a multivalued function, we just consider those domains in which the functions are single valued. Branch By branch of a multivalued function f(z) defined on a domain D1 we mean a single valued function g(z) which is analytic in some sub-domain D ⊂ D1 at each point of which g(z) is one of the values of f(z). Branch Points and Branch Lines A point z = α is called a branch point of a multivalued function f(z) if the branches of f(z) are interchanged when z describes a closed path about α. We consider the function w = z1=2. Suppose that we allow z to make a complete circuit around the origin in the anticlockwise sense starting from the point P. Thus we p p have z = reiθ; w = reiθ=2; so that at P, θ = β and w = reiβ=2: After a complete p p circuit back to P, θ = β + 2π and w = rei(β+2π)=2 = − reiβ=2: Thus we have not achieved the same value of w with which we started.
    [Show full text]
  • Riemann Surfaces
    Riemann Surfaces See Arfken & Weber section 6.7 (mapping) or [1][sections 106-108] for more information. 1 Definition A Riemann surface is a generalization of the complex plane to a surface of more than one sheet such that a multiple-valued function on the original complex plane has only one value at each point of the surface. Briefly, the notion of a Riemann surface is important whenever considering functions with branch cuts. Perhaps the easiest example is the Riemann surface for log z. Write z = r exp(iθ), then log z = (log r)+ iθ. However, the θ in that expression is not uniquely defined – it is only defined up to a multiple of 2π. As you walk around the origin of the complex plane, the function log z comes back to itself up to additions of 2πi. We can construct a cover of the complex plane on which log z is single-valued by taking a helix over the complex plane – going 360◦ about the origin takes from you from sheet of the cover to another sheet. Put another way, every time you go through the branch cut, you go to a different sheet. We can construct such a cover as follows. (Our discussion here is verbatim from [1][section 106].) Consider the complex z plane, with the origin deleted, as a sheet R0 which is cut along the positive half of the real axis. On that sheet, let θ range from 0 to 2π. Let a second sheet R1 be cut in the same way and placed in front of the sheet R0.
    [Show full text]
  • Full Text (PDF Format)
    Methods and Applications of Analysis © 1998 International Press 5 (4) 1998, pp. 425-438 ISSN 1073-2772 ON THE RESURGENCE PROPERTIES OF THE UNIFORM ASYMPTOTIC EXPANSION OF THE INCOMPLETE GAMMA FUNCTION A. B. Olde Daalhuis ABSTRACT. We examine the resurgence properties of the coefficients 07.(77) ap- pearing in the uniform asymptotic expansion of the incomplete gamma function. For the coefficients cr(T7), we give an asymptotic approximation as r -» 00 that is a sum of two incomplete beta functions plus a simple asymptotic series in which the coefficients are again Cm(rj). The method of this paper is based on the Borel-Laplace transform, which means that next to the asymptotic approximation of CrC??), we also obtain an exponentially-improved asymptotic expansion for the incomplete gamma function. 1. Introduction and summary Recently it has been shown in many papers that the incomplete gamma function is the basis function for exponential asymptotics. The asymptotic expansion of r(a, z) as z —> 00 and a fixed is rather simple. However, in exponential asymptotics, we need the asymptotic properties of r(a, z) as a —> 00 and z = Aa, A 7^ 0, a complex constant. In [9] and [11], the uniform asymptotic expansions of the (normalized) incomplete gamma function as a —> 00 are given as r(« + ll^-rf-a,2e-") ~ ^erfc(-itijia) + >-,=; ^MlX-or' (Lib) where 77 is defined as 1 77=(2(A-l-lnA))t (1.2) and z x= -. (1.3) a The complementary error function in (1.1b) will give the Stokes smoothing in expo- nentially improved asymptotics for integrals and solutions of differential equations.
    [Show full text]
  • "According to Abramowitz and Stegun" Or Arccoth Needn't Be Uncouth
    "According to Abramowitz and Stegun" or arccoth needn't be uncouth Robert M. Corless James H. Davenport* David J. Jeffrey Dept. Mathematical Sciences Stephen M. Watt University of Bath Ontario Research Centre Bath BA2 7AY for Computer Algebra England . orcca.. On. Ca. jhd@maths, bath. ac, uk Abstract functions, as well as relations between these functions, and relations between them and the forward functions. This paper ~ addresses the definitions in OpenMath of the We discuss the implications of accurate translation on elementary functions. The original OpenMath definitions, the design of the phrasebooks [5] translating between Open- like most other sources, simply cite [2] as the definition. We Math and actual systems. Note that OpenMath does not show that this is not adequate, and propose precise defi- say that one definition is 'right' and another 'wrong': it nitions, and explore the relationships between these defini- merely provides a lingua francs for passing semantically ac- tions. curate representations between systems. Semantically cor- In particular, we introduce the concept of a couth pair rect phrasebooks would deduce that of definitions, e.g. of arcsin and arcsinh, and show t:hat the pair axccot and arccoth can be couth. z= Maple Derive 1 Introduction Notation. Throughout this paper, we use arccot etc to Definitions of the elementary functions are given in many mean the precise function definitions we are using, and vari- textbooks and mathematical tables, such as [2, 7, 16]. How- ants are indicated by annotations as in the equation above. ever, it is a sad fact that these definitions often require a Throughout, z and its decorations indicate complex vari- great deal of common sense to interpret them, or, to be ables, x and Y real variables.
    [Show full text]
  • Zeros of the Dilogarithm
    Zeros of the dilogarithm Cormac O’Sullivan∗ August 21, 2018 Abstract We show that the dilogarithm has at most one zero on each branch, that each zero is close to a root of unity, and that they may be found to any precision with Newton’s method. This work is motivated by ap- plications to the asymptotics of coefficients in partial fraction decompositions considered by Rademacher. We also survey what is known about zeros of polylogarithms in general. 1 Introduction In the recent resolution of an old conjecture of Rademacher, described below in Section 1.2, the location of a particular zero, w0, of the dilogarithm played an important role. It has been known since [LR00] that the only zero of the dilogarithm on its principal branch is at 0. The zero w0 is on the next branch. Zeros on further branches were also needed in [O’Sa] and in this paper we locate all zeros on every branch. The dilogarithm is initially defined as ∞ zn Li (z) := for z 6 1, (1.1) 2 n2 | | n=1 X see for example [Max03, Zag07], with an analytic continuation given by z du log(1 u) . (1.2) − − u Z0 The principal branch of the logarithm has π < arg z 6 π with a branch cut ( , 0]. From (1.2), the − −∞ corresponding principal branch of the dilogarithm has branch points at 1, and branch cut [1, ). Crossing ∞ ∞ this branch cut from below, it is easy to show that the dilogarithm increases by 2πi log(z) over its principal value. On this new sheet there is now an additional branch point at 0 coming from the logarithm.
    [Show full text]
  • Gamma Function and Bessel Functions - Lecture 7
    Gamma Function and Bessel Functions - Lecture 7 1 Introduction - Gamma Function The Gamma function is defined by; ∞ t z 1 Γ(z) = dt e− t − R0 Here, z can be a complex, non-integral number. When z = n, an integer, integration by parts produces the factorial; Γ(n)=(n 1)! − In order for the integral to converge, Rez > 0. However it may be extended to negative values of Re z by the recurrence relation; zΓ(z) = Γ(z + 1) This diverges for n a negative integer. Another representation of the Gamma function is the infinite product; γz ∞ z/n 1/Γ(z) = ze (1 + z/n) e− nQ=1 γ = Γ′(1)/Γ(1) = 0.5772157 − The term γ is the Euler-Mascheroni constant. 2 Analytic properties of the Gamma function For Rez > 0, Γ(z) is finite and its derivative is finite. Thus Γ(z) is analytic when z > 0. Suppose Re (z + n + 1) > 0. For the smallest possible value of the integer, n ℜ Γ(z + n + 1) Γ(z) = (z + n)(z + n 1) z − ··· Thus Γ(z) is expressed in terms of an analytic function Γ(z + n + 1) and a set of poles when z =0, 1, Therefore Γ(z) is analytic except for poles along the negative real axis. Near the pole− set···z = n + ǫ and consider ǫ 0. The residue is ( 1)n/n! which is similar to the poles of 1/sin−(πz) except that this→ function has poles along− the positive real axis as 1 well as the negative real axis.
    [Show full text]
  • 2.4 Branches of Functions
    2.4 I Branches of Functions 79 2.4BRANCHESOFFUNCTIONS In Section 2.2 we defined the principal square root function and investigated some of its properties. We left unanswered some questions concerning the choices of square roots. We now look at these questions because they are similar to situations involving other elementary functions. 80 Chapter 2 I Complex Functions Inour definitionofa functioninSection2.1,we specified that each value of the independent variable in the domain is mapped onto one and only one value in the range. As a result, we often talk about a single-valued function, which emphasizes the “only one” part of the definition and allows us to distinguish such functions from multiple-valued functions, which we now introduce. Let w = f (z) denote a function whose domain is the set D and whose range is the set R.Ifw is a value inthe range,thenthere is anassociated inverse relation z = g (w) that assigns to each value w the value (or values) of z in D for which the equation f (z)=w holds. But unless f takes onthe value w at most once in D, thenthe inverse relation g is necessarily many valued, and we say that g is a multivalued function. For example, the inverse of the function 2 1 w = f (z)=z is the square root function z = g (w)=w 2 . For each value z other than z = 0, then, the two points z and −z are mapped onto the same point w = f (z); hence g is, in general, a two-valued function.
    [Show full text]
  • Properties and Computation of the Inverse of the Gamma Function
    Western University Scholarship@Western Electronic Thesis and Dissertation Repository 4-23-2018 11:30 AM Properties and Computation of the Inverse of the Gamma function Folitse Komla Amenyou The University of Western Ontario Supervisor Jeffrey,David The University of Western Ontario Graduate Program in Applied Mathematics A thesis submitted in partial fulfillment of the equirr ements for the degree in Master of Science © Folitse Komla Amenyou 2018 Follow this and additional works at: https://ir.lib.uwo.ca/etd Part of the Numerical Analysis and Computation Commons Recommended Citation Amenyou, Folitse Komla, "Properties and Computation of the Inverse of the Gamma function" (2018). Electronic Thesis and Dissertation Repository. 5365. https://ir.lib.uwo.ca/etd/5365 This Dissertation/Thesis is brought to you for free and open access by Scholarship@Western. It has been accepted for inclusion in Electronic Thesis and Dissertation Repository by an authorized administrator of Scholarship@Western. For more information, please contact [email protected]. Abstract We explore the approximation formulas for the inverse function of Γ. The inverse function of Γ is a multivalued function and must be computed branch by branch. We compare ˇ three approximations for the principal branch Γ0. Plots and numerical values show that the choice of the approximation depends on the domain of the arguments, specially for small arguments. We also investigate some iterative schemes and find that the Inverse Quadratic Interpolation scheme is better than Newton's scheme for improving the initial approximation. We introduce the contours technique for extending a real-valued function into the complex plane using two examples from the elementary functions: the log and the arcsin functions.
    [Show full text]