導航:首頁 > 編程語言 > 運用列主元消去法編程

運用列主元消去法編程

發布時間:2024-11-20 01:10:21

Ⅰ 用列主元消去法,用c++編個程序

程序代碼如下:

double* Gauss(double **a,double *b)
{
int i,j,k,flag;
double temp;
double *x;
x=new double[n];
for(i=0;i<n-1;i++)
{
max=0;
flag=0;
for(j=i;j<n;j++)
{
if(a[j][i]<0)
temp=-a[j][i];
else
temp=a[j][i];
if(max<temp)
{
max=temp;
row=j;
flag=1;
}
}
if(max==0)
{
cout<<"此題不能用列主元消去法做"<<endl;
exit(0);
}
if(flag==1)
{
for(j=0;j<n;j++)
{
temp=a[i][j];
a[i][j]=a[row][j];
a[row][j]=temp;
}
temp=b[i];
b[i]=b[row];
b[row]=temp;
}
for(j=i+1;j<n;j++)
{
temp=1.0*a[j][i]/a[i][i];
for(k=i;k<n;k++)
a[j][k]=a[j][k]-temp*a[i][k];
b[j]=b[j]-temp*b[i];
}
if(a[n-1][n-1]==0)
{
cout<<"此題不能用列主元消去法做"<<endl;
exit(0);
}
}
for(i=0;i<n;i++)
{
b[i]=b[i]/a[i][i];
for(j=i+1;j<n;j++)
a[i][j]=a[i][j]/a[i][i];
a[i][i]=1;
}
for(i=n-2;i>=0;i--)
for(j=i;j>=0;j--)
{
b[j]=b[j]-b[i+1]*a[j][i+1];
a[j][i+1]=0;
}
for(i=0;i<n;i++)
x[i]=b[i];
return(x);
}

使用程序時在main函數中調用x=Gauss(a,b),a為系數矩陣頭指針,b為常數項頭指針,程序返回一個double類型的數組x[n]即為解集.

閱讀全文

與運用列主元消去法編程相關的資料

熱點內容
程序員格子襯衫過聖誕 瀏覽:222
插入單片機檢測不到埠 瀏覽:467
svn文件夾如何刪除 瀏覽:621
編譯過程和解析過程 瀏覽:819
德雲天團app如何參與活動 瀏覽:740
編譯系統的硬體原理 瀏覽:330
黑色玫瑰什麼時候出的伺服器 瀏覽:636
如何打開列印共享伺服器 瀏覽:867
手機熱點密碼的加密方式 瀏覽:328
貴州大數據伺服器閑置雲伺服器 瀏覽:527
文件夾加密微軟 瀏覽:383
蘋果突然提示無法驗證伺服器身份怎麼解決 瀏覽:988
調度命令發布的原則 瀏覽:679
dos提示不是內部或外部命令 瀏覽:683
騰訊100g編程源碼 瀏覽:209
機械設計和編程 瀏覽:830
八角單片機紅外遙控led 瀏覽:94
血液凈化pdf 瀏覽:383
烏市停車APP哪裡下載 瀏覽:904
單片機中fg是什麼意思 瀏覽:883