function v = BR_2circles(x,c,r1,r2)
%

v = - norm(x)*norm(c-x)/(r1*r2);
v = 1*(norm(x)/r1<=1)*1*(norm(c-x)/r2<=1);
