矩阵函数

******************************/

void plus(double *A, double *B, double *C, int m, int n, int R) { int i,j,k; for(i=0;i/****************************** 矩阵转置

矩阵a 为m*n的矩阵， 返回矩阵a_t n*m

******************************/

void transpose(double *a, double*a_t,int m,int n) { int i,j; for(i=0;i*(a_t+i+j*m)=*(a+j+i*n); } return; }

/****************************** 矩阵求逆

矩阵a 为m*n的矩阵， 返回矩阵e 为a的逆矩阵

矩阵a可逆的充要条件是a是方阵且a为非奇异矩阵 ******************************/

void inverse(double *a,double*e,int m,int n) { int i,j; double temp; double temp2; /*若该矩阵不是方阵，则不能求逆*/ if(m!=n) { cout<<\"a不是方阵，不能求逆！\"<{ temp2=*(a+i+k*m); for(j=0;j0;i--) { for(j=0;j=0;k--) { temp2=*(a+i+k*m); for(j=0;j