[nr,nc]=size(xi);
## triangulate data
tri=delaunay(x,y);
tri=delaunay(x,y,"QJ");
zi=zeros(size(xi));
if strcmp(method,'cubic')