There are several options, likely the easiest being:
xv = linspace(min(x), max(x), numel(x));
yv = linspace(min(y), max(y), numel(y));
[X,Y] = ndgrid(xv,yv);
Z = griddata(x, y, z, X, Y);
That will create a uniform grid in ‘Z’ that interpolates the original ‘z’ values and can then be plotted as desired.
