[X,Y,Z] = peaks(30);
surf(X,Y,Z)  % same plot as above
shading flat
xlabel('X-axis'), ylabel('Y-axis'), zlabel('Z-axis')
title('Figure 26.10: Surface Plot with Flat Shading')