function z = Humps2D(x,y);
   z = 1/( (.2*(x-3)^2 +.1) + (.3*(y-1)^2 + .1) );
