scipy.special.beta

scipy.special.beta(x1, x2[, out]) = <ufunc 'beta'>

beta(a, b)

Beta function.

beta(a,b) =  gamma(a) * gamma(b) / gamma(a+b)