首先,我在c.i上尝试使用此代码获取nxn矩阵的决定因素
double Determinant(double *A, int N){
int i;
int j;
int k;
double y[100][100];
double x[100][100];
int sign = 1;
double result = 0;
int arrayRow;
int arrayColumn;
for(i = 0 ; i < N ; i++){
for(j = 0; j < N ;j++){
x[i][j] = A[(N*i)+j];}}
if(N == 2){
result =x[0][0]*x[1][1] - x[0][1]*x[1][0];
}
else
{
for(i = 0; i < N ; i++){
arrayRow = 0;
arrayColumn = 0;
for(j = 1 ;j < N; j++){
for( k = 0; k < N ; k++){
if(k == i){
continue;}
y[arrayRow][arrayColumn] =x[j][k];
arrayColumn++;
}
arrayColumn = 0;
arrayRow++;
}
result += sign*x[0][i]*Determinant(*y,N-1);
sign = -sign;}
}
printf("%d",result);
return result;
} /* end-Determinant */
但它显示了一些非常意外的数字,如-85899 ....这是测试代码,行列式不会返回正确的结果。
int main(){
double A[] = {5, 2, 4, 2, 4, 6, 1, 2, 8};
int a ;
if(Determinant(A, 3) == 80){
a = Determinant(A,3);
printf("%d",a);};
getchar();
} /* end-main */
答案 0 :(得分:0)
printf中的格式说明符错误。
尝试:printf("%f",a);
和printf("%f",result);
请参阅此处的printf格式说明符http://www.cplusplus.com/reference/cstdio/printf/(%d表示不是双精度的整数)。