
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',1.5)
   end
   
axis square
colormap([1 1 1])
%showgridlines
clear mapc2p

hold off
