求解函数的零点
- 问题描述:数值分析中一个重要的问题是对某个任意的函数f找出方程f(x)=0的一个解。如果该函数是连续的并有两个点low和high使得f(low)和f(high)符号相反,那么在low和high之间必然存在一个根,并且这个根可以通过折半查找求得。写出一个函数,以f、low、high为参数,并且解出一个零点。为保证能够正常终止,你必须要做什么?
- 代码实现
#include<iostream>
using namespace std;
const double tol = 0.00001;
double f(double x)
{
return x*x - x;
}
double bisection( double (*f)(double), double low, double high)
{
double mid = (high + low)/2.;
double value;
if ((high - low) < tol)
return (mid);
value = f(mid);
if (value*f(low) < 0)
return bisection(f, low, mid);
else
return bisection(f, mid, high);
}
int main()
{
cout<<bisection(f, .1 , 10)<<endl;
return 0;
}