pts[k][1]=(double)(i);
y[k]=cos(pts[k][0]/20.0)*cos(pts[k][1]/20.0);
}
}
Powvargram vgram = new Powvargram(pts,y);
Krig krig = new Krig(pts,y,vgram);
for (i=0;i<M;i++) {
ppt[0]=pt[i][0];
ppt[1]=pt[i][1];
estim[i]=krig.interp(ppt);
}
System.out.printf(" Discrepancy: %f\n", maxel(vecsub(actual,estim)));
localflag = maxel(vecsub(actual,estim)) > sbeps;
globalflag = globalflag || localflag;
if (localflag) {