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);
}