矩阵求逆,求解行列式的起源代码,关于开发矩阵类很有帮助的文件

#include

#include

#define MC NR

#define NR 5 //只需在此修改矩阵的行数，然后准备好矩阵文件，结果在weng_out.txt文件中

void ArrayOut(long double *p,int m,int n);

void MatDiv(long double *b,int m,long double *c);

void main()

{

long double A[NR][MC]={0.0},AT[MC][NR]={0.0},C[NR][MC]={0.0};

int i=0,j=0;

FILE *fp=NULL;

//----------------------------------读取文件--------------------

fp=fopen("weng.txt","r+");

for(i=0;i

fscanf(fp,"%lf",(*A+i) );

if(fgetc(fp)==EOF) break;

}

fclose(fp);

printf("读取的文件：\n"); ArrayOut(*A,NR,MC); printf("---------------------\n"); MatDiv(*A,NR,*C); ArrayOut(*C,NR,MC); //------------------------------------------------------------------- fp=fopen("weng_out.txt","w+"); fprintf(fp,"\n\n---------result 矩阵的逆 ，精确结果10-8以上------------\n"); for(i=0;i