矩阵代码大全(矩阵的逆、乘法、加法、行列式)(c++程序)(5)

矩阵的相关操作

double fun(double array[n][n]);

double a[n][n],b[n][2*n],c[n][n],det1,yinzhi;

double bb;

int i,j,kk=0,k,u;

for(i=0;i<n;i++)

for(j=0;j<2*n;j++)

b[i][j]=0;

cout<<"请输入一个方阵"<<endl;

for(i=0;i<n;i++)

for(j=0;j<n;j++)

cin>>a[i][j];

for(i=0;i<n;i++)

for(j=0;j<n;j++)

b[i][j]=a[i][j];

for(j=0;j<n;j++)

b[j][n+j]=1;

// cout<<"a所对应的at矩阵b为："<<endl;

// for(i=0;i<n;i++)

// for(j=0;j<2*n;j++)

// { cout<<setw(6)<<b[i][j];

// kk=kk+1;

// if(kk%(2*n)==0) cout<<endl;

// }

det1=fun(a);

for(i=0;i<n;i++)

{ if(b[i][i]==0)

for(j=i;j<n;j++)

{ if(b[j][i]!=0)

temp(b[i],b[j],n);

}

for(k=i+1;k<n;k++)

{ yinzhi=-1*b[k][i]/b[i][i];

for(u=0;u<2*n;u++)

{ b[k][u]=b[k][u]+b[i][u]*yinzhi;

}

}

}

det1=fun(a);

if(det1==0)

cout<<"此矩阵不可逆："<<endl;

if(det1!=0)

{

for(i=0;i<n;i++)

{bb=b[i][i];