02:二分法求函数的零点

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

02:二分法求函数的零点

总时间限制: 

1000ms

 

内存限制: 

65536kB

描述

有函数：

f(x) = x⁵ - 15 * x⁴+ 85 * x³- 225 * x²+ 274 * x - 121

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

输入

无。

输出

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

样例输入

无

样例输出

不提供

#include<stdio.h>
#include<math.h>
double f(double x)
{
    double ans=0;
    ans=pow(x,5)-15*pow(x,4)+85*pow(x,3)-225*pow(x,2)+274*x-121;
    return ans;
}
int main(int argc, char *argv[])
{
    double x,y,mid,fx,fy,fm;
    x=1.5;
    y=2.4;
    mid=(x+y)/2;
    fx=f(x);
    fy=f(y);
    fm=f(mid);
    while( fabs(fm-0) > 1e-6)
    {
        if(fx*fm<0) { y=mid; fy=f(y); }
        else { x=mid; fx=f(x); }
        mid=(x+y)/2;
        fm=f(mid);
    }
    printf("%.6lf\n",mid);
    return 0;
}