% given N points and a triangulation of these points
% find a radius of the sector associated to each vertex
%  so that the omitted central region has intervals on 
% one each edge of the triangle, and so the the disks
% with these intervals as diameters do not to touch 
% the other two side of the triangle.


function radii=find_radii2(x,y,tri)
   N=size(tri,1);
   radii=zeros(N,3);
   for n=1:N
       a=x(tri(n,1))+i*y(tri(n,1));
       b=x(tri(n,2))+i*y(tri(n,2));
       c=x(tri(n,3))+i*y(tri(n,3));
       theta_a=abs(imag(log( (b-a)/(c-a) ) ));
       theta_b=abs(imag(log( (a-b)/(c-b) ) ));
       theta_c=abs(imag(log( (a-c)/(b-c) ) ));
       radii(n,1)=find_sec_radius(a,b,c,theta_a,theta_b,theta_c);
       radii(n,2)=find_sec_radius(b,c,a,theta_b,theta_c,theta_a);
       radii(n,3)=find_sec_radius(c,a,b,theta_c,theta_a,theta_b);
   end
return

