Veja a documentação da função rootSolve::uniroot(). CMR

f1 <- function(x, y){ x^2-3*y }
f2 <- function(x, z){ 4*z-x^2 }
x <- seq(1,10,0.5); y <- 3; z <- 15

curve(f1(x, y=y), 1, 10)
curve(f2(x, z=z), add=TRUE, col=2)
abline(v=6)

f <- function(x, y, z){
  x^2-3*y-(4*z-x^2)
}

curve(f(x, y=y, z=z), 0, 10); abline(h=0, lty=2)
abline(v=6)

require(rootSolve)

ur <- uniroot(f=f, y=y, z=z, interval=c(0,10))
str(ur)

curve(f(x, y=y, z=z), 0, 10); abline(h=0, lty=2)
abline(v=ur$root)

curve(f1(x, y=y), 1, 10)
curve(f2(x, z=z), add=TRUE, col=2)
abline(v=ur$root)

À disposição.
Walmes.
​
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código 
mínimo reproduzível.

Responder a