16 Out
2014
16 Out
'14
21:12
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.