
global zout cout ncirc x0circ y0circ r1circ r2circ zin cin;

fid = fopen('setprob.data');
data = fscanf(fid,'%g',inf);
zout = data(1);
cout = data(2);
ncirc = data(3);
for i=1:ncirc
   x0circ(i) = data(4+6*(i-1));
   y0circ(i) = data(5+6*(i-1));
   r1circ(i) = data(6+6*(i-1));
   r2circ(i) = data(7+6*(i-1));
   zin(i) = data(8+6*(i-1));
   cin(i) = data(9+6*(i-1));
   end
  
fclose(fid);

