CrigFunct f1 = new CrigFunct();
f1.setA(a);
f1.setB(Stat.logGammaFunction(a));
// Instantiate Integration
Integration intgn1 = new Integration(f1);
double xx = x;
double yy = x + incr;
intgn1.setLimits(xx, yy);
// Perform quadrature
sum = intgn1.gaussQuad(64);
boolean test2 = true;
for(int i=1; i<nIncr; i++){
xx = yy;
yy = xx + incr;
intgn1.setLimits(xx, yy);
sum += intgn1.gaussQuad(64);
}
return sum;
}