noi.openjudge 02:二分法求函数的零点

http://noi.openjudge.cn/ch0111/02/

描述

有函数：

$f(x) = x^5 - 15 * x^4+ 85 * x^3- 225 * x^2+ 274 * x - 121$f(x)=x5−15∗x4+85∗x3−225∗x2+274∗x−121

已知 f(1.5) > 0 , f(2.4) < 0 且方程 f(x) = 0 在区间 [1.5,2.4] 有且只有一个根，请用二分法求出该根。

输入

无。

输出

该方程在区间[1.5,2.4]中的根。要求四舍五入到小数点后6位。

代码

#include<iostream>
#include<math.h>
#include<stdio.h>
using namespace std;

double f(double x)
{
	return x*(x*(x*(x*(x-15)+85)-225)+274)-121;
}
int main()
{
//	cout<<f(1.84901590)<<" "<<f(1.84901591)<<endl;
	double mid,L,R;
	double y;
	double eps=1e-6;
	L=1.5;
	R=2.4;
	mid=L+(R-L)/2;
		
	do{
		if(f(mid)>0)
			L=mid;
		else 
			R=mid;
			
		y=mid;
		mid=L+(R-L)/2;
		y=mid-y;
	}while(fabs(y)>eps);
	
	printf("%.6lf",mid);
	return 0;
}