function contour(x,y,z)
 %
 xlin=linspace(min(x),max(x),33);
 ylin=linspace(min(y),max(y),33);
 [X,Y]=meshgrid(xlin,ylin);
 Z=griddata(x,y,z,X,Y,'cubic');
 %
 contour(X,Y,Z)