install.packages("nleqslv")
library("nleqslv")

soustava<-function(x) {
  y<-x[1]^2+x[1]-5
}
pocatek<-c(2)
nleqslv(pocatek,soustava)[[1]]