Метод Зейделя

[True_girl]

Ув. программисты! Есть код решающий СЛАУ методом Зейделя, но ответ выдается странный. Что не так? Как исправить? Помогите исправить код.

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <iostream.h>
#define N               3     // size of matrix
#define N1              N+1
float matrix[N][N1]=
      {{4,1,-5,11},
       {3,4,2,1},
       {1,5,5,-10}
      };
	
float z[N]={11.0/4.0,4.0/1.0,5.0/-10.0};    // first iteration
double epsilon=0.001;  // required accuracy
            

/*void GetMatrix (void)
{
	float r;
	for (int i=0; i<N; i++)
		for (int j=0; j<N1; j++)
		{
			scanf("%f", &r);
			matrix[i][j]=r;

		}
}*/
void ShowMatrix(void)
{
  for (int i=0;i<N;i++)
  {
       for (int j=0;j<N;j++)
            printf("%+f*x%d",matrix[i][j],i+1);
       printf("=%f\n",matrix[i][N]);
  }
}

void main(void)
{
  // Variables declaration

  float x[N];
  short int i,j;
  int iterations=0,finish=0;
  // Printing given matrix
//GetMatrix();
  ShowMatrix();
  while (!finish)
  {
     finish=1;
     for (i=0;i<N;i++)
    {
        x[i]=-matrix[i][N];
        for (j=0;j<N;j++) x[i]+=matrix[i][j]*z[j];
       // don't stop iterations until required accuracy is reached
       if (fabs(x[i]/matrix[i][i])>=epsilon) finish=0;
       x[i]=z[i]-x[i]/matrix[i][i];
      // next iteration
      z[i]=x[i];
    }
    iterations++;
	if (iterations==20) printf("20 iter"); 
	if (iterations==100) printf("100 iter");
	if (iterations==2000) printf("2000 iter");
  }
  // Printing solution
  printf("\nSolution:\n");
  for (i=0;i<N;i++)
      printf("x%d=%f\n",i+1,x[i]);
  printf("%d iterations were made",iterations);
}

[True_girl]

ПОМОГИТЕ! ПОЖАЛУЙСТА!