加速迭代.cpp

来自「数值分析程序。。。。。好多程序下载 喜欢的朋友可以下 随时欢迎」· C++ 代码 · 共 34 行

#include<iostream.h>
#include<math.h>
double f(double x)
{     double y;
      y=exp(-x);
	  return  y;
}
void main()
{    double x0,x1,x2,z;
     int N;
	 int k=1;
     cout<<"请输入初值x0, 误差限z,最大迭代次数N的值"<<endl;
	 cin>>x0>>z>>N;
   L:
	 x1=f(x0);
	 x2=f(x1);
	 x2=x2-(x2-x1)*(x2-x1)/(x2-2*x1+x0);
	 if(fabs(x2-x0)>=z)
	 {
		 if(k!=N)
		 {
			 k=k+1;
			 x0=x2;
			 goto L;
		 }
		 else
			 cout<<"迭代失败"<<endl;

	 }
	 else
		 cout<<x2<<endl;

}