function  draw_circle(Nodes,y1)


global d;

n1 = length(y1);

d1 = sqrt(8/pi - d^2)/2 - d/2;
d2 = d1 + d;

a = zeros(n1,2);
b = zeros(n1,2);


for i=1:n1
    if y1(i,1)==1,
         a(i,:) = Nodes(i,:);
    else
         b(i,:) = Nodes(i,:);
    end
end



t = 0:.001:2*pi;
u1 = d1*sin(t);
v1 = d1*cos(t);
u2 = d2*sin(t);
v2 = d2*cos(t);


axis image, plot(a(:,1),a(:,2),'k+');
% title('n = 200, radii are {0.7356, 0.8556}');
title('Example 1: n = 400')
hold on;
axis image, plot(b(:,1),b(:,2),'ro');
hold on;
axis image, plot(u1,v1,'--');
hold on;
axis image, plot(u2,v2,'--');
hold on;
%gtext('+ denotes +1, o denotes -1 in data');
%gtext('A');
%gtext('B');
%gtext('C');

%print data.wc400.ps;
