C语言中的矩阵演算:行列式演算
#include
#include
#define d 2
double det (double A[d][d], int N)
{
double c, r=1跨度>;
for(int i = 0跨度>;我
for(int k = i+1; k
c = A[k][i] / A[i][i];
for(int j = i; j
A[k][j]= A[k][j] - c*A[i][j];
}
}
for (int i = 0; i
r *=A[i][i];
返回 r;
}
int main()
{
双 M[d][d];
M[0][0]=9跨度>;
M[0][1]=5跨度>;
M[1][0]=3跨度>;
M[1][1]=7跨度>;
printf(Det(M) = %f\n”,det(M,d));
返回 0;
}