解n元线性方程组的C语言实现

#include<stdio.h>
int main ()
{
    int n,i,j,p,a,b,q;

    printf("Enter value n:\t");//输入阶数n// 
	scanf("%d",&n);
    double A[n][n],B[n-1][n-1],C[n][n],D[n][n];//A为原始数组，B为A的余子式，C为A的伴随矩阵（未转置），D为A的伴随矩阵A//
    double X[n],Q[n];//X    ,Q   //
	double det(double*p,int n);

    printf("Enter the array A:\n");//读入数组A// 
	for(i=0;i<n;i++)
	{
		for(j=0;j<n;j++)
			scanf("%lf",&A[i][j]);
	}
    printf("Enter the array A:\n");//读入数组Q// 
	for(i=0;i<n;i++)
	{
			scanf("%lf",&Q[i]);
	}
	
    printf("\n\nThe array A you enter:\n");//print original array/
	for(i=0;i<n;i++)
	{
		for(j=0;j<n;j++)
			printf("%10lf\t",A[i][j]);
		printf("\n\n");
	}
	 printf("\n\nThe array Q you enter:\n");//print original array/
	for(i=0;i<n;i++)
	{
			printf("\t%10lf\n",Q[i]);
		printf("\n\n");
	}
	
	
	
//(n-1)array to a new array//
    for(i=0;i<n;i++)
    {
	    for(j=0;j<n;j++)
	    {
		    for(p=0,a=0;p<n;p++)
		    {
			    if(p!=i)
			    { 
				    for(q=0,b=0;q<n;q++)
				    {
					    if(q!=j)
					    {
						    B[a][b]=A[p][q];

							if((i+j)%2==0)
								C[i][j]=det(B,n-1);//forget *(-1)//
							else
								C[i][j]=-det(B,n-1);

						    b++;
					    }
					
				    }
				    a++; 
			    }	
		    }	
	    }	
    }// end//


double result;
result=det(A,n);


if(result==0)
{
	printf("The array X is not exist.\n\n\n");
}

else                                        //error?//
{

    for(i=0;i<n;i++)
    {
	    for(j=0;j<n;j++)
		    D[j][i]=C[i][j]/result;
    }	

//计算A逆(D)*Q==X//
    for(i=0;i<n;i++)
    {
        X[i]=0;
        for(j=0;j<n;j++)
        {
            X[i]+=D[i][j]*Q[j];
        }
    }

    printf("\n\nThe array X :\n");
	for(i=0;i<n;i++)
	{
			printf("%10lf\t",X[i]);
	}
	
}
    return 0;
}



//det A//
double det(double*p,int n)
{

	int x,k,i,j;
	double t,h,det;
	double A[n][n],C[n][n];

for(i=0;i<n;i++)
	{
		for(j=0;j<n;j++)
            A[i][j]= *(p+i+n*j);
	}

	//核心代码// 
	for(j=0;j<n-1;j++)//第一次化归//
	{ 
		
		for(i=n-1,t=0;i>=j;i--)
		{
			if(A[i][j]!=0) 
			{
				t=A[i][j];
				x=i;
			}
		}
		if(t!=0)
		{
		for(h=0,i=n-1;i>=j+1;i--)
		{
			if(i!=x)
			{   h=A[i][j];
				for(k=0;k<n;k++)
				{
					A[i][k]-=h*A[x][k]/t;//error?// 
				}
			}
		}
		}
	}
	
	
	for(i=0;i<n;i++)//旋转数组// 
	{
		for(j=0;j<n;j++)
		{
			C[i][j]=A[n-j-1][n-i-1];
		}
	}
	for(j=0;j<n-1;j++)//第2次化归//
	{ 
		
		for(i=n-1,t=0;i>=j;i--)
		{
			if(C[i][j]!=0) 
			{
				t=C[i][j];
				x=i;
			}
		}
		if(t!=0)
		{
		for(i=n-1;i>=j+1;i--)
		{
			h=C[i][j];
			if(i!=x)
			{
				for(k=0;k<n;k++)
				{
					C[i][k]-=h*C[x][k]/t;//error// 
				}
			}
		}
		}
		
	}
	for(i=0;i<n;i++)//旋转回去// 
	{
		for(j=0;j<n;j++)
		{
			A[i][j]=C[n-j-1][n-i-1];
		}
	}
	for(det=1,i=0;i<n;i++)
	{
		det*=A[i][i];
	}
	return (det);
		//end det//
}