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

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

描述

有函数：

$f(x)=x^5-15\ast x^4+85\ast x^3-225\ast x^2+274\ast 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;
}