final BinomialDistribution binomdist = new BinomialDistribution(p,n);
for(int i=0;i<n;i++){
final int z = i;
final double expected = testvalues[i];
final double computed = binomdist.op(z);
// double tolerance = (z<6 ) ? 1.0e-15: 1.0e-10;
final double tolerance = 1.0e-15; // try to to get 1.0e-15 accuracy whenever possible
//assertEquals(expected, computed, tolerance);