my matrix.c

行列式的计算

/* Note:Your choice is C IDE */
#include "stdio.h"
int matrix(int b[40][40],int n)
{
	int temp[40][40];
	int i,j,c,c1; 
    int result=0,t=1; 
    for(i=0;i<n;i++) 
    {//构造临时矩阵，用来计算行列式 
      for(j=0;j<2*n;j++) 
       { 
         temp[i][j]=b[i][j%n]; 
       } 
    } 
    if(n==1) result=b[0][0];
    else
    {
    	if(n==2)
        {
    	  result=b[0][0]*b[1][1]-b[0][1]*b[1][0];
        }
        else 
        {
    	  for(c1=0;c1<n;c1++) 
          {//计算正值 
              i=0; 
              j=c1; 
              t=1; 
              for(c=0;c<n;c++) 
              { 
                 t*=temp[i][j]; 
                 i++; 
                 j++; 
              } 
              result+=t; 
          } 
          for(c1=0;c1<n;c1++) 
          {//计算负值 
             i=n-1; 
             j=c1; 
             t=1; 
             for(c=0;c<n;c++) 
             { 
                 t*=temp[i][j]; 
                 i--; 
                 j++; 
             } 
             result-=t; 
          }
         }
      } 
       return(result); 
     } 
main()
{
    int h=1,i,j,n,r=0,a[40][40];
    char temp;
    printf("XuKai,2008.All rights reserved.\n该工具只能用于1~3阶的行列式或矩阵计算\n按Enter键继续");
    scanf("%c",&temp);
    printf("阶数: ");/*阶数*/
    scanf("%d",&n);
    for(i=0;i<n;i++)
    {
       printf("第 %d 行:\n",i+1);
       for(j=0;j<n;j++)
       scanf("%d",&a[i][j]);
    }
    printf("行列式为:  \n");
    	for(i=0;i<n;i++)
    	{
    		printf("|");
    		for(j=0;j<n;j++)
    		printf("%5d",a[i][j]);
    		printf("|\n");
    	}
    r=matrix(a,n);
    printf("结果为:%d\n         按Enter键退出",r);
    scanf("%c",&temp);
    scanf("%c",&temp);
}