scipy.special.beta

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

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