pts[k][0]=(double)(j);
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);