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

/*************** 表4.1の数値解とその誤差を求めるプログラム **************

	・二分法を用いて3次代数方程式の実根を求めます.

*************************************************************************/

#define delta pow(10,-10)	//δ(textに対応)
#define epsilon 0.5*pow(10,-5)	//ε(textに対応)
#define a0 -3.			//aの初期値
#define b0  0			//bの初期値
#define kmax 1000		//収束判定条件との比較回数の上限

double f(double x){
	double y;
	y=pow(x,3)+6*pow(x,2)+21*x+32;
	return y;
}

/***** "区間"を引数にとり，根(必ず見つかる)を返す関数 *****/
double NIBUN(double a, double b){
	int k;
	double c;
	for(k=1; k<kmax+1; k++){
		c = (a+b)/2;
		if(fabs( f(c) )<delta || fabs(a-b)<epsilon)
			return c;	//根を返す
		if(f(c)*f(a)<0){
//					a=a;
					b=c;
		}else{
					a=c;
//					b=b;
		}
		printf("| %3d回 | %13.5e | %9.7lf | %13.5e |\n"						, k, fabs(a-b), c, f(c));
	}
	printf("error!!\n");
	exit(1);			//エラーによる強制終了
}

int main(void){
	int i;		
	double a, b;			//区間[a,b]
	double root;			//根を入れる変数
	a = a0;
	b = b0;
	if(f(a)*f(b)>0)
		printf("区間が不適当です\n");
	else{
		root = NIBUN(a,b);
		printf("\n root=%8.7lf\n",root);
	}
	return 0;
}