double det(double*M,int n)//M为矩阵的首地址，你为它的行数（列数）
{
    double a1,a2,add1=0,add2=0,det_M=0;
    int i,j,k;
    if (n==1)
    return det_M=M[0];
    if (n==2)
    return det_M=M[0]*M[3]-M[1]*M[2];
    for(k=0;k<n;k++)
    {
        a1=1;
        for(i=0,j=k;j<n;i++,j++)
            a1*=M[i*n+j];
        if(i!=j)
            for(i=n-1,j=k-1;j>=0;i--,j--)
                a1*=M[i*n+j];
        add1+=a1;
    }
    for(k=0;k<n;k++)
    {
        a2=1;
        for(i=n-k-1,j=0;i>=0;i--,j++)
            a2*=M[i*n+j];
        if((i+j)!=n+1)
            for(i=n-k,j=n-1;i<n;i++,j--)
                a2*=M[i*n+j];
        add2+=a2;

    }
    det_M=add1-add2;
    return det_M;
}
求4阶行列式就开始出错，为什么？该怎样写，求大神么指教。。。

#include<stdio.h>

int oddeven(int i)
{
  if (i&1) return(-1);
  else return(1);
}


int ok(int j,int i,int size)
{
  if (j%size!=i)    return 1; else return 0;
}


double det(double *p,int size ) //work out the determinant
{
  double sum=0;
  int i,j,k;
  double newMatrix[49]={0};
  double *newp=newMatrix;

  if (size==2) {
     return (*p)*(*(p+3))-(*(p+1))*(*(p+2)); }

  if (size==1) {return *p; }

  for (i=0;i<size;i++)
      {
//generate the new matrix

        k=0;
        j=size;

        while (k<=size*size-size-size){

            if (ok(j,i,size))
            newMatrix[k++]=*(p+j);
              j++;
           }



 

//add it to sum by recursive function

          sum+=(*(p+i))*oddeven(i)*det(newp,size-1);
   }

    
    return sum;

}

void main()
{

  double  Matrix[][4]=
  {
   1,2,3,5,
   4,5,6,44,
   7,8,5,41,
   74, 55, 66, 4
  };
  int n,m;
  double *p =Matrix;

   printf("%lf",det(p,4));



}

无聊蛋疼地写了个行列式计算

[ 本帖最后由 zjsxwc 于 2012-5-1 20:46 编辑 ]