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.