と を解く連立方程式が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 ライブラリの multiroot 関数を使用して C で解決しようとしました。マルチルートではシステムを解決できなかったようです。非線形同時多項式関数 (指数部なし) に対して multiroot の同じ実装を試みたところ、うまく機能し、正しい結果が得られました。このようなシステムを解決する方法について何か提案はありますか? または、そのような方程式に使用できる他の強力なライブラリはありますか?