次のような L+2 変数を持つ連立方程式があります。
zL = f(z0, ..., z(L-1)), z0 = g(zL), zl = h_l(z(l-1)) for 0<l<L, and 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) のみに制限する方法はありますか? または、正の制約を持つ非線形連立方程式を解く他の方法はありますか?