% hump2.m
function hump2
  diary('hump2_m.out')
  format compact
  disp('hump2.m running');
  a=[0:0.1:1.0]; 
  b=[0:0.1:1.0];
  x=zeros(11,11);
  y=zeros(11,11);
  z=zeros(11,11);
  for i=1:11
    for j=1:11
      x(i,j) = a(i);
      y(i,j) = b(i);
      z(i,j)=16*a(i)*(1-a(i))*b(j)*(1-b(j));
    end
  end
  disp('x=');
  disp(x);
  disp('y=');
  disp(y);
  disp('z=');
  disp(z);
  diary off
  figure
  mesh(z)
  title('hump2.m')
  xlabel('x');
  ylabel('y');
  zlabel('z');

end

