我有一个方程组来解决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中解决它。似乎multiroot无法解决系统问题。我为非线性同时多项式函数（没有指数部分）尝试了相同的多根实现，并且它运行良好并给出了正确的结果。 关于如何解决这样一个系统的任何建议？或者可以用于此类方程的任何其他强大的库？