function hump
  % hump.m
  diary('hump_m.out')
  format compact
  disp('hump.m running');
  a=[0:0.1:1.0];
  b=[0:0.1:1.0];
  c=zeros(11,11);   % like Z above
  for i=1:11
    for j=1:11
      c(i,j)=16*a(i)*(1-a(i))*b(j)*(1-b(j));
    end
  end
  disp('a=');
  disp(a);
  disp('b=');
  disp(b);
  disp('c=');
  disp(c);
  diary off
  figure;
  mesh(c);
  title('hump.m')
  xlabel('a');
  ylabel('b');
  zlabel('c');
  
end % hump
