global y;
global n;
global d;
global Nodes1;
global K;

n = 400;    % sample size 
H = eye(n);
K = eye(n);

d = 0.12;   % d2 - d1, radius difference
%d1 is the inside circel radius, and d2 is the outside radius, d2>d1 
%satisfy:  delta_2 - delta_1 = 0.04; pi * (delta_2 ^2 + delta_1 ^2) = 4
%here d1 = 0.7356, d2 = 0.8556
d1 = sqrt(8/pi - d^2)/2 - d/2;  
d2 = d1 + d;


Nodes1 = zeros(n,2);
y = zeros(n,1);

rand('state', 100);
for i=1:n
  Nodes1(i,:) = 2*rand(1,2)-1;
end

y = gen_y(d1,d2);

