二分法求解方程的解

运用二分法求解方程，方程必须是单调的，这样才能保证在寻找的区间内找到解

f(x)=x*x*x-5*x*x+10*x-80;在0到100之间的根，|f(x)|<1e-6;

#include<iostream>

#include<algorithm>
#include<cstdio>
#include<cmath>
using namespace std;
#define EPS 0.000001 //认为方程的值小于这个精确度即为其解
double y(double x)
{
    return x*x*x-5*x*x+10*x-80;
}
int main()
{
    double x1=0,x2=100; //方程的解所在的区间范围
    double root=x1+(x2-x1)/2; //防止想x1+x2超出范围
    while(fabs(y(root))>EPS) //只要不小于精确度，便继续循环
    {
        if(y(root)>0) //根据值的大小，不断调整查找区间
        {
            x2=root;
        }
        else
        {
            x1=root;
        }
        root=x1+(x2-x1)/2;
    }
    printf("%.8lf\n",root);    //保留八位小数，输出结果
    return 0;
}