Logarithm of the beta function based on the gammaLn function.
double betaLn(num x, num y) => x <= 0 || y <= 0 ? double.nan : gammaLn(x) + gammaLn(y) - gammaLn(x + y);