printf("\t\t%4.1f",b[i]);

printf("\n");

    }

printf("\n");

do

    {

Далі іде ітераційний процес:

norma=0.0;

for(i=0;i < N;i++)

        {

xn[i]=-b[i];

for(j=0;j < N;j++)

            {

if(i!=j)

xn[i]+=a[i][j]*x[j];

            }

xn[i]/=-a[i][i];

        }

for(i=0;i < N;i++)

        {

Далі іде обчислення норми вектора:

if(fabs(x[i]-xn[i]) >norma)

norma=fabs(x[i]-xn[i]);