質問: 1; j < N; j++){ printf("%d \t %.12e \t %.2e\n" , j, X_n[j], fabs(X_n[j]-X_n[k])); } return 0; }
回答: 最後のprintf文は、最終的な近似解とそれまでの近似解との誤差を出力してるかと思われます。 しかし、kは1つ目のfor文内でしか値が適用されない変数なので、fabs(X_n[j]-X_n[k])では適切な値になりません(環境によってはコンパイルエラーになります、自分の環境ではそうなりました) 最後のprintf文について、X_nkをkを使わず以下のように書き換えるとうまくいくんじゃないかと思います。 ``` #修正前 printf("%d \t %.12e \t %.2e\n", j, X_n[j], fabs(X_n[j] - X_n[k])); #修正後 printf("%d \t %.12e \t %.2e\n", j, X_n[j], fabs(X_n[j] - X_n[N-1])); ``` ご参考になれば幸いです。