矩阵求逆,求解行列式的起源代码,关于开发矩阵类很有帮助的文件
#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