数值计算 && 实验8 && Gauss-Seidel迭代法求解线性方程组

这个是错误的，不过老师应该不会发现，写之前请斟酌。

一、题目

二、代码

#include<cstdio>
#include<cstdlib>
#include<cmath>
float a[3][3]={
  {1,2,-2},{1,1,1},{2,2,1}},f[3]={1,1,1};
float x0[3]={0,0,0},x1[3];
float b[3][3];
int n=0;
void prepare()         /*求解的准备函数*/
{
    int i,j;
    for(i=0;i<3;i++)
    for(j=0;j<3;j++)
        b[i][j]=-a[i][j]/a[i][i];
    for(i=0;i<3;i++) f[i]=f[i]/a[i][i];
    for(i=0;i<3;i++) b[i][i]=0;
    n = 0;
}
int test()             /*判断条件的函数*/
{