x=-0.5:0.1:0.5;       % Creamos los intervalos x= [-1/2,1/2] y=[0,2]
y=0:0.4:4;             
[xx,yy]=meshgrid(x,y); %Creamos la matriz a partir del intervalo 
mesh(xx,yy,0*xx) % Dibujamos malla
hold on
axis ([-0.5 0.5, 0 4, 0 10])      % seleccionamos región
% temperatura
T=exp(-xx.^2-2*xx-1).*(8-yy.^2+2*yy); 
surf(xx,yy,T)          % Dibujamos la superficie 
axis  ([-0.5 0.5, 0 4, 0 10])     % Seleccionamos la región
%Como ya tenemos definida la malla y el campo T, calculamos directamente el gradiente
Tx=-2*(xx+1)*(8-yy.^2+2*yy).*exp(-xx.^2-2*xx-1);
Ty=exp(-xx.^2-2*xx-1).*(-2*yy+2);
%Dibujamos el Vector gradiente (figura 3)
quiver(xx,yy,Tx,Ty)
contour(xx,yy,T)       % lineas de nivel
 
plot(x,x-x,'m','linewidth',1);
plot(x,2+x-x,'m','linewidth',1);
plot(-0.5+y-y,y,'m','linewidth',1);
plot(0.5+y-y,y,'m','linewidth',1);
axis([5,5,5]) %Numeramos los ejes
hold off