N! AND THE GAMMA FUNCTION

Similar documents
SUMMATION OF INFINITE SERIES REVISITED

Calculus Limits. Limit of a function.. 1. One-Sided Limits...1. Infinite limits 2. Vertical Asymptotes...3. Calculating Limits Using the Limit Laws.

1. Solve by the method of undetermined coefficients and by the method of variation of parameters. (4)

12 Getting Started With Fourier Analysis

Moment Generating Function

K3 p K2 p Kp 0 p 2 p 3 p

An interesting result about subset sums. Nitu Kitchloo. Lior Pachter. November 27, Abstract

Extremal graph theory II: K t and K t,t

th m m m m central moment : E[( X X) ] ( X X) ( x X) f ( x)

Ideal Amplifier/Attenuator. Memoryless. where k is some real constant. Integrator. System with memory

ECE-314 Fall 2012 Review Questions

Solutions to Problems 3, Level 4

CLOSED FORM EVALUATION OF RESTRICTED SUMS CONTAINING SQUARES OF FIBONOMIAL COEFFICIENTS

King Fahd University of Petroleum & Minerals Computer Engineering g Dept

UNIT 1: ANALYTICAL METHODS FOR ENGINEERS

Using Linnik's Identity to Approximate the Prime Counting Function with the Logarithmic Integral

ECE 350 Matlab-Based Project #3

Lecture 9: Polynomial Approximations

Extended Laguerre Polynomials

L-functions and Class Numbers

Math 6710, Fall 2016 Final Exam Solutions

S n. = n. Sum of first n terms of an A. P is

Comparison between Fourier and Corrected Fourier Series Methods

The Central Limit Theorem

A TAUBERIAN THEOREM FOR THE WEIGHTED MEAN METHOD OF SUMMABILITY

Review Exercises for Chapter 9

Section 8 Convolution and Deconvolution

Math 2414 Homework Set 7 Solutions 10 Points

Department of Mathematical and Statistical Sciences University of Alberta

Lecture 15 First Properties of the Brownian Motion

Notes 03 largely plagiarized by %khc

Sampling Example. ( ) δ ( f 1) (1/2)cos(12πt), T 0 = 1

BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS

The Eigen Function of Linear Systems

Big O Notation for Time Complexity of Algorithms

2 f(x) dx = 1, 0. 2f(x 1) dx d) 1 4t t6 t. t 2 dt i)

David Randall. ( )e ikx. k = u x,t. u( x,t)e ikx dx L. x L /2. Recall that the proof of (1) and (2) involves use of the orthogonality condition.

6.003: Signals and Systems Lecture 20 April 22, 2010

Research Article On a Class of q-bernoulli, q-euler, and q-genocchi Polynomials

Fermat Numbers in Multinomial Coefficients

LIMITS OF FUNCTIONS (I)

LIMITS OF SEQUENCES AND FUNCTIONS

Hadamard matrices from the Multiplication Table of the Finite Fields

Solutions to selected problems from the midterm exam Math 222 Winter 2015

In this section we will study periodic signals in terms of their frequency f t is said to be periodic if (4.1)

1 Notes on Little s Law (l = λw)

A Generalization of Hermite Polynomials

Some identities related to reciprocal functions

ODEs II, Supplement to Lectures 6 & 7: The Jordan Normal Form: Solving Autonomous, Homogeneous Linear Systems. April 2, 2003

Problems and Solutions for Section 3.2 (3.15 through 3.25)

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 4 9/16/2013. Applications of the large deviation technique

Mathematical Statistics. 1 Introduction to the materials to be covered in this course

Lecture 15: Three-tank Mixing and Lead Poisoning

Exercise 3 Stochastic Models of Manufacturing Systems 4T400, 6 May

Energy Density / Energy Flux / Total Energy in 1D. Key Mathematics: density, flux, and the continuity equation.

A Note on Random k-sat for Moderately Growing k

The Connection between the Basel Problem and a Special Integral

BE.430 Tutorial: Linear Operator Theory and Eigenfunction Expansion

6.003: Signal Processing Lecture 2a September 11, 2018

Actuarial Society of India

On Another Type of Transform Called Rangaig Transform

F.Y. Diploma : Sem. II [AE/CH/FG/ME/PT/PG] Applied Mathematics

The analysis of the method on the one variable function s limit Ke Wu

SOLUTION SET VI FOR FALL [(n + 2)(n + 1)a n+2 a n 1 ]x n = 0,

On The Eneström-Kakeya Theorem

ANALYSIS OF THE CHAOS DYNAMICS IN (X n,x n+1) PLANE

STK4080/9080 Survival and event history analysis

A Study On (H, 1)(E, q) Product Summability Of Fourier Series And Its Conjugate Series

Enumeration of Sequences Constrained by the Ratio of Consecutive Parts

(C) x 3 + y 3 = 2(x 2 + y 2 ) (D) x y 2 (A) (10)! (B) (11)! (C) (10)! + 1 (D) (11)! 1. n in the expansion of 2 (A) 15 (B) 45 (C) 55 (D) 56

International journal of Engineering Research-Online A Peer Reviewed International Journal Articles available online

arxiv: v1 [math.nt] 13 Dec 2010

Dynamic h-index: the Hirsch index in function of time

MATH 507a ASSIGNMENT 4 SOLUTIONS FALL 2018 Prof. Alexander. g (x) dx = g(b) g(0) = g(b),

A note on deviation inequalities on {0, 1} n. by Julio Bernués*

Electrical Engineering Department Network Lab.

Zhi-Wei Sun and Hao Pan (Nanjing)

F D D D D F. smoothed value of the data including Y t the most recent data.

Applying the Moment Generating Functions to the Study of Probability Distributions

CSE 241 Algorithms and Data Structures 10/14/2015. Skip Lists

AN EXTENSION OF LUCAS THEOREM. Hong Hu and Zhi-Wei Sun. (Communicated by David E. Rohrlich)

Calculus BC 2015 Scoring Guidelines

ON THE n-th ELEMENT OF A SET OF POSITIVE INTEGERS

Outline. simplest HMM (1) simple HMMs? simplest HMM (2) Parameter estimation for discrete hidden Markov models

Linear Time Invariant Systems

INVESTMENT PROJECT EFFICIENCY EVALUATION

CHAPTER 2 Quadratic diophantine equations with two unknowns

λiv Av = 0 or ( λi Av ) = 0. In order for a vector v to be an eigenvector, it must be in the kernel of λi

Chapter 2: Time-Domain Representations of Linear Time-Invariant Systems. Chih-Wei Liu

Solution. 1 Solutions of Homework 6. Sangchul Lee. April 28, Problem 1.1 [Dur10, Exercise ]

Samuel Sindayigaya 1, Nyongesa L. Kennedy 2, Adu A.M. Wasike 3

Closed-Form Solution for the Nontrivial Zeros of the Riemann Zeta Function

Supplement for SADAGRAD: Strongly Adaptive Stochastic Gradient Methods"

Stochastic Processes Adopted From p Chapter 9 Probability, Random Variables and Stochastic Processes, 4th Edition A. Papoulis and S.

ES 330 Electronics II Homework 03 (Fall 2017 Due Wednesday, September 20, 2017)

Some inequalities for q-polygamma function and ζ q -Riemann zeta functions

BINOMIAL COEFFICIENT AND THE GAUSSIAN

GAUSSIAN CHAOS AND SAMPLE PATH PROPERTIES OF ADDITIVE FUNCTIONALS OF SYMMETRIC MARKOV PROCESSES

6.003 Homework #5 Solutions

Online Supplement to Reactive Tabu Search in a Team-Learning Problem

Transcription:

N! AND THE GAMMA FUNCTION Cosider he produc of he firs posiive iegers- 3 4 5 6 (-) =! Oe calls his produc he facorial ad has ha produc of he firs five iegers equals 5!=0. Direcly relaed o he discree! fucio oe has he coiuous gamma fucio Γ(). Boh are defied by he same iegral Oe iegraio by pars yields- ( )! 0 exp( ) d from which follow he ideiies- 0 0 exp( ) d exp( ) d (+)!=!(+) ad Γ(+)=Γ() We also have ha 0!=!= ad Γ(0)=, Γ(/)=sqr(π) ad Γ()=. A plo of! ad Γ(+) follow-

The blue dos are! values while he red curve represes he coiuous Γ(+) fucio. The curve reaches a miimum value of Γ()= 0.88560 a =.46. Oe ypically fids he values for Γ() are abulaed oly i he rage <<, sice he res ca be quickly geeraed via he above recurrece formula. Alhough here is o rue facorial for egaive iegers, oe ca exed he Γ() fucio o egaive ad obai o ifiie values whe is a egaive o-ieger. We ca also express he derivaive of he gamma fucio as he iegral- d( ) d 0 exp( ) d ( ) d d 0 l( ) exp( ) d This derivaive goes oward - as ->0, has zero value ear =.46 ad akes o progressively larger posiive values as heads oward plus ifiiy. Noice ha his las iegral is jus he Laplace rasform of (-) l() afer s is se equal o uiy. A fucio relaed o Γ(x) is he digamma fucio- ( x) d( x) / dx ( x) d[l ( x)] dx I has he value ψ()=-γ where γ=0.57756649.. is he Euler cosa. Oe ca also sum he reciprocals of various combiaios of!. We have, amog may oher examples, ha- exp() I () 0 0 cosh() sih() 0! (!) 0 0.78888459045...795853033607...543080634854.. ()!.7509364380.. ( ) Noice also ha- ()!=[ 3 5 (-)][ 4 6 ]= [ 3 5 (-)]! Thus we have ha- 3 5 (-)=()!/(!)

From his i follows ha he produc of he firs five odd umbers equals 0!/(3*5!)=945. This las form also allows oe o wrie cerai ifiie series i compac form. For example, we have ha- x x ( 3) x ( 35) x... (m)! 4 6 3 m!! 3! m 0 ( m! ) As we leared i our earlier discussios o Legedre polyomials P (x), hese ca be geeraed by he geeraig fucio - P ( x) x 0 So o seig ε=x- we ca wrie- x m + 0 P ( x) (m)! m! [ (x )] m m0 ( m!) (!) 3 x [(3x ) / ] O( ) which produces he Legedre polyomials. 4! [ (x )] 4 (!) Oe ca also use he gamma fucio Γ() o evaluae 3 5 -. We have Γ(/)=sqr(π) so ha Γ(3/)=sqr(π)/, Γ(5/)=( 3)sqr(π)/ ad Γ(7/)=( 3 5)sqr(π)/ 3. From his i follows ha- [ 3 5 (-)]= Γ(+/)/sqr(π) [ (x )] Combiig his resul wih he form of ()! give earlier, oe obais he Legedre Duplicaio Formula-! ()! ( ) Tryig his ou for =5, we fid- 0!=( 0 5!)Γ(/)/sqr(π)=368800... Aoher combiaio of facorials which ofe arises is he famous biomial coefficie- C m =!/[(m!(-m)!] I is produced by he followig biomial expasio-

( a b) a a b /! ( ) a b /! m0 C m a m b m Noe ha he C m for a fixed jus represes he umbers i he h row of a Pascal riagle. Thus he 4 h row has he coefficies C 4m =4!/[(m!(4-m)!] which are -4-6-4-. Aoher exesio of he facorial is he produc of squares which read- F()= 4 9 6 This is easy o evaluae by oig F() is jus he produc of! wih iself. Tha is- F()=( 3 )( 3 )=(!) I also follows ha he produc of he firs ph powers of he iegers equals (!) p. Thus 8 7 64 5=(5!) 3 =78000 Nex we examie he value of Γ(+/). Usig he Legedre Duplicaio Formula ad he form for ()! give earlier, we have- ()! ( )! This allows oe o fid he half-ieger gamma fucio. I says Γ(3/)=30!sqr(π)/( 30 5!)= (6908335369375/3768)sqr(π) As expeced his value lies bewee 4! Ad 5!. I is also possible o develop gamma fucio ideiies o foud i exisig mahemaical hadbooks. Oe of hese is- G()=Γ(+/) Γ(-/) We develop he geeral value for G() by sarig wih = where G()=π/. Nex a = we have G()=+3π/8 ad a =3 we fid G(3)=45π/3. This suggess ha ( ) ( ) G( ) { 3 5.. ( 3) } [ ( )] This ideiy checks for all values of ried for of oe or greaer.

Aoher variaio is he gamma produc fucio- P(x,y)=Γ(x+y) Γ(x-y) which reduces o (x+y-)! (x-y-)! whe x ad y are iegers. A coour plo of his fucio for x>0 ad -4<y<4 looks like his- The coours form closed curves ad P(x,y) goes o ifiiy whe y=±(x+) sice gamma for ay egaive ieger i ubouded. Alerae srip regios bewee he ui slope curves also show fiie valued coours. The miimum coour value occurs ear x=.46 ad y=0 ad has he value P=0.7844. Fially we look a he gamma fucio Γ(z) whe z=x+iy is a complex umber. Here he bes approach is o use he iegral defiiio- ( x iy) 0 xiy exp( ) d 0 x exp( ) cos( y l ) i si( y l ) d From he iegral we see ha Γ(x+iy) has a real ad imagiary par represeed by wo differe iegrals. We fid- Γ(+i)= 0.49805668-i0.54949883 so ha Γ(+i) Γ(-i)= Γ(+i) =0.7090550... Also we have ha-

0 cos[l( )]exp( )) d 0.49805668.. We ca also plo Γ(x+iy)=u+iv i he u-v plae. This ca produce some ieresig figures such as he followig- I he firs we plo Γ(z) for z=x+i o ge a ru-away spiral paer. For he secod figure we have se z=+iy. I produces a closed double loop. May oher plos are possible by jus seig x or y o differe cosa values. March 03