%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Generating problem data 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
randn('state',1)
N = 4;                  % number of balls 
n = 4;                  % variable size 
% centers of the N balls 
c = 10*randn(n,N/2); 
c = [c -c];
% radii of the N balls 
R = norms(c) + .1; 
% translation
offset = 20*randn(n,1);
c = c - repmat(offset, 1,N);

niter = 500; 
tol =1e-6;
