我有以下方程组（简化版）。
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的具体值
当非线性方程与CDF结合时，似乎很难获得明确的解决方案。
我在这里查看过各种帖子，包括：
和
How to solve a pair of nonlinear equations using Python?
这些帖子没有具体解决我的问题。如果有任何在python中解决非线性方程组和统计系统经验的人都能回答我的问题，我将不胜感激。