Метод итераций и метод Зейделя

[prikolist]

Это мы обсуждали вопрос решения методом итераций, а методом Зейделя, это как можно осуществить, вот я попробывал так, но ошибка со скобкой.

Вот код, метода Зейделя:

#include "stdafx.h"
#include <iostream>
#include <cmath>
#define eps 0.0001 //требуемая точность вычислений
#define N 4
using namespace std;

//Метод Зейделя
void iteration(double a[N][N],double b[N],double x[N])
{
 int i,j;
 double norma = 0; //чебышевская норма вектора
 double xn[N]={0};//вектор для текущей итерации, начальное значение
       //должно быть равно начальному приближению

 for(i=0;i < N;i++)
  {
   x[i]=-b[i];  
   for(j=0;j < N;j++)
   {
	if(i!=j)   
    x[i]+=a[i][j]*x[j];  
   }
   x[i]/=-a[i][i];  
 } 
 for(i=0;i < N;i++)
}
if(fabs(x[i]-xn[i]) > norma) 
norma=fabs(x[i]-xn[i]); 
xn[i]=x[i];  
}
}

int main()
{

double a[N][N],b[N],x[N];
int i;

a[0][0] = -1.01; a[0][1] = 0.02; a[0][2] = -0.62; a[0][3] = 0.08; b[0] = 1.3;
a[1][0] = -1.03; a[1][1] = 0.28; a[1][2] = 0.33; a[1][3] = -0.07; b[1] = -1.1;
a[2][0] = -1.09; a[2][1] = 0.13; a[2][2] = 0.42; a[2][3] = 0.28; b[2] = 1.7;
a[3][0] = -1.19; a[3][1] = -0.23; a[3][2] = 0.08; a[3][3] = 0.37; b[3] = 1.5;
m1:
for(i=0; i<N; i++)
 x[i] = 0;

iteration(a,b,x);
cout<<"------Metod Zeydela---------"<<endl<<endl;
for(i=0; i<N; i++)
 cout<<x[i]<<endl;

}