this.q = new QFunction(new NFunction(alpha1, beta1),
new NFunction(alpha2, beta2), p);
}
public double eval(int n, double e) {
GammaDistribution g1 = new GammaDistributionImpl(alpha1 + n, beta1 + e);
GammaDistribution g2 = new GammaDistributionImpl(alpha2 + n, beta2 + e);
PiFunction pi = new PiFunction(q.eval(n, e), g1, g2);
PiFunctionIntegral ipi = new PiFunctionIntegral(pi, target);