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

運用列主元消去法編程

發布時間: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]即為解集.

閱讀全文

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

熱點內容
電視端APP如何開會員劃算 瀏覽:887
布爾運算的命令 瀏覽:466
蝦皮app商家版怎麼翻譯成中文 瀏覽:550
移動網加密不可上網 瀏覽:698
工作組伺服器有什麼特點 瀏覽:869
求算術平方根演算法設計思想 瀏覽:363
核彈打擊加密文件 瀏覽:820
虛擬機安卓9怎麼安裝toos 瀏覽:578
安卓怎麼下載joox 瀏覽:912
pdf天空 瀏覽:763
php讀取二維數組 瀏覽:348
php編譯安裝參數 瀏覽:278
其實壓力沒那麼大程序員圖片 瀏覽:416
如何查看app內訪問的網頁地址 瀏覽:757
安卓手機信號旁邊的漢字怎麼設置 瀏覽:304
nrf2401單片機 瀏覽:713
清除電腦文件夾垃圾的方法 瀏覽:227
天河程序員 瀏覽:192
成都程序員公積金 瀏覽:768
程序員為什麼叫程序猿 瀏覽:484