% This script plots the temperature and surface
clf;
patch('Vertices',xy','Faces',tri2nod','FaceVertexCData',T, ...
      'FaceColor','interp','EdgeColor','none');
axis('equal');
title('Temperature');
HC = colorbar;

% Loop over boundary edges and plot them

hold on;
for ii = 1:Nbc,
  kn1 = bedge(1,ii);
  kn2 = bedge(2,ii);
  plot([xy(1,kn1), xy(1,kn2)],[xy(2,kn1), xy(2,kn2)]);
end
hold off;