algdiv function

Compute log(gamma(b)/gamma(a+b)) when b >= 8