这段二分法的代码只要迭代9步，但是运行之后一直busy，不知道哪里出了问题

blci<-function(r,n,alpha=.05,eps=1e-3){
  if(r<=0) pl<-0 else{
    plu<-r/n
    pll<-0
    pl<-(pll+plu)/2
    u<-1-pbinom(r-1,n,pl)
    i<-1
    while((plu-pll)/plu>eps){
      i<-i+1
      if(u>alpha){
        plu<-pl
        pl<-(pll+plu)/2
        u<-1-pbinom(r-1,n,pl)
      }
        else{
        pll<pl
        pl<-(pll+plu)/2
        u<-1-pbinom(r-1,n,pl)
      }
    }
  }
    c(lcl=pl,neval=i)
}

blci(5,15,alpha=.05,eps=1e-2)