  function [w1,w2] = Close(z1,z2,delta)
%
% z1,z2 are complex scalars that define a 2-root segment
% delta  is positive real
%
% w1,w2 are complex scalars that define a 2-root segment
%            with |w1 - w2| < delta|z1 - z2|.

bound = delta*abs(z1-z2);
w1 = z1;
w2 = z2;
[steps1,j1] = WhichRoot(w1);
[steps1,j2] = WhichRoot(w2);
while abs(w1-w2) > bound
   m = (w1+w2)/2;
   [stepsm,jm] = WhichRoot(m);
   if jm~=j1
      w2 = m;
      j2 = jm;
   else
      w1 = m;
      j1 = jm;
   end	  
end