Semi-colon % prevents output of command.
% Define a collection of x's & y's.
x=0:.05:1;y=0:.05:1;
% Define a mesh of points.
[xx,yy]=meshgrid(x,y);
% Compute the T distribution at each point
zz=some function of xx and yy
% Plot surface
surf(zz);
view(120.,30.);
%shading interp;
%colormap(jet);
pause
% Compute gradient field
[px,py] = gradient(zz,.05,.05);
% contours
% Plot colormap(jet);
contour(x,y,zz,20);
% keep contours on the screen & plot vectors
hold on;
quiver(x,y,px,py,2);
hold off;
pause
% Get gradient at (.5,.5)==> point (11,11)
gTp(1)=px(11);...