
% plot the grid

setprob  % reads in data defining the circlular inclusions

n = 60;
[xc,yc] = meshgrid(linspace(0,1,n+1), linspace(0,1,n+1));
[xp,yp]=mapc2p(xc,yc);
pcolor(xp,yp,0*yp);   colormap([1 1 1]);   daspect([1 1 1])

axis([0 1 0 1])

% plot circles:
hold on
theta = linspace(0,2*pi,200);
for i=1:ncirc
   xcirc = x0circ(i) + r1circ(i)*cos(theta);
   ycirc = y0circ(i) + r1circ(i)*sin(theta);
   plot(xcirc,ycirc,'r','linewidth',2.0)
   end

hold off

