public static double Log(double a, double b) { return Gamma.Log(a) + Gamma.Log(b) - Gamma.Log(a + b); }