\[\gamma\left(a,z\right)\sim-z^{a}e^{-z}\sum_{k=0}^{\infty}\frac{(-a)^{k}b_{k}(%
\lambda)}{(z-a)^{2k+1}},\]