[X,Y,Z] = peaks;
contourf(X,Y,Z,12) % filled contour plot with 12 contours
xlabel('X-axis'), ylabel('Y-axis')
title('Figure 26.23: Filled Contour Plot')