Binomialdist normc=new Binomialdist(n,p);
sbeps=5.0e-14;
localflag=false;
for (i=1;i<10;i++) {
a=normc.cdf(i);
jb=normc.invcdf(a);
// if (abs(i-jb) > sbeps) {
// System.out.printf(setprecision(15) << i << " %f\n", jb << " %f\n", abs(i-jb));
// }
localflag = localflag || (i-jb) != 0 && (i-jb-1) != 0;
}