c语言---运用高斯消去法解方程

c语言—运用高斯消去法解方程

前段时间学习了高斯消去法，今天才把代码写出，小白第一次发博，若有不妥之处请大家指正。
本算法原理如下（摘自黄健元主编的河海大学出版社的《计算方法》）

#include<stdio.h>
#define n 3
void main()
{
  int k,i,j,q,p；
  double a[n][n+1]={2,2,3,3,4,7,7,1,-2,4,5,-7},m；
  double x[n];
  printf("消元过程\n");
  for(k=0;k<=n-2;k++)
  {
      if(a[k][k]==0)
          break;
      for(i=k+1;i<=n-1;i++)
   {
      m=a[i][k]/a[k][k];
      a[i][k]=0;
      for(j=k+1;j<=n;j++)
     {
      a[i][j]=a[i][j]-m*a[k][j];
     }
   }
  }
  for(q=0;q<n;q++)
  {
   for(p=0;p<=n;p++)
    printf("%.4lf\t",a[q][p]);
   printf("\n");
  }
  printf("回代过程\n");
  x[n-1]=a[n-1][n]/a[n-1][n-1];
  for(i=n-2;i>=0;i--)
  {
    int m;
    double s[n-1]={0};
    for(m=n-1;m>i;m--)
    {
      s[i]=s[i]+x[m]*a[i][m];
    }
    x[i]=(a[i][n]-s[i])/a[i][i];
  }
  for(i=1;i<=n;i++)
  printf("x%d=%lf\n",i,x[i-1]);
}