c言語 ニュートン法

Question

#include <stdio.h>
#include <math.h>

#include <stdio.h>
#include <math.h>

double g(double x)
{
    return exp(-x)-sin(x);
}
double dg(double x) {
    return -exp(-x)-cos(x);
}

double newton2(double xn)
{
    int i = 0;
    double xo;

    while (fabs(xn-xo)>1e-10){
        xo = xn;
        xn = -1 * g(xo) / dg(xo) +xo;
        i++;
    };
    printf("result:%lf,試行回数:%d\n",xn,i);
}
int main()
{
   
    printf("ニュートン法\n");
    newton2(2.0);
   
    return 0;
}

結果
ニュートン法
result:25.132741,試行回数:8
が出てきます。
実際、出力したいのは、result:0.588533です。
どこが間違っているのか教えていただければ幸いです。

Answer 1

蛇足ではありますが、理由としては以下になります。

Comments

[fana]

（縦軸を"y", 横軸を"x"と呼ぶとして…）
x=2 における接線もプロットする，且つ，グラフの横軸範囲をその接線とy=0の交点が見れるくらいの範囲にしてやる，
…みたいな感じになっていればより親切(?)だった気がします．

[vel3633]

詳しく解説助かりました。納得できました！

Answer 2

初期値を 2.0 ではなく、1.0 にしてみてください。

Comments

[vel3633]

コメントありがとうございます。解決しました。