ITPub博客

首页 > 应用开发 > IT综合 > 完善后的二分法(C++版)

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

原创 IT综合 作者:chenjinglzz 时间:2007-10-07 15:18:33 0 删除 编辑

本程序为本人参考一些资料后编写,在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<<"请输入一元函数f(x)的定义域下界:a=";cin>>x1;
cout<<"请输入一元函数f(x)的定义域上界:b=";cin>>x2;
cout<<"请输入精度要求(推荐取err=1E-6):err=";cin>>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<<"函数值f(root)="<}
//////////
//////////

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

上一篇: 二分法(c++版)
请登录后发表评论 登录
全部评论

注册时间:2007-12-10

  • 博文量
    17
  • 访问量
    335010