n = 14;
k = 4;
lim = 10;
Cgiven = [-lim, 0; 0, -lim; 0, lim; lim, 0];
Rgiven = ones(k, 1) * 2;

Gindexes = [0.000000000000000000e+00,1.300000000000000000e+01; ...
1.000000000000000000e+00,7.000000000000000000e+00; ...
1.000000000000000000e+00,1.200000000000000000e+01; ...
2.000000000000000000e+00,1.100000000000000000e+01; ...
3.000000000000000000e+00,1.200000000000000000e+01; ...
4.000000000000000000e+00,5.000000000000000000e+00; ...
4.000000000000000000e+00,9.000000000000000000e+00; ... 
5.000000000000000000e+00,8.000000000000000000e+00; ...
5.000000000000000000e+00,9.000000000000000000e+00; ... 
5.000000000000000000e+00,1.100000000000000000e+01; ...
6.000000000000000000e+00,1.000000000000000000e+01; ... 
6.000000000000000000e+00,1.200000000000000000e+01; ...
7.000000000000000000e+00,1.300000000000000000e+01; ...
8.000000000000000000e+00,1.300000000000000000e+01; ...
1.000000000000000000e+01,1.100000000000000000e+01; ... 
1.000000000000000000e+01,1.200000000000000000e+01] + 1;

%% Feel free to use this code for plotting
%%% centers: n-by-2 matrix of circle centers
%%% radii: length-n vector of radii
% figure;
% viscircles(centers(1:k,:), radii(1:k), 'EdgeColor', 'r'); hold on
% viscircles(centers(k+1:n,:), radii(k+1:n), 'EdgeColor', 'b');
% scatter(centers(1:k, 1), centers(1:k, 2), 100, 'ro', 'fill');
% scatter(centers(k+1:n, 1), centers(k+1:n, 2), 100, 'bo', 'fill');
% for i = 1 : size(Gindexes, 1)
%     x1 = centers(Gindexes(i, 1),:); a = x1(1); b = x1(2);
%     x2 = centers(Gindexes(i, 2),:); c = x2(1); d = x2(2);
%     plot([a, c], [b, d], 'k-', 'LineWidth', 2);
% end
% xlim([-12.5, 12.5]); ylim([-12.5, 12.5]);
% pbaspect([1, 1, 1]); hold off