二分法

二分法

任取两点x1和x2,判断(x1,x2)区间内有无一个实根。如果f(x1)和f(x2)符号相反，说明（x1，x2）之间有一个实根。取（x1，x2）的中点x，检查f(x)与f(x1)是否同符号，如果不同号，说明实根在（x,x1）区间，这样就已经将寻找根的范围减少了一半了。然后用同样的办法再进一步缩小范围。再找x1与x2（x2=x）的中点“x”，并且再舍弃其一半区间。如果f（x）与f(x1)同号，则说明根在（x,x2）区间，再取x与x2的中点，并舍弃其一半区间。用这个办法不断缩小范围，直到区间相当小为止。

float fun(float x) // 计算f(x)
{...}

int main()
{
float x1, x2， x; //需要给定初始值，即一个大概的实根区间，要求异号
...
while(fabs(x2-x1)>1e-6) //假设 fun(x1)>0, fun(x2)<0
{
x=(x1+x2)/2;
if(fun(x)>1e6)x1=x;
else x2=x;
}
...
return 0;
}