localflag=false;
for (i=0;i<N;i++) {
Midpnt qq2 = new Midpnt(dist2,0.0,x[i]);
integral=qromo(qq2);
c[i]=integral;
d[i]=normcdf.cdf(x[i]);
// System.out.printf(c[i]-d[i]);
localflag = localflag || abs(c[i]-d[i]) > sbeps;
}
globalflag = globalflag || localflag;
if (localflag) {