r=a*a/(sqrt(a*a-b*b))