回答編集履歴
1
コードがC++だった→Cなコードに修正
answer
CHANGED
|
@@ -19,10 +19,11 @@
|
|
|
19
19
|
//デバッグ用に経過を表示しよう
|
|
20
20
|
printf( "[i=%d] x=%.12f, f(x)=%.12f, df(x)=%.12f\n", i, X_n, F, dF );
|
|
21
21
|
|
|
22
|
-
//収束判定．
|
|
22
|
+
//収束判定．
|
|
23
|
+
//dFについても別途必要かも（極値付近でまともに解を更新計算できないと思うので）．
|
|
23
|
-
if(
|
|
24
|
+
if( fabs(F)<epilon )
|
|
24
25
|
{
|
|
25
|
-
|
|
26
|
+
printf"(Finish by Thresh)\n";
|
|
26
27
|
break;
|
|
27
28
|
}
|
|
28
29
|
//解の更新
|