+ " calculated: " + calculated + "\n"
+ " expected: " + expected);
for (int i=2; i<9000; i++) {
expected += Math.log(i);
calculated = gfn.logValue(i+1);
if (Math.abs(calculated-expected)/expected > 1.0e-9)
fail("GammaFunction(" + i + ")\n"
+ " calculated: " + calculated + "\n"
+ " expected: " + expected + "\n"
+ " rel. error: " + Math.abs(calculated-expected)/expected);