SciPy

scipy.special.beta

scipy.special.beta(a, b) returns gamma(a) * gamma(b) / gamma(a+b) = <ufunc 'beta'>