scipy.special.beta
scipy.special.betainc
y=betaln(a,b) returns the natural logarithm of the absolute value of beta: ln(abs(beta(x))).