我有以下方程组（简化版） .
y = x*scipy.special.gdtr(a, b, f(x,y))
x = y*scipy.speical.gdtr(c, d, f(x,y))
这里 scipy.special.gdtr 指的是伽玛分布的CDF， a ， b 是伽马CDF所采用的相应的两个参数 .
我的问题如下：
-
在这种情况下，是否有首选求解器？
fsolve或brentq或其他？换句话说，选择求解器的标准是什么？ -
到目前为止我使用了
fsolve，我认为它有两个很大的缺点，
i），每当我输入不同的初始猜测时，它总是给我不同的解决方案 .
ii），它经常告诉我迭代不会收敛到单个值（有几个案例有效，它似乎取决于a和b的具体值） .
_999_当非线性方程与CDF结合时，似乎很难获得一个明确的解决方案 .
我看过这里的各种帖子，包括：
和
How to solve a pair of nonlinear equations using Python?
这些帖子没有具体解决我的问题 . 如果有任何在python中解决非线性方程和统计系统经验的人能够回答我的问题，将不胜感激 .