同じ出力の値を二度表示しないようにするプログラミング

Question

以下の近似値を求めるプログラムを実行し
x0=1,xn=2を入力すると

範囲[x0,xn]
x0=1
xn=2
x0 xn
1.000000 2.000000
1.500000 2.000000
1.500000 1.750000
1.500000 1.750000
1.625000 1.750000
1.625000 1.687500
1.625000 1.687500
1.656250 1.687500
1.671875 1.687500
1.671875 1.679688
1.671875 1.679688
1.671875 1.675781
1.671875 1.675781
1.671875 1.673828
1.671875 1.673828
1.672852 1.673828

と出力されますが下記太字のように同じ値が出力されています。

x0 xn
1.000000 2.000000
1.500000 2.000000
1.500000 1.750000
**1.500000 1.750000 **
1.625000 1.750000
1.625000 1.687500
1.625000 1.687500
1.656250 1.687500
1.671875 1.687500
1.671875 1.679688
1.671875 1.679688
1.671875 1.675781
**1.671875 1.675781 **
1.671875 1.673828
**1.671875 1.673828 **
1.672852 1.673828

同じ値が２度以上表示されないようにするにはどう変更すればよいか
ご存じの方ご教示お願い致します。

#include<stdio.h>
#include<math.h>
#include<stdlib.h> //exit関数のために必要

#define IMAX 1000 //計算回数
#define EPS 1.0e-5 //収束条件

double f(double x);

int main(void){

    double x0,xn,xc;
    int i;

     printf("範囲[x0,xn]\nx0=");
    scanf("%lf",&x0);

    printf("xn=");
    scanf("%lf",&xn);

    printf("   x0      xn\n");
    printf("%.6f %.6f\n",x0,xn);


    xc=(x0+xn)/2;

    if(f(x0)*f(xn)<0.0){    //範囲の中に解があれば計算開始
        for(i=0;i<IMAX;i++){
            if( f(x0)*f(xc)<0.0 ){ //範囲を縮める方向の判定
                xn=xc;
                printf("%.6f %.6f\n",x0,xn);

            }else x0=xc;
               printf("%.6f %.6f\n",x0,xn);
               
            xc=(x0+xn)/2;

            if( fabs(x0-xn) < 0.001) {
                
                exit(1);                           
            }
        }printf("解が収束しませんでした\n");
    }printf("[x0,xn]には解が存在しません\n");

    return 0;
}
double f(double x){
    return x*x*x-3*x*x+7*x-8;
}    

Answer 1

中括弧の付け忘れのようですね。
（x0=xc; の後のprintfが毎回実行されていますよ。）

            if( f(x0)*f(xc)<0.0 ){ //範囲を縮める方向の判定
                xn=xc;
                printf("N:%.6f <%.6f\n",x0,xn);

            }else{ 
                x0=xc;
                printf("0:%.6f> %.6f\n",x0,xn);
            }

ついでながら、同じフォーマットの表示には「見分けるための目印」を付けておくと
デバッグがはかどりますよ。

   x0      xn
s:1.000000  2.000000
0:1.500000> 2.000000
N:1.500000 <1.750000
0:1.625000> 1.750000
N:1.625000 <1.687500
0:1.656250> 1.687500
0:1.671875> 1.687500
N:1.671875 <1.679688
N:1.671875 <1.675781
N:1.671875 <1.673828
0:1.672852> 1.673828

どちらから寄せたか判りやすい…んじゃないかと思います。
元のロジック（括弧なし）で実行すると、
「（小さい方から）寄せてないのに表示が出た」と気付けます（たぶん）

Comments

[KeiD]

ありがとうございます。