{14.0, 7.665444025768806E-9},
{15.0, 2.659774281674882E-9}};
// Test for a = 1.0 (alpha)
double a = 1.0;
final GammaDistribution gammDistribution = new GammaDistribution(a);
for (final double[] testvalue : testvalues) {
final double expected = testvalue[1];
final double x = testvalue[0];
final double computed = gammDistribution.op(x);
// QL.info(computed); // for testing
final double tolerance = 1.0e-15;
if (Math.abs(expected-computed)>tolerance) {
fail("x: " + x + " expected: " + expected + " realised: " + computed);
}
}
// Test for a = 0.1 (alpha)
a = 0.1;
final GammaDistribution gammDist = new GammaDistribution(a);
for (final double[] element : testvalues2) {
final double expected = element[1];
final double x = element[0];
final double computed = gammDist.op(x);
// QL.info(computed); // for testing
final double tolerance = 1.0e-15;
if (Math.abs(expected-computed)>tolerance) {
fail("x: " + x + " expected: " + expected + " realised: " + computed);
}