如何在e函数中避免NaN？

Question

我正在尝试为q1和q2求解这些方程。然而，输出给了我＆＃34; NaN＆＃34;。有没有人知道如何避免这些？

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

提前致谢!!!