# una delle soluzioni (a 9 cifre) di x^2-x=3.2 usando la bisezione
g <- function (x) x^2-x-3.2
plot(g,-5,5); abline(h=0,col="blue")
a <- 1; b <- 3
for(i in 1:30){h <- (b-a)/2; ifelse(g(a)*g(a+h)>0,a <- a+h,b <- b-h); print(c(a,b),9)}
points(a,g(a))