我有一个带有 L+2 变量的方程组，看起来像：
zL = f(z0, ..., z(L-1)), z0 = g(zL), zl = h_l(z(l-1)) 对于 0<l<L, 并且 m(c)= 0 .
这里的函数 f,g,h_l (0<l<L) 和 m 是非线性的，我已经有一个方法fun_eq(z0,z1,...,zL,c)返回：
c(z0-g, z1-h_1,z2-h_2 ..., z(L-1)- h_(L-1), zL-f, m(c))
我想确定所有变量都严格为正的c(z0,z1,...,zL,c)的值。我尝试使用：
library(nleqslv)
nleqslv(some_initial_conditions, fun_eq)
但一些解决方案被证明是否定的。有没有办法将我的解决方案空间限制为仅正 c(z0,z1,...,zL,c)？还是有其他方法可以解决具有正约束的非线性方程组？