[R-br] Ponto de intercepção entre duas funções no R

walmes . walmeszeviani em gmail.com
Quinta Outubro 16 18:12:26 BRT 2014


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.
​
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141016/38ce41dd/attachment.html>


Mais detalhes sobre a lista de discussão R-br