for i=1:41,
   for j=1:41,
      x(i) = (i-1)/10 - 2;
      y(j) = (j-1)/10 - 2;
      z(j,i) = exp(x(i)).*(4*x(i)*x(i)+ ...
         2*y(j)*y(j)+4*x(i)*y(j)+2*y(j)+1);
   end;
end;
contour(x,y,z,[.2 .4 .7 1 1.7 1.8 2 3 4 5 6 20]);
grid on;