Following from last posting [baiscs](https://steemit.com/math/@beoped/math-computation-the-expansion-of-gamma-function-1-basic), at this time we are going to derive 

![](https://steemitimages.com/DQma9VR4N2R9X8FYXyMLEYfaVmuSD6xGhfvSB5jZfbErRSf/image.png)
this equation and do some application. 

-------------------------------------------------------------



First from recursive relation of gamma function, i.e., Gamma(z+1) = z Gamma(z), we have 

![](https://steemitimages.com/DQmUSs9cSP4XFWTcPKbzBxgReXskjHMCtQwgSA8SRc96EVW/image.png)

This is a usual textbook approach for some negative factorials. Using the formula in the previous post
![](https://steemitimages.com/DQmYXJmTq658cbYsUxuHkQJ7tMNQtRJesN7XcWUL4Ntvcw9/image.png)
with the gemoetric series, ![](https://steemitimages.com/DQmPQcyx9QuSqWncPtsPmhQjgNrdbidTv5GreoA1zhg1Sg1/image.png) by collecting terms of epsilon^2 in the denomiantor, we can easily obtain above equation. 


The next thing to do is taylor expansion about Gamma(n+1+epsilon) about the parameter epsilon. 

![](https://steemitimages.com/DQmUrVwe7TiLG2itQNvsyGXb9i9AQ77Qd3VATPM5KvaYsvT/image.png)

In the process we used
![](https://steemitimages.com/DQmQF1s1qDFwHHxLKWJdcDXkgJw5YxXvSsTgAkfQq2cDF4s/image.png)
which we derived in the previous post. 


Now just plugging these results we have

![](https://steemitimages.com/DQmd3DeezJk8tad7aAEoS36KJaknF4ZkHxrNNDB1WLmZwTf/image.png)
Note that upon here we computed upto order 4 of epsilon, but you can see that there is another epsilon in the denoimnator so in this computation we obtain at most order 3 of epsilon. To compute more higer terms we need to compute more higher orders. 

Now we are left with plugging n=0 and re-ordering. Using the values of psi(1), psi'(1), psi''(1) which we derive in previous post
![](https://steemitimages.com/DQmP6bAP1A2vwo2kwr5aS7xbyrjFgtoKyH8nrrLKZj6kcQv/image.png)

by plugging these, we have

![](https://steemitimages.com/DQma9VR4N2R9X8FYXyMLEYfaVmuSD6xGhfvSB5jZfbErRSf/image.png) 

Which is desired results. 


Application
--------------------

With the combination of complex integration we do more than that. 

First state some formula about complex integration

![](https://steemitimages.com/DQmYdWXtgYvfjinnd565hAZJjkDssTptTySVQsZ1bJgAagV/image.png)

From 
![](https://steemitimages.com/DQma9VR4N2R9X8FYXyMLEYfaVmuSD6xGhfvSB5jZfbErRSf/image.png)  we have

![](https://steemitimages.com/DQmRu7xaoKc16Fh3dBUqxJ6nN1YYYMG2Ep5xDFGjpBF4Ax3/image.png)
 

In case of no pole, contour integral vanishes, so for arbitrary function A(z), we can do 

![](https://steemitimages.com/DQmc8Y973LFsRZLk7fuUJF7J9KDDQmh6RovHJJKTeUDoJsd/image.png)

In this way we can use gamma fucntion and do the integration more easily. For another type of special function, like beta function or generalized gamma function, etc., you can do the similar things!.