私は、パウエルの方法を使用してf(x)=0; の解を見つけるEigenhybrd1のminpackのポートを使用しています。この場合、ヤコビアンは数値的に計算されます。
ソルバーはxtol、パラメーター、相対許容誤差を受け入れます。これは、デフォルトsqrt(eps)で、特定の数値型 (double私の場合) に対応しています。
ほとんどの場合、 のような残差ベクトル ノルムで満足のいく解が得られ1.446575e-9ます。ただし、523.543 などの規範を与える解が受け入れられることもありますが、これはxtol実際には何を意味するのか疑問に思います。現在の残差ノルムが比較される基準値は何ですか?