¿Cómo resolver una ecuación no lineal usando el paquete nleqslv en R?

Question

Jul 31, 2016, 03:08 PM

¿Cómo resolver una ecuación no lineal usando el paquete nleqslv en R?

Tengo que resolver esta ecuación (por ejemplo, f (x, y) = 0):

library(nleqslv)
target <- function(x)
{
  z = x[1]/(x[1]+x[2])
  y = numeric(2)
  y[1] <- z*exp(-x[2]*(x[2]+z*(1-exp(-x[1]/z))))-0.00680
  y[2] <- z/x[2]*(1-exp(-x[2]))-exp(-x[2])*z/x[1]*(1-exp(-x[1]))-3.43164
  y
}

# Usage
xstart <- c(1,1)
target(xstart)
nleqslv(xstart, target, control=list(ftol=.0001, allowSingular=TRUE),jacobian=TRUE,method="Newton")

usando R con nleqslv u otro tienes otros :)

Gracias