gauss.cpp

用Gauss法来直接解方程组的c程序

/***************************************************************
 **********      Gauss法直接解线性方程组      *****************
 **************************************************************/
#include "stdio.h"
#define n 12
void main()
{
  FILE *fp;
  double a[n][n+1],s;
  int i,j,k;

 /*****  打开t1.txt文件,并读入数据 *****/
  fp=fopen("D:\\t1.txt","r");      
  if(!fp) {printf("file t1.txt not found.\n");return;}
  for(i=0;i<n;i++) for(j=0;j<n+1;j++) fscanf(fp,"%lf",&a[i][j]);
  fclose(fp);

  /***** 开始进行正向消元 *****/
  for(k=0;k<n-1;k++)
  { for(i=k+1;i<n;i++)
	 {
	  s=-a[i][k]/a[k][k];a[i][k]=0;
	  for(j=k+1;j<n+1;j++) a[i][j]+=s*a[k][j];
	 }
  }

  /*****  反向迭代的过程  *****/
  for(k=n-1;k>=0;k--)
  {
    for(j=k+1;j<n;j++) 
	 a[k][n]-=a[k][j]*a[j][n];
	 a[k][n]/=a[k][k];
  
  }

  /*****  输出所得结果 *****/
  for(k=0;k<n;k++)
  {	printf("%12lf\n",a[k][n]);}
}