[X,Y,Z] = peaks(15);
surfnorm(X,Y,Z)
xlabel('X-axis'), ylabel('Y-axis'), zlabel('Z-axis')
title('Figure 26.15: Surface Plot with Normals')