It s a simply gamma function for multiply input 

To calculate of Pe(a, b) = gamma(a+0.5)*gamma(b+0.5)/(pi*gamma(a+b+1)), which in this case. The output requires log.