本程序为本人参考一些资料后编写,在visual c++6.0下调试通过。
[@more@]//////////////
////牛顿迭代法,又称作切线法,其具体原理为
////求方程f(x)=0的近似根x_k,f(x)=f(x_k)+f'(x_k)(x-x_k)+o(x-x_k)
////从而转化为求f(x_k)+f'(x_k)(x-x_k)=0的根
////牛顿迭代公式为:x_(k+1)=x_k-f(x_k)/f'(x_k)
////相应的迭代函数为:g(x)=x-f(x)/f'(x),x=g(x)
////改进:输出每一次迭代的结果,将迭代模块做全面
//////////////
#include
#include
using namespace std;
///////////
double fun(double x);
double dfun(double x);
void NewtonIterative(double(*pfun)(double x),double (*pdfun)(double x));
//////////
//////////
int main()
{
NewtonIterative(fun,dfun);
return 0;
}
//////////
//////////
double fun(double x)
{
return (x*x*x-x-1);
}
//////////
double dfun(double x)
{
return (3*x*x-1);
}
//////////
void NewtonIterative(double(*pfun)(double x),double (*pdfun)(double x))
{
int maxflag,k=1;
double x0,x;
double err;
cout<>maxflag;
cout<<" 请输入精度要求:err=";cin>>err;
cout<<" 请输入初始值x0=";cin>>x;
cout< do
{
x0=x;
x=x0-pfun(x0)/pdfun(x0);
cout<<" 第"< ++k;
}while(fabs(x-x0)>err&&k if(k>=maxflag)
cout< else
{
cout< cout<<" 牛顿迭代结果:近似根为root="< cout<<" 函数值近似为:f(root)="< }
}
//////////
//////////
来自 “ ITPUB博客 ” ,链接:https://blog.itpub.net/11355887/viewspace-974946/
,如需转载,请注明出处,否则将追究法律责任。