高斯賽德爾迭代算法C_第1頁
高斯賽德爾迭代算法C_第2頁
高斯賽德爾迭代算法C_第3頁
高斯賽德爾迭代算法C_第4頁
全文預(yù)覽已結(jié)束

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)

文檔簡介

1、 高斯賽德爾迭代算法C+實現(xiàn)依據(jù)公式 (i=1,2,n)#include "stdafx.h"#include "iostream"#include "iomanip"using namespace std;int main()int i,j,k=0,m,n;double t1,t2,e1,e2=0.0; cout<<"請輸入精度e:"cin>>e1;cout<<"請輸入系數(shù)矩陣行數(shù):"cin>>m;cout<<"請輸入系數(shù)矩

2、陣列數(shù):"cin>>n;cout<<endl;double (*a)=new double *m;/生成二維動態(tài)數(shù)組for(i=0;i<=m;i+)ai=new doublen;double (*b)=new double m;double (*x)=new double n;cout<<"請輸入系數(shù)矩陣:"<<endl;cout<<"-"<<endl;for(int num1=0;num1<m;num1+)for(int num2=0;num2<n;nu

3、m2+)cin>>anum1num2;cout<<endl;cout<<"輸入的系數(shù)矩陣為:"<<endl;for (int num3=0;num3<m;num3+)for(int num4=0;num4<n;num4+)cout<<anum3num4<<" "cout<<endl;cout<<"-"<<endl;cout<<"請輸入矩陣b:"<<endl;cout<

4、<"-"<<endl;for(int num5=0;num5<m;num5+)cin>>bnum5;cout<<"輸入的矩陣b為:"<<endl;for(int num6=0;num6<m;num6+)cout<<bnum6<<" " cout<<endl; cout<<"-"<<endl;for(int num7=0;num7<n;num7+)xnum7=0.0000;do cout

5、<<"第"<<k<<"次迭代值:"e2=0.0;for(i=0;i<m;i+) double sum=0.0;for(j=0;j<n;j+) if(j!=i)sum+=aij*xj;t1=xi;t2=e2;xi=(bi-sum)/aii;e2=(xi)-t1>=0?(xi)-t1:t1-(xi);e2=(e2>=t2?e2:t2);cout<<setprecision(8)<<xi<<" "cout<<endl; k+; while(e2>=e1&&k<30);cout<<

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論