質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

89.52%

Newton法と誤差について

解決済

回答 1

投稿 編集

  • 評価
  • クリップ 1
  • VIEW 296

Co-Min

score 5

Newton法でf(x)=0の解を求めるプログラムをC++で作っています。
f(x)=10sin(x)+e^x
反復回数が6回になると思うのですが、そうすると誤差の計算(最後のfor文)がおかしくなってしまいます。5回だと上手くいく気もするのですが、、
どっちが正しいのでしょうか?間違っている箇所があれば教えてください。
恐らく近似解は-9.139306700409e-02
誤差は最後の反復で1.85e-11になると思います。

include<stdio.h>

include<math.h>

double f(double x){
return 10*sin(x) + exp(x);
}
double df(double x){
return 10*cos(x) + exp(x);
}

int main(){

int k;
const int Max = 50;//最大反復回数
int N = 0;//反復回数
double x0 = 1.0;//初期値
double eps = 1e-12;//許容誤差
double X_n[Max + 1];//近似解を保存する数列
X_n[0] = x0;

for(int k = 0; k < Max; k++){
if(fabs(f(X_n[k])) < eps){
N = k + 1 ;
printf("反復回数は%d回\n", N);
printf("近似解は%.12e\n", X_n[k]);
break;
}
X_n[k + 1] = X_n[k] - (f(X_n[k]) / df(X_n[k]));
}
for(int j = 1; j < N; j++){
printf("%d \t %.12e \t %.2e\n" , j, X_n[j], fabs(X_n[j]-X_n[k]));
}
return 0;
}

  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 過去に投稿した質問と同じ内容の質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 1

checkベストアンサー

0

最後のprintf文は、最終的な近似解とそれまでの近似解との誤差を出力してるかと思われます。

しかし、kは1つ目のfor文内でしか値が適用されない変数なので、fabs(X_n[j]-X_n[k])では適切な値になりません(環境によってはコンパイルエラーになります、自分の環境ではそうなりました)

最後のprintf文について、X_n[k](最終的な近似解)を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]));

ご参考になれば幸いです。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 89.52%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

同じタグがついた質問を見る