SciPy

scipy.special.beta

scipy.special.beta(a, b) = <ufunc 'beta'>

Beta function.

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