程序如下：
fd<-function(x){
1/x/(1 + x) - log(x)/(1 + x)^2}
fzero<-function(f,a,b,eps=1e-5){
if(f(a)*f(b)>0)
   list(fail="failed")
else{
   repeat{
     if(abs(b-a)<eps) break
     x<-(a+b)/2
     if(f(a)*f(b)<0) b<-x else a<-x
   }
list(root=(a+b)/2,fun=f(x))
}
}
fzero(fd,0,10)
得到的结果是2.5与标准的3.59112相去甚远，望解答