scipy.special.betaln

scipy.special.betaln(x1, x2[, out])

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