Beta function.
beta(a,b) = gamma(a) * gamma(b) / gamma(a+b)
scipy.special.gammainccinv
scipy.special.betaln