G-S迭代 高斯迭代

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define M 10
#define N 100
void main()
{ int i,j,k,n;
  double a[M][M],b[M],x[M]={0},max,XI,eps,s;
  clrscr();
  printf("PLease input n:\n");
  scanf("%d",&n);
  printf("Please input xishu:\n");
  for(i=1;i<=n;i++)
   for(j=1;j<=n;j++)
    scanf("%lf",&a[i][j]);
  printf("please input b:\n");
  for(i=1;i<=n;i++)
    scanf("%lf",&b[i]);
  printf("PLease input eps:\n");
    scanf("%lf",&eps);
  for(k=1;k<=N;k++)
  { max=0;
    for(i=1;i<=n;i++)
    { s=0;
      XI=x[i];
      for(j=1;j<=n;j++)
       if(j!=i)s+=a[i][j]*x[j];
      x[i]=(b[i]-s)/a[i][i];
      printf("The %d time is:",k);
      printf("x[%d]=%lf\n",i,x[i]);