function plothillshade(Easting, Northing, Hillshds)    
%Given previously computed Easting and Northing matrices (from LoadFile function)
%and the Hillshds hillshade matrix (from hillshademe function), plot a
%hillshade map.
%
%Ramon Arrowsmith, 2009

plot(Easting(1,1),Northing(1,1), 'b.');
hold on
imagesc(Easting,Northing,Hillshds);
colormap('bone'); 
%colormap('summer'); 
axis equal;