不动迭代算法.cpp

来自「不动迭代算法 计算方法」· C++ 代码 · 共 37 行

#include<iostream.h>
#include<math.h>
const int K=10;
void main()
{
	//double absmm(0.0000);
	 
    double min=0.01;
	double mm;
	double x[K];
	for(int i=0;i<K;i++)
	{
		x[i]=0;
	}
	
    x[0]=1.5;
	for(int k=1;k<=K;k++)/////x1=f(x0)
	{
		x[1]=sqrt(10-x[0]*x[0]*x[0])/2;
	
       // cout<<x[1]<<endl;
    if(x[1]>x[0])
	   mm=x[1]-x[0];
	else
	   mm=x[0]-x[1];
//	cout<<mm<<endl;
	//absmm=abs(mm);
	if(mm<=min)
	{
		cout<<"x[1]="<<x[1]<<"\tk="<<k<<endl;
		k=K+1;
	}
	else 
		x[0]=x[1];
//cout<<x[0];
	}
}