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.