0

首先,我是 ci 的新手,尝试使用此代码获取 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 */
4

1 回答 1

0

您在 printf 中有错误的格式说明符。

尝试:printf("%f",a);printf("%f",result);

请参阅此处的 printf 格式说明符http://www.cplusplus.com/reference/cstdio/printf/(%d 用于整数而不是双精度数)。

于 2013-11-27T01:25:24.803 回答