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;
}