割線法（セカント法）

前提・実現したいこと

C++で割線法のプログラムを作りたいのですが、何も出力されません。おそらくwhile文の中のif条件のところに問題があると思っていますが、どのように直すべきでしょうか。f(x) = x^5 + ax^4 + bx^3 + cx^2 + dx + e = 0で誤差0.000001以下の時にwhile文を抜けて答えを出力したいです。確認宜しくお願いします。

該当のソースコード

C++
1#include<iostream>
2#include<cmath>
3using namespace std;
4int a, b, c, d, e;
5float f(float x){ //f(x) = x^5 + ax^4 + bx^3 + cx^2 + dx + e
6    return pow(x, 5) + a * pow(x, 4) + b * pow(x, 3) + c * pow(x, 2) + d * x + e;
7}
8float xpoint(float x1, float x2){
9    return f(x2) * ((x2 - x1) / (f(x2) - f(x1)));
10}
11float root(float x1, float x2){
12    return x2 - xpoint(x1, x2);
13}
14int main()
15{
16    cin >> a; getchar(); cin >> b; getchar(); cin >> c; getchar(); cin >> d; getchar(); cin >> e;
17    float x1, x2, x3;
18    cin >> x1;
19    getchar();
20    cin >> x2;
21    while(1){
22        x3 = root(x1, x2);
23        if(fabs(xpoint(x1, x2)) < 1e-6) break;
24        x2 = x3; x1 = x2;
25    };
26    printf("%.3f", x3);
27    return 0;
28}

入力

1,1,1,1,1
-2,0

出力の目標

-1.000

行動規範の内容に同意します

回答2件

ベストアンサー

x2 = x3; x1 = x2; を
x1 = x2; x2 = x3; にしないと、
x1 と x2 が同じ値になって、ゼロ除算を実行してしまいますよ。

追記
速くする方法ですか。
・pow を使わない。
・同じ値 x で何度も f(x) を計算するのをやめる。
・内部で float から double への変換があって、double のほうが速いかも。
・int と double を混ぜると、int から double への変換が余計に入る。

C++
1#include <iostream>
2#include <cmath>  // fabs
3using namespace std;
4
5double a, b, c, d, e;
6
7double f(double x) { //f(x) = x^5 + ax^4 + bx^3 + cx^2 + dx + e
8    return ((((x + a)*x + b)*x + c)*x + d)*x + e;
9}
10
11int main()
12{
13    double x1, x2, x3;
14	char p;
15    cin >> a >> p >> b >> p >> c >> p >> d >> p >> e >> x1 >> p >> x2;
16	double f1 = f(x1), f2 = f(x2);
17    while (1) {
18		double xp = f2 * (x2 - x1) / (f2 - f1);
19        x3 = x2 - xp;
20        if (fabs(xp) < 1e-6) break;
21		x1 = x2, x2 = x3, f1 = f2, f2 = f(x3);
22    }
23    printf("%.3f", x3);
24    return 0;
25}