localflag=false;
sbeps=5.e-10;
for (i=0;i<1000;i++) {
z=0.3+3.0*myran.doub();
a=norm.qks(z);
b=norm.invqks(a);
// if (abs(z-b) > sbeps) {
// System.out.printf(setprecision(15) << z << " %f\n", b << " %f\n", abs(z-b));
// }
localflag = localflag || abs(z-b) > sbeps;
}