#include<stdio.h>
#include<math.h>
/*
x^3-4x^2+3x+2=0 の 2.1 と 3.0 の間ある解の数値計算
期待される答: 1+2^0.5〜2.414213562373095 
gcc -o root.exe root.c でコンパイル
 */
#define LEFT 2.1 
#define RIGHT 3.0
#define N 20  /* 精度 */
extern double f(double x);
main()
{
	double a,b,c, fa, fb, fc;
	int i;
	a=LEFT; b=RIGHT;
	for(i=0;i<N;i++)
		{
			c=(a+b)/2.0;
			if (f(a)*f(c)<0.0)
				b=c;
			else 
				a=c;
		}
	printf("%f と %f の間にある方程式の数値解（の1つ）は x = %lf です\n", LEFT, RIGHT, c); 
	printf("計算の誤差は ± %lf 未満です\n", pow(2,-N)*(RIGHT-LEFT));
}

double f (x)
		 double x;
{ return (x*x*x-4*x*x+3*x+2);
}