Hi,you can use uniroot function. It's included in stats packages for R.> f<-function(x){ 4.512338-log(x/(x+1310.6))/ log(x/(x+185.6))}> uniroot(f,c(0.0001,100000)) > uniroot(f,c(0.001,10000))$root[1] 715.032$f.root[1] 1.365308e-11$iter[1] 9$estim.prec[1] 6.103516e-05