完善后的二分法(C++版)

本程序为本人参考一些资料后编写，在visual c++6.0下调试通过。

[@more@]
/二分法求一个未知数方程的根f(x)=0,x属于[a,b]
/改进：除了显示每次计算的小区间外，还根据给定的精度计算了所需的次数k
/其中k>[log(b-a)-log(2*err)]/log(2);log=ln;
/完善二分法的模块

#include
#include
#include
using namespace std;
//
double fun(double x);
void RootByBisection(double(*pf)(double x));
//
//
int main()
{
RootByBisection(fun);
return 0;
}
//
///可以在下面修改函数
double fun(double x)
{
return (x*x*x-x-1);
}
//
void RootByBisection(double(*pf)(double x))
{
double x1,x2,err;
int flag;///记录迭代次数
cout<>x1;
cout<>x2;
cout<>err;
flag=static_cast((log(x2-x1)-log(2*err))/log(2))+2;
cout< cout.precision(10);
cout< < flag=0;
while(fabs(x2-x1)>err)
{
if(pf(x1)*pf((x1+x2)/2)<0)
x2=(x1+x2)/2;
else
x1=(x1+x2)/2;
flag=flag+1;
cout< < }
cout< cout<}
//
//

来自 “ ITPUB博客 ” ，链接：http://blog.itpub.net/11355887/viewspace-974943/，如需转载，请注明出处，否则将追究法律责任。

转载于:http://blog.itpub.net/11355887/viewspace-974943/