
x = rand(100, 1);
y = rand(100, 1);
z = x.^2 + y.^2;
tri = delaunay(x, y);

figure(2); clf; hold on;
trisurf(tri, x, y, z);
shading interp; axis equal; box on; colorbar;