高斯消元法求解线性方程组的解集

问题描述：

    运用高斯消元法求解线性方程组的解集，其中线性方程组的表示方式为：

例如：2x+4y=5

         2x+y=2

存储在二维数组中形式为

a[2][3]={{2,4,5,},{2,1,2}}

具体代码实现如下：

#include<malloc.h>
#include<stdio.h>
#include<iostream.h>
void Gauss(double **A,int n)
{
 int i,j,k;
 for(k=0;k<n-1;++k)
 for(i=k+1;i<n;++i)
 for(j=k+1;j<=n;++j)
 A[i][j]-=A[k][j]*A[i][k]/A[k][k];
 for(i=n-1;i>=0;--i)
 {
  for(j=i+1;j<n;++j)
     A[i][n]-=A[i][j]*A[j][n];
     A[i][n]/=A[i][i];//此处A[i][j]改为A[i][i]
 }
}
void main()
{
    int i,j,count;
    printf("请输入维数：");
    scanf("%d",&count);

    //c中建立二维动态数组的方法
    double **A=(double**)malloc(sizeof(double*)*count);
    for (i=0;i<count;i++)
        A[i]=(double*)malloc(sizeof(double)*(count+1));

    //输入线性方程组的参数值，
    for(i=0;i<count;i++)
    {
        printf("请输入第%d行的值:",i+1);
        for(j=0; j<count+1; j++)
             scanf("%lf",&A[i][j]);
    }

    //调用高斯消元公式求解
    Gauss(A,count);

    //打印最终结果（结果保留在每一行的最后一个值中）
    for(i=0;i<count;i++)
         printf("%lf\n",A[i][count]);
    
    //最后释放空间
    for(i=0;i<count;i++)
        free(A[i]);
    free(A);

}

代码中采用了动态数组，其实定义全局数组也可以很好的实现。。。。