
% plot circular inclusions (assumes setprob.m has already been called)
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
yrbcolormap
colormap([1 1 1])  % for white background
daspect([1 1 1])
caxis auto
%showgridlines

return

% for mesh plot of solution:
[xc,yc] = meshgrid(xcenter,ycenter);
[xp,yp] = mapc2p(xc,yc);
figure(2)
mesh(xp,yp,q')
colormap([0 0 1])
figure(1)

