简单迭代法的代码实现.doc

《简单迭代法的代码实现.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、/*简单迭代法的代码实现*/#include#include#includeusingnamespacestd;doublee=2.718281818284;doublef(doublex){returnpow(e,-1*x);}voidSimpleDiedai(doublex,doubled){doublea=x;doubleb=f(a);intk=0;//记录循环的次数while(((a-b)>d)

2、

3、((a-b)<-1*d)){cout<

4、++;if(k>100){cout<<"迭代失败！（可能是函数不收敛）"<>x>>d;SimpleDiedai(x,d);return0;}/*牛顿迭代法的代码实现*/#include#include#includeusingnamespacestd;doublee=2.718281818284;doublef(

5、doublex){doublea=pow(e,-1*x);returnx-(x-a)/(1+a);}voidNewtonDiedai(doublex,doubled){doublea=x;doubleb=f(a);intk=0;//记录循环的次数while(((a-b)>d)

6、

7、((a-b)<-1*d)){cout<100){cout<<"迭代失败！（可能是函数不收敛）"<

8、"请输入初始值x0和要求得结果的精度：";doublex,d;cin>>x>>d;NewtonDiedai(x,d);return0;}/*雅可比算法的代码实现*/#include#include#include#includeusingnamespacestd;//函数求数组中的最大值doubleMaxOfList(vectorx){doublemax=x[0];intn=x.size();for(inti=0;ima

9、x)max=x[i];returnmax;}//雅可比迭代公式voidJacobi(vector>A,vectorB,intn){vectorX(n,0);vectorY(n,0);vectorD(n,0);intk=0;//记录循环次数do{X=Y;for(inti=0;i

10、i][i];cout<100){cout<<"迭代失败！（可能是函数不收敛）"<0.00001

11、

12、MaxOfList(D)<-0.00001);return;}intmain(){intn;cout<<"请输入方程组未知数的个数n：";cin>>n;cout<

13、le>>A(n,vector(n,0));vectorB(n,0);cout<<"请输入方程组的系数矩阵："<>A[i][j];}}cout<>B[k];}cout<

14、b++){cout<