二分法求方程的根

【题目】来源：xmuoj

【二分法关键原理】

对于在闭区间[a, b]上连续不断且满足f(a) * f(b) < 0的函数y=f(x)，通过不断地把函数f(x)的零点所在的区间一分为二，使区间的两个端点逐步逼近零点，进而得到零点近似值。由于每次迭代都将搜索区间减半，因此二分法的收敛速度相对较快。

【二分法实现步骤】

1.确定初始区间[a,b]，满足f(a)f(b)<0;

2.区间中点c=(a+b)/2，计算f(c)

        ·若f(c)=0，则c为所求解；

         ·若f(c)与f(a)同号，则零点在[c,b]之间，更新a=c；

         ·若f(c)与f(b)同号，则零点在[a,c]之间，更新b=c；

重复步骤2，直到满足预设精度要求或达到最大迭代次数；

【核心代码】

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

const double eps=1e-9;//定义我们计算的精度 

double f(double x){
	return x*x*x-5*x*x+10*x-80;//定义函数方程 
}

int main(){
	double a=0,b=10,c;//先粗略确定区间两端a,b的值，c为中点
	while(fabs(a-b)>eps){
		c=(a+b)/2;
		if(fabs(f(c))<eps)break;//c为零点 
		else if(f(c)*f(a)>0)a=c;//更新 
		else if(f(c)*f(b)>0)b=c;//更新 
	}
	c=(a+b)/2;//最终二分法得到的零点
	printf("%.9f",c); 
}