
function gcplot(g, xy, bndry, c, s)
   if (nargin==4)
      s=50;
   end
   fig = figure;
   set(fig, 'Color', [1 1 1]);
   gplot(g, xy);
   axis off square;
   hold on;
   scatter(xy(:,1), xy(:,2), s, c, 'filled');
   caxis([0 1])
   % Now make the labeled nodes larger
   for i = 1:length(bndry)
      scatter(xy(bndry(i),1), xy(bndry(i),2), 2.5*s, c(bndry(i)));
   end

    
