#include <iostream.h>
#define DIM 200
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#ifndef EPSILON
#define EPSILON 0.000000001
#endif
int gelim( int process, double A[DIM][DIM], double xx[DIM] )
{
int k,i,j;
if( process == 1 ) printf("The process of elimination\n");
/* elimination step */
for(k=0;k<DIM;k++)
  {
  if( fabs(A[k][k]) < EPSILON ) return(1);
  for(i=k+1;i<DIM;i++)
    {
    A[i][k]=A[i][k]/A[k][k];
    for(j=k+1;j<DIM;j++) A[i][j]=A[i][j]-A[i][k]*A[k][j];
    xx[i]=xx[i]-A[i][k]*xx[k];
    }
  if(process == 1 )
    {
    for(i=0;i<DIM;i++)
      {
      for(j=0;j<DIM;j++) printf("%10.6f",A[i][j]);
      printf("   | %10.6f\n",xx[i]);
      }
    printf("\n");
    }
  }
/* backward step */
xx[DIM-1]=xx[DIM-1]/A[DIM-1][DIM-1];
for(i=DIM-2;i>=0;i--)
   {
   for(j=i+1;j<DIM;j++) xx[i]=xx[i]-A[i][j]*xx[j];
   xx[i]=xx[i]/A[i][i];
   }
return(0);
}
int main()
{
 float a[4][3]={0,2,1,3,1,0,1,3,2,1,5,0};
 float b[4]={40,100,100,80};
 gelim(0,a,b);
 for(i=0;i<4;i++)
  cout<<b[i]<<endl;
}

评论