高斯消元法

#include "stdafx.h"
#include<stdio.h>
#include<math.h>
int main()
{
   double m,p,A[10][10],X[10];
   int n,i,j,k,q;
   char c;
   printf("请输入方程的阶数（小于等于8）：\n");
   scanf("%d",&n);
   for(i=1;i<=n;i++)//方便起见从A[1][1]开始存入数据
   {
     printf("请输入增广矩阵第%d行：\n",i);
     for(j=1;j<=n+1;j++)
     scanf("%lf",&A[i][j]);
   }//for
   for(i=1;i<n;i++)
   {
      m=fabs(A[i][i]);j=i;   //m表示矩阵第i列中的最大值，j用来标记最大值所在的行数
      for(k=i+1;k<=n;k++)
      {
        if(fabs(A[k][i])>m)
        {m=fabs(A[k][i]);j=k;}//易错
      }//for
      for(q=i;q<=n+1;q++)//交换第i行和第j行
      {
        p=A[i][q];
        A[i][q]=A[j][q];
        A[j][q]=p;
      }//for
      //*************************消去过程*********
      for(k=i+1;k<=n;k++)
      {
         m=-A[k][i]/A[i][i];
         for(j=i;j<=n+1;j++)
         A[k][j]=A[k][j]+A[i][j]*m;//易错
      }//for
       }//for
      //*************回代过程*************
      X[n]=A[n][n+1]/A[n][n];
      for(i=n-1;i>=1;i--)
      {
        p=0.0;
        for(j=i+1;j<=n;j++)
        p=p+A[i][j]*X[j];
        X[i]=(A[i][n+1]-p)/A[i][i];
      }//for
   printf("线性方程组的解为：\n");
   for(i=1;i<=n;i++)
   printf("x[%d]=%lf\n",i,X[i]);
   c=getchar();
   return 0;
}