x=0:0.1:6;
y=0:0.1:10;
[X,Y]=meshgrid(x,y);
figure(1);
hold on
T=log(X+1).*(exp(-(Y-2).^2));
pcolor(X,Y,T);
shading flat
grid on
axis([0,6,0,10]);
colorbar;
title('Campo de temperaturas')
hold off