gamma function 

KT estimation: Pe(a, b) = gamma(a+0.5)*gamma(b+0.5)/(pi*gamma(a+b+1))