The Gamma function. Marco Bonvini. October 9, dt e t t z 1. (1) Γ(z + 1) = z Γ(z) : (2) = e t t z. + z dt e t t z 1. = z Γ(z).

Similar documents
1 6 = 1 6 = + Factorials and Euler s Gamma function

APPM 4360/5360 Exam #2 Solutions Spring 2015

Physics 116A Solutions to Homework Set #9 Winter 2012

Chapter 8. Euler s Gamma function

The Gamma function Michael Taylor. Abstract. This material is excerpted from 18 and Appendix J of [T].

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

Chapter 8. Euler s Gamma function

j=1 dz Res(f, z j ) = 1 d k 1 dz k 1 (z c)k f(z) Res(f, c) = lim z c (k 1)! Res g, c = f(c) g (c)

... and realizing that as n goes to infinity the two integrals should be equal. This yields the Wallis result-

Notes 8 Singularities

Topics in Probability Theory and Stochastic Processes Steven R. Dunbar. Stirling s Formula Derived from the Gamma Function

Lecture 2 Appendix B: Some sample problems from Boas, Chapter 1. Solution: We want to use the general expression for the form of a geometric series

The natural exponential function

Math 2784 (or 2794W) University of Connecticut

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

Dirichlet s Theorem on Arithmetic Progressions

COMPUTING THE EULER S CONSTANT: A HISTORICAL OVERVIEW OF ALGORITHMS AND RESULTS

SCORE. Exam 2. MA 114 Exam 2 Fall 2017

SOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.

18.440, March 9, Stirling s formula

Introductions to HarmonicNumber2

1 Approximating Integrals using Taylor Polynomials

AP Calculus Chapter 9: Infinite Series

Definition of z-transform.

of the Barnes function.

PROPERTIES OF THE POSITIVE INTEGERS

f(w) w z =R z a 0 a n a nz n Liouville s theorem, we see that Q is constant, which implies that P is constant, which is a contradiction.

Worksheet on Generating Functions

The Riemann Zeta Function

ExpIntegralEi. Notations. Primary definition. Specific values. Traditional name. Traditional notation. Mathematica StandardForm notation

6.3 Testing Series With Positive Terms

f(x) dx as we do. 2x dx x also diverges. Solution: We compute 2x dx lim

Bernoulli Numbers. n(n+1) = n(n+1)(2n+1) = n(n 1) 2

Solutions for Math 411 Assignment #8 1

REGULARIZATION OF CERTAIN DIVERGENT SERIES OF POLYNOMIALS

Dupuy Complex Analysis Spring 2016 Homework 02

s = and t = with C ij = A i B j F. (i) Note that cs = M and so ca i µ(a i ) I E (cs) = = c a i µ(a i ) = ci E (s). (ii) Note that s + t = M and so

(A sequence also can be thought of as the list of function values attained for a function f :ℵ X, where f (n) = x n for n 1.) x 1 x N +k x N +4 x 3

Additional Notes on Power Series

INFINITE SEQUENCES AND SERIES

Carleton College, Winter 2017 Math 121, Practice Final Prof. Jones. Note: the exam will have a section of true-false questions, like the one below.

THE LEGENDRE POLYNOMIALS AND THEIR PROPERTIES. r If one now thinks of obtaining the potential of a distributed mass, the solution becomes-

The Z-Transform. (t-t 0 ) Figure 1: Simplified graph of an impulse function. For an impulse, it can be shown that (1)

Sequences and Limits

Sequences of Definite Integrals, Factorials and Double Factorials

Advanced Analysis. Min Yan Department of Mathematics Hong Kong University of Science and Technology

Series: Infinite Sums

Quadratic Transformations of Hypergeometric Function and Series with Harmonic Numbers

Analytic Continuation

Midterm Exam #2. Please staple this cover and honor pledge atop your solutions.

Chapter 7: The z-transform. Chih-Wei Liu

Harmonic Number Identities Via Euler s Transform

Fourier Series and their Applications

Solutions to Homework 1

TR/46 OCTOBER THE ZEROS OF PARTIAL SUMS OF A MACLAURIN EXPANSION A. TALBOT

Lecture 4 Conformal Mapping and Green s Theorem. 1. Let s try to solve the following problem by separation of variables

A New Method to Order Functions by Asymptotic Growth Rates Charlie Obimbo Dept. of Computing and Information Science University of Guelph

Discrete Orthogonal Moment Features Using Chebyshev Polynomials

3 Gauss map and continued fractions

PAPER : IIT-JAM 2010

1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and second-order partial derivatives for ( x0, y

Notes 12 Asymptotic Series

ENGI Series Page 6-01

Chapter 6 Infinite Series

The Arakawa-Kaneko Zeta Function

An application of the Hooley Huxley contour

CHAPTER 10 INFINITE SEQUENCES AND SERIES

COMM 602: Digital Signal Processing

NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE

Topics in Probability Theory and Stochastic Processes Steven R. Dunbar. Evaluation of the Gaussian Density Integral

Some integrals related to the Basel problem

Informal Notes: Zeno Contours, Parametric Forms, & Integrals. John Gill March August S for a convex set S in the complex plane.

Evaluation of Some Non-trivial Integrals from Finite Products and Sums

Math 155 (Lecture 3)

Series with Central Binomial Coefficients, Catalan Numbers, and Harmonic Numbers

FUNCTIONAL EQUATIONS WITH PRIME ROOTS FROM ARITHMETIC EXPRESSIONS FOR. BARRY BRENT Efmhurst, New York 11373

Bernoulli Polynomials Talks given at LSBU, October and November 2015 Tony Forbes

Topics in Probability Theory and Stochastic Processes Steven R. Dunbar. Stirling s Formula from the Sum of Average Differences

Appendix: The Laplace Transform

A PROOF OF THE TWIN PRIME CONJECTURE AND OTHER POSSIBLE APPLICATIONS

De la Vallée Poussin Summability, the Combinatorial Sum 2n 1

1 Statement of the Game

Complex Numbers Solutions

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 2 9/9/2013. Large Deviations for i.i.d. Random Variables

C. Complex Numbers. x 6x + 2 = 0. This equation was known to have three real roots, given by simple combinations of the expressions

Q-BINOMIALS AND THE GREATEST COMMON DIVISOR. Keith R. Slavin 8474 SW Chevy Place, Beaverton, Oregon 97008, USA.

Subject: Differential Equations & Mathematical Modeling-III

(I.C) THE DISTRIBUTION OF PRIMES

Product measures, Tonelli s and Fubini s theorems For use in MAT3400/4400, autumn 2014 Nadia S. Larsen. Version of 13 October 2014.

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

Assignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.

Solutions to Exercises 4.1

Beyond simple iteration of a single function, or even a finite sequence of functions, results

A MODIFIED BERNOULLI NUMBER. D. Zagier

7 Sequences of real numbers

Series III. Chapter Alternating Series

Physics 116A Solutions to Homework Set #1 Winter Boas, problem Use equation 1.8 to find a fraction describing

We are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n

Lecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods

MATH 1080: Calculus of One Variable II Fall 2017 Textbook: Single Variable Calculus: Early Transcendentals, 7e, by James Stewart.

Transcription:

The Gamma fuctio Marco Bovii October 9, 2 Gamma fuctio The Euler Gamma fuctio is defied as Γ() It is easy to show that Γ() satisfy the recursio relatio ideed, itegratig by parts, dt e t t. () Γ( + ) Γ() : (2) Γ( + ) dt e t t e t t + dt e t t Γ(). (3) For iteger eq. (2) is the recursio relatio of the factorial, ad thus we have Γ( + ) ; (4) because i additio Γ() (easly derived from the defiitio), we have the idetificatio Γ( + ), (5) ad i this sese the Gamma fuctio is a complex extesio of the factorial. The sequece of Gamma fuctio computed i all half-itegers ca be obtaied usig subsequetly the recursio relatio (2) ad kowig that ( Γ 2) π (6) that is easy to compute: ( Γ 2) 2 2 dt e t t 2 dx e x2 π 2 π. (7)

5 5 2 2 4 5 Figure : Re Γ(x) for real x.. Aalytical structure First, from the defiitio (), we see that Γ( ) Γ(), (8) that is to say that Gamma is a real fuctio, ad i particular Im Γ(x) for x R. Ivertig eq. (2) we have Γ() Γ( + ), (9) ad whe it diverges (because Γ() is fiite). It is evidet that correspods to a simple pole of order ad residue : ideed Res Γ() Γ() Γ( + ). () Formally, we ca discover that the Gamma fuctio has simple poles i all egative itegers simply iteratig the recursio (2): Γ( ) Γ( + )... ( ) Γ() ( N). () This brigs to the coclusio that Γ( ) has a simple pole of order ad residue ( ). More rigorously, this ca be show i the usual way: Res Γ( ) Γ( )... Γ( + ) 2

( )( + ) ( ) Γ() ( ) ( N). (2) I all the complex plae (except the egative real axis) the Gamma fuctio is well defied. Fo o-iteger egative real values the Gamma fuctio ca be aalytically cotiued (as we have see for example for half-itegers, positives ad egatives). We ca prove these results i a simpler way. Startig from the defiitio () we rewrite Γ() dt e t t + dt e t t dt e t t + dt e t t + dt t ( ) ( ) + where the first itegral is a aalytic fuctio of ad the secod term gives the positio of the poles ad their residues. Moreover, the Gamma fuctio has a essetial sigularity to complex ifiity, because Γ ( ) has a o-defied limit for. This is to say that the Gamma fuctio is ot well defied i the compactified complex plae. Now we wat to show that ear a pole i oe has the expasio Γ( ) ( ) [ + ψ( + ) + O() ] where ψ() is the logarithmic derivative of the Gamma fuctio, defied i (4); ote that ψ( + ) has a explicit expressio give i eq. (44). First, ote that where we have defied Γ( ) Γ( + )... Γ( + + ) ( )( + ) ( ) ( + ) ( + ) G(, ) t (3) (4) G(, ) (5) Γ( + + ) ( ) [ ( 2 2 ) ( 2 Γ( + + ) + O( 2 ) () 2 ) ]. (6) By Taylor expadig G(, ) aroud i (5) we get that is the claimed result. Γ( ) ( ) [ G(, ) + G (, ) + O( 2 ) ] [ + Γ ( + ) + O()] 3, (7)

.2 Alterative defiitios There are some alterative defiitios of the Gamma fuctio. Oe, due to Euler, is Γ() ( + ) + ( + ) ( + ) ( + ). (8) It is easy to prove (usig the expressio with the limit) that with this defiitio the recursio relatio (2) is satisfied, ad also that for gives Γ(). Aother oe, due to Weierstrass, is Γ() e γ e / + e γ lim [ ( exp + 2 +... + ( + ) ( + ) )], (9) where γ is the Euler-Mascheroi umber defied i Sectio 2.2 The equivalece to the Euler defiitio ca be see by takig the logarithm of this two forms: the first gives ad the secod gives log Γ() log + log Γ() log γ + [ ( log + ) ( log + )] [ ( log + )] (2). (2) Comparig this two expressios we ote that they are equal because of eq. (48). Note that this two expressio are useful i order to umerically approximate log Γ(), by trucatig the series to a certai value of..3 Euler reflectio formula A useful formula is Γ() Γ( ) π si(π) called the Euler reflectio formula. To demostrate it, cosider the alterative defiitio (8) ad ote that where we have used Γ() Γ( ) Γ() Γ( ) ( ) ( + ) ( ) ( + ) ( + ) ( si(π) π 2 2 ) (22) (23) ( ) si(π) π 2 2. (24) 4

.4 Stirlig approximatio The Stirlig approximatio for the factorial is 2π e (25) which is valid for large. To demostrate it, we perform a saddle poit approximatio to the itegral (): Γ( + ) dt e t t t+ log t dt e ow we expad the expoet ear its maximum ad get that is the claimed result..5 Derivative recursio relatio dt exp ( + log ) 2 t2 e dt e 2 t2 Cosider the Taylor expasio ear of Γ( + + ): 2π e (26) Γ( + + ) k Γ (k) ( + ) k ; (27) usig eq. (2) this is equal to ( + ) Γ( + ), that expaded is ( + ) Γ( + ) ( + ) k Γ( ) + Γ (k) ( ) Comparig the two expasios, we get the equality (k > ) k k k [ Γ (k) ( ) + k Γ (k ) ( ) ]. (28) Γ (k) ( + ) Γ (k) ( ) + k Γ (k ) ( ). (29).6 Iverse Gamma fuctio We ca defie a iverse Gamma fuctio () Γ() (3) 5

2..5..5 2 2 4 6.5. Figure 2: Re (x) for real x. which is a etire fuctio (aalytical i all the complex plae). The recursio relatio (2) becomes () ( + ). (3) We are able to produce a itegral represetatio of the fuctio: cosider the Laplace trasform of the fuctio t with respect to t: dt e st t s Γ() (32) ad the we have that the iverse Laplace traform of s is c+i ds e ts s () t (33) 2πi c i where c has to be choose to the right of all the sigularities of s, that is c >. Now, by choosig t, we get () c+i ds e s s, (34) 2πi c i where the itegratio path ca be deformed to be a cotour of the egative real axis. Note that for positive itegers + the brach cut vaishes ad oly a pole i s remais, ad the the path ca be closed aroud s to obtai 2πi that ca be also easly demostrated usig the residue theorem. Now cosider the Taylor expasio ear of ( + ): ( + ) k ds e s s, (35) (k) ( ) usig eq. (3) this is equal to ( + ) ( + + ), that expaded is ( + ) ( + + ) ( + ) k (k) ( + ) 6 k k ; (36)

k [ ( + ) + (k) ( + ) + k (k ) ( + ) ]. (37) Comparig the two expasios, we get the equality (k ) k (k) ( ) (k) ( + ) + k (k ) ( + ) (38) ad, for, (k) () k (k ) (). (39).6. Saddle poit approximatio The aalogous of the Stirlig approximatio ca be doe for the () fuctio. The result is: () 2π e. (4) 7

2 Gamma related fuctios 2. Log derivative ψ() Defie the logarithmic derivative of the Gamma fuctio as follows Form the recursio relatio (2) it follows that ψ() d log Γ(). (4) d ψ( + ) ψ() +, (42) as we ca easly demostrate: ψ( + ) d log Γ( + ) d d log [ Γ()] d + ψ() (43) (ote that the same result ca be obtaied from eq. (29) for k ). For a positive iteger, we ca iterate the recursio (42) to obtai ψ( + ) ψ() + + 2 + 3 +... +, (44) ad the value of ψ() is ψ() γ (45) where γ is the Euler-Mascheroi umber, as demostrated i Sectio 2.2. The fuctios ψ () ψ () () d+ log Γ() (46) d+ are also defied. Iteratig the derivative i (42) we obtai the recursio relatio 2.2 Euler-Mascheroi umber The Euler-Mascheroi umber is defied as ψ ( + ) ψ () + ( ) +. (47) γ [ ( log + )], (48) ad its value is γ.57726.... (49) 8

Note that ( log + ) N N N ( log + ) N [ ] log( + ) log N log( + N) N log N, (5) where we oted that the sum is a telescopic sum ad that i the limit N log( + N) ad log N are equivalet. Thus we ca rewrite eq. (48) as γ N N log N. (5) Now we prove that ψ() γ. To do this, let us cosider eq. (2), ad take its derivative: ψ() [ ( + log + ) ] +. (52) For we get ψ() + i which we recogise the defiitio (48). 2.3 Geeralied Gamma fuctios [ ( log + ) ] + [ ( log + )] You ca also defie a icomplete Gamma fuctio (or plica fuctio) ad a trucated Gamma fuctio Γ(, a) γ(, a) a a, (53) dt e t t (54) dt e t t, (55) that have a brach cut o the egative real axis i the a complex plae. Obviously Γ() Γ(, a) + γ(, a). (56) 9

For iteger k +, itegratig repeatedly by parts we get 2.4 Beta fuctio The Beta fuctio is defied as Γ(k +, a) e a k γ(k +, a) e a a (57) k a. (58) B(a, b) dx x a ( x) b Γ(a)Γ(b) Γ(a + b). (59) To show the last equality, cosider the product Γ(a)Γ(b) 4 dt e t t a du e u u b dx 2π 2 Γ(a + b) Γ(a + b) dx e x2 x 2a dy e y2 y 2b dy e (x2 +y 2) x 2a y 2b dθ cos 2a θ si 2b θ π/2 dθ cos 2a θ si 2b θ dr e r2 r 2a+2b dx x a ( x) b. (6)