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