function [flag,xi,yi] = f_interbis(a,b,xa,ya,rg)


delta=4*( ((rg^2)*(1+(a^2))) - (b^2) )

if (delta > 0)
   
	flag=1;
	solx_1=( -2*a*b-sqrt(delta) ) / ( 2*((a^2)+1) );
	solx_2=( -2*a*b+sqrt(delta) ) / ( 2*((a^2)+1) );


	soly_1=a*solx_1+b;
	soly_2=a*solx_2+b;

	if ( abs(xa-solx_1) < abs(xa-solx_2) )
   xi=solx_2;
   yi=soly_2;
	else
   xi=solx_1;
   yi=soly_1;
   end

else
   xi=xa;
   yi=ya;
   flag=0;
   
end