x=linspace(-.5,1.5);
y=humps(x);
plot(x,y)
grid on
title('Figure 22.1: Humps Function')
 
%(mm2201.m plot)