http://noi.openjudge.cn/ch0111/02/
02:二分法求函数的零点
- 总时间限制:
- 1000ms
- 内存限制:
- 65536kB
- 描述
-
有函数：
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位。
- 样例输入
-
无
- 样例输出
-
不提供
1 #include<stdio.h> 2 #include<math.h> 3 double f(double x) 4 { 5 double ans=0; 6 ans=pow(x,5)-15*pow(x,4)+85*pow(x,3)-225*pow(x,2)+274*x-121; 7 return ans; 8 } 9 int main(int argc, char *argv[]) 10 { 11 double x,y,mid,fx,fy,fm; 12 x=1.5; 13 y=2.4; 14 mid=(x+y)/2; 15 fx=f(x); 16 fy=f(y); 17 fm=f(mid); 18 while( fabs(fm-0) > 1e-6) 19 { 20 if(fx*fm<0) { y=mid; fy=f(y); } 21 else { x=mid; fx=f(x); } 22 mid=(x+y)/2; 23 fm=f(mid); 24 } 25 printf("%.6lf\n",mid); 26 return 0; 27 }
