Regula falsi

f <- function(x) (x^x)-100
F1x <- function(x) 1-exp(-x)
regfals <- function(a,b) {
error<-1
while (error > .1) {
x<-(f(b)*a-f(a)*b)/(f(b)-f(a))
if (f(x)*f(a) < 0) {b <- x}
else {a <- x}
error<-abs(f(x)/F1x(x))
print(x)
}
}
regfals(3,4)

## [1] 3.318777
## [1] 3.474996
## [1] 3.545438
## [1] 3.575675
## [1] 3.588346
## [1] 3.593599
## [1] 3.595767
## [1] 3.596661
## [1] 3.597028