Rc:=sqrt((Ro**2+Ro*r+r**2)/3);