x=0:0.1:2;
y=0:0.1:10;
[X,Y]=meshgrid(x,y);
rho=sqrt(X.^2+Y.^2);
hold on
T=log(rho+1).*(exp(-(Y-2).^2));
surf(X,Y,T);
view
axis([0,3,0,10]);
colorbar
title('Campo de temperaturas')
hold off