sbeps=5.0e-14;
localflag=false;
for (i=0;i<1000;i++) {
u=myran.doub();
a=normc.cdf(u);
b=normc.invcdf(a);
// if (abs(u-b) > sbeps) {
// System.out.printf(setprecision(15) << u << " %f\n", b << " %f\n", abs(u-b));
// }
localflag = localflag || abs(u-b) > sbeps;
}