m=0.5;s=1.5;
func_Normaldist dist2 =new func_Normaldist(m,s);
Normaldist normcdf = new Normaldist(m,s);
for (i=0;i<N;i++) {
c[i]=qsimp(dist2,-20.0,x[i]);
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) {