6阶矩阵求逆 c语言,矩阵的LU分解以及求逆矩阵--c语言

matrix.h

#ifndef MATRIX_H_

#define MATRIX_H_

#include

#include

#include

#include

#include

#include

#include

#ifdef N

#define PRINTMATRIX(A, S) \

for (int i = 0; i < N; i++) \

{ \

for (int j = 0; j < N; j++) \

{ \

printf("%lf\t", A[i][j]); \

} \

printf("\n"); \

} \

printf("=================================\n");

bool floatCompare(double f1, double f2)

{

return fabs(f1 - f2) < 1.0e-05;

}

//求对角矩阵的det

double DiagDet(double Matrix[N][N])

{

double res = Matrix[0][0];

for (int i = 0; i < N; i++)

{

res = res * Matrix[i][i];

}

return res;

}

void MatMul(double lhs[N][N], double rhs[N][N], double res[N][N])

{

for (int row = 0; row < N; row++)

{

for (int col = 0; col < N; col++)

{

double ret &#