public double evaluate(double gam) { if (gam <= 0) return 1.0e100; double sum = 0.0; for (int j = 0; j < m; j++) sum += Num.digamma(gam + x[j]); return sum / m + Math.log(p) - Num.digamma(gam); }