ji.c

华中科技大学计算机学院《计算方法》课件（含C源代码）。通过它可以学习求解复杂数学问题的基本方法以及编程实现。

#include <stdlib.h>
#include <stdio.h>
static int seed =10000+1000*1.016;
unsigned int MAX=32767;
unsigned int  M=6,N=7;
static float A[6][7];
void SetMatrixWith(int seed)
{ int i,j,k;
  float x;
  srand(seed);
  for(i=0;i<1000;i++) k=rand();
  for(i=0;i<M;i++)
  for(j=0;j<N;j++)
     { k=rand();
       x=(float)k;
       A[i][j]=x*20/MAX-10;
     }
  return;
}
PrintMatrix()
{ int i,j;
  printf("\n");
  for(i=0;i<M;i++)
  { printf("\n");
    for(j=0;j<N;j++)
    printf("%10.4f",A[i][j]);
  }
  getch();
}
void solution(float A[6][7],int M,int N)
{
 int i,j,k;
 float x;
 for(k=0;k<M;k++)
    { x=A[k][k];
      for(j=k;j<N;j++) A[k][j] /= x;
      for(i=0;i<M;i++)
	 { if(i==k) continue;
	   x=A[i][k];
	   for(j=k;j<N;j++)A[i][j] -=x*A[k][j];
	 }
 }

}
void MatrixOneColumnOutput(float A[6][7],int M,int N)
{
 int i,j;
 for(i=0;i<M;i++)
  { printf("\n");
    for(j=0;j<N;j++)
    printf("%10.4f",A[i][j]);
  }
 //for(i=0;i<M;i++)
 //printf("\nx[%d]=%f",i+1,A[i][N-1]);
}

main()
{ SetMatrixWith(seed);
  PrintMatrix();
  printf("\nOutputsolution\n");
  solution(A,M,N);
  MatrixOneColumnOutput(A,M,N);
  getch();
}