利用Cramer法则求具有唯一解的方程组的解

//求解存在唯一解的方程组

#include <stdio.h>
#include <stdlib.h>
double count(double D[][100],int n);//计算行列式
int toDj(double D[][100],double Dj[][100],double *b,int n,int j);//转换成Dj

int main(void)
{
    double D[100][100]={},Dj[100][100],b[100],x[100],d;
    int n,i,j;//n为未知数个数
    
//    输入未知数数量
    printf("Input the number of unknown digits.\n");
    scanf("%d",&n);
    
//    输入方程
    printf("Enter the equation.\n");
    for (i=0; i<n; i++)
    {
        printf("? x%d",i+1);
        if (i!=n-1)
        {
            printf(" + ");
        }
    }
    printf("=?\n");
    for (i=0; i<n; i++)
    {
        for (j=0; j<n; j++)
        {
            scanf("%lf",&D[i][j]);
        }
        scanf("%lf",&b[i]);
    }
    
    d = count(D, n);//计算D
    
//    D=0时不可以Cramer法则，解不唯一
    if (d == 0)
    {
        printf("More than one solution.\n");
    }
    else
    {
//        D不等于0时，解唯一,利用Cramer计算
        for (j=0; j<n; j++)
        {
            toDj(D, Dj, b, n, j);
            x[j] = count(Dj, n) / d;
            printf("x%d = %.2lf  ",j+1,x[j]);
            
            if (j%5==0 && j!=0 && j!=n-1)
            {
                printf("\n");
            }
        }
        printf("\n");
    }
    
    return 0;
}

//利用第一列展开求行列式
double count(double D[][100],int n)
{
    int i,p,q,s;//i为原行列式行数，p、q为余子式的行列号，s为符号表示
    double M[100][100] = {}, result = 0.0;//M[][]为余子式，result为结果
    
    //    令拆到行列数为0时，余子式为1
    if (n==0)
    {
        return 1;
    }
    
    //    第一列展开
    for (i=0,s=1; i<n; i++,s*=(-1))
    {
        //        求a[i][0]的余子式M[i][0],即删去第i+1行第1列
        for (p=0; p<n-1; p++)
        {
            if (p<i)
            {
                for (q=0; q<n-1; q++)
                {
                    M[p][q] = D[p][q+1];
                }
            }
            else if (p>=i)
            {
                for (q=0; q<n-1; q++)
                {
                    M[p][q] = D[p+1][q+1];
                }
            }
        }
        //        利用递归求解
        result += s*D[i][0]*count(M, n-1);
    }
    
    return result;
}

//将第j+1列与b列交换
int toDj(double D[][100],double Dj[][100],double *b,int n,int j)
{
    int p,q;//p为行，q为列
    
    for (p=0; p<n; p++)
    {
        for (q=0; q<n; q++)
        {
            if (q == j)
            {
                Dj[p][q] = b[p];
            }
            else
            {
                Dj[p][q] = D[p][q];
            }
        }
    }
    
    return 0;
}