02:二分法求函数的零点

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 }

 

posted on 2017-05-20 18:30  华山青竹  阅读(1703)  评论(0编辑  收藏  举报

导航