R/falsaposi.R

In luisxsuper/metodNum: Metodos Numericos

#' @export
falsa.pos <- function(fx, a, b, tol) {
   #' @author Luis Fernando Maron I.
   curve(fx,a,b,type="l",col="cyan4", ylab = "y")
   abline(h=0,col= c("#333333C7"), lty = 2)
   if (fx(a) * fx(b) > 0)
      cat("Por favor intente con otro intervalo [a,b]\n")
   else{
      x <- 0
      i <- 0
      while (abs(fx(x)) > tol) {
         i <- i + 1
         x <- (a * fx(b) - b * fx(a)) / (fx(b) - fx(a))
         if (fx(a) * fx(x) < 0)
            b <- x
         else
            a <- x
         cat(sprintf("\n\tIteracion: %2d", i))
         cat(sprintf("\tRaiz: %f \n", x))
      }
   }
   cat("\n\t\t\t\t\t-->By Luis Maron<--\n")
}