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

矩阵的相关操作

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

b[i][j]=b[i][j]/bb;

}

for(i=n-1;i>0;i--)

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

{bb=b[k][i];

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

b[k][u]=b[k][u]-bb*b[i][u];

}

}

// cout<<"变化后的at矩阵"<<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;

// }

// cout<<endl;

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

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

c[i][j]=b[i][j+n];

kk=0;

if(det1!=0)

{ cout<<"其可逆且其行列式的值det为："<<det1<<endl;

cout<<"a矩阵的逆矩阵为c矩阵:"<<endl;

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

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

{ cout<<setw(15)<<c[i][j];

kk=kk+1;

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

}

}

return 0;

}

void temp(double aa[],double bb[],int n)

{ int i;

double temp1;

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

{temp1=aa[i];aa[i]=bb[i];bb[i]=temp1; }

}

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

{ int ii,jj,k,u;

double det1=1,yin;