m=0.5;s=1.5;
Normaldist normc = new Normaldist(m,s);
Ran myran = new Ran(17);
sbeps=5.0e-14;
for (i=0;i<1000;i++) {
u=m-3.0*s+6.0*s*myran.doub();
a=normc.cdf(u);
b=normc.invcdf(a);
// System.out.printf(setprecision(15) << u << " %f\n", b << " %f\n", abs(u-b));
localflag = localflag || abs(u-b) > sbeps;
}