我有一个方程组来求解x1和x2：
2*erf(A*x2)*exp(B*x1^2+C*x2^2)-D*erf(x1)=0
2*erf(D*x2)*exp(E*x1^2+F*x2^2)-G*erf(x1)-erf(H*x2)*exp(I*x1^2)=0
其中B和C是负数，而其余常数是正数。
我尝试使用 GSL 库中的多根函数在 C 中解决它。似乎多根无法解决系统问题。我为非线性同时多项式函数（其中没有指数部分）尝试了相同的多根实现，它运行良好并给出了正确的结果。关于如何解决这样一个系统的任何建议？或者任何其他可用于此类方程的强大库？