Metoda Gaussa Seidla

http://pl.wikipedia.org/wiki/Metoda_Gaussa-Seidla

Tutaj we wzorze mam x_j do k+1, ale przecież w 0 iteracji ten wektor nie istnieje to znaczy, że mam zaimplementować dwa warianty?

w którym?

 //Iteracje
for(int iter=0;iter<2;iter++)
{
for(i=0;i<W;i++)
{
for(j=0;j<K;j++)
{
double s=0;
s+=A[i][j]*X[j];
if(i!=j)XN[i]=(B[i]-s)/A[i][i];
}
cout<<"X:"<<i<<"="<<XN[i]<<endl;
}
}

Ten kawałek kodu dobrze przedstawia implementacje numeryczne Metody Gaussa Seidla z wiki? http://pl.wikipedia.org/wiki/Metoda_Gaussa-Seidla