导航:首页 > 编程语言 > 运用列主元消去法编程

运用列主元消去法编程

发布时间: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怎么改密码 浏览:72
魔兽世界表情命令 浏览:985
智能还款信用卡源码 浏览:553
zoo文件夹 浏览:762
安卓2k21如何下载 浏览:648
某年某月的天数python 浏览:912
广度优先算法的复杂度 浏览:399
系统重装网站源码 浏览:152
相册加密相片 浏览:297
美国正常化行政命令 浏览:277
中级审计师教材pdf 浏览:696
wps中pdf旋转 浏览:600
getex命令 浏览:190
云闪付和农行卡app怎么授权 浏览:123
羁绊命令 浏览:51
解压视频怪兽大全 浏览:964
最新中华人民共和国预算法 浏览:861
电信dns服务器地址杭州 浏览:506
java内存泄露内存溢出 浏览:815
笔记本系统映像恢复网络文件夹 浏览:932