#include
#include
#include
#define N 3
void main()
{
double a[N][N]={
{1,2,3},{2,4,5},{3,5,6}};//线性方程组的系数矩阵
double b[N]={1,0,0};
double m[N];//辅助乘数
double T,BT=0;//存储换行临时变量
double temp=0;
double t;//最大列主元
int tap1=0,tap2=0;//最大列主元下标
cout<
for(int i=0;i
{
for(int j=0;j
{
cout<
}
cout<
}
cout<
for( i=0;i
{
cout<
}
cout<
for(int s=0;s
{
t=a[s][s];//赋初值
for(int p=s;p
{
if(fabs(a[p][s])>t)
{
t=a[p][s];
tap1=p;
tap2=s;
}
}
//cout<
if(t==0)
{cout<
if(tap1!=tap2)//换行
{
for(int i=s;i
{
T=a[tap1][i];
a[tap1][i]=a[s][i];
a[s][i]=T;
}
BT=b[tap1];//给B换行
b[tap1]=b[s];
b[s]=BT;
}
tap1=tap2=0;//置零
/*cout<
for(int j=s;j
{
a[s][j]=a[s][j]/t;
}
b[s]=b[s]/t;
for(int i=0;i
{
if(i!=s){
m[i]=a[i][s]/a[s][s];
for(int j=s;j
{
a[i][j]=a[i][j]-m[i]*a[s][j];
}
b[i]=b[i]-m[i]*b[s];}
}
/*cout<
for(i=0;i
{
for(int j=0;j
{
cout<
}
cout<
}
cout<
}
cout<
for(i=0;i
{
for(int j=0;j
{
cout<
}
cout<
}
cout<
for( i=0;i
{
cout<
}
cout<
}