r - 如何在 e 函数中避免 NaN？

Question

我正在尝试为 q1 和 q2 求解这些方程。但是，输出给了我“NaN”。有谁知道如何避免这些？

k2<-2
k1<-2
a<- 20
b<- -0.03
w <- 1.1

model <- function(q) {
pi1 <- -exp(-a*(w*(100*q[1]-0.5*q[1]^2-0.5*q[1]*q[2])-k1)*(-a*(w*(100-q[1]-0.5*q[2])-k1)))
pi2 <- -exp(-a*(w*(100*q[2]-0.5*q[2]^2-0.5*q[1]*q[2])-k2)*(-a*(w*(100-q[2]-0.5*q[1])-k2)))
c(pi1 = pi1, pi2 = pi2)
}

ss1 <- rootSolve::multiroot(f = model, start = c(1, 1))
ss1

提前致谢！！！