float rma(float a, float b,float r){
float m;
m=(a+b)/2;
if ((b-a)<0.001){return m;}
else
if(m*m>r){return rma(a,m,r);}
else
{return rma(m,b,r);}



}