Ajuda: não consigo fazer Gauss Seidel (c++)

[User F]

Estou tentando fazer um programa em gauss seidel pra faculdade e não esta funcionando, gostaria que alguém me desse um auxilio...

a matriz A ficaria

A00 A01 A02 x0 B0

A10 A11 A12 * x1 = B1

A20 A21 A22 x2 B2

aí ficaria considerando x0=0 , x1=0 e x2=0 pelo que o usuario definirá ao ser perguntado

com o valor encontrado de x0 encontrado considerando x1=0 e x2 =0, eu iria jogar pra achar x1 da segunda linha, e com o valor encontrado de x0 e x1 eu iria encontrar x 2 da terceira linha e assim sucessivamente por uma quantidade de interações (repetições) até chegar na aproximação destes...

O programa tá dando looping, gostaria de saber qual o problema.

a00 * x0 + a01*x1 + a02*x2 = b0

x0 = ( b0 - a01*x1 - a02*x2 )/ x0

#include<iostream>
#define N 3 // tamanho da matriz
using namespace std;
int main(){
   double A[N][N]; // MAtriz
   double B[N]; // Resultado
   double X[N];  // coluna X
   int max;
   cout<<"\n Entre com os valores da Matriz A \n";
   for(int a=0;a<N;a++){
           for(int b=0;b<N;b++)
           cin >> A[a][b];
   }
   cout <<"\nDigite os valores da Matriz B \n ";
   for(int c=0;c<N;c++)
           cin>>B[c];
           
   cout <<"\n qual o numero de interacoes q você deseja? ";
   cin >> max;
   cout << "\n Digite os valores iniciais de X ";
   cin >> X[N]; //tem que ser iguais a zero

 //parte que calcula numa matriz 3 por 3 e imprime

 /*  for(int k=0;k<max;k++){ //aqui calcula pela matriz 3x3 , mas não se eu mudar o valor de N
            
           X[0]=1/A[0][0] * (B[0]-A[0][1]*X[1]-A[0][2]*X[2]);
           X[1]=1/A[1][1] * (B[1]-A[1][0]*X[0]-A[1][2]*X[2]);
           X[2]=1/A[2][2] * (B[2]-A[2][0]*X[0]-A[2][1]*X[1]); 
         
           cout<<"\n"<<X[0]<<"\t"<<X[1]<<"\t"<<X[2];
           
           }*/
   
  
  
  //parte que deve calcular qualquer tamanho de matriz escolhida:
   
   
   int p=0, w=0, k=0 ,z=0;
   for(int k=0;k<max;k++){
           //vai fazer k iterações até max quantidade
           
   for(int i=0;i<N;i++)
   for(int j=0;j<N;j++){ //varrer a matriz A[i][j] e B[i][j] segundo a operacao
                                                    
                           
           double Soma_e=0,Soma_d=0;      //Soma_e : soma de termos pela esquerda 
                           
                                                           //Soma_d : soma de termos pela direita
                       if(i==j){                                    //somente os a[0][0] , a[1][1] e a[2][2] possuem os X´s interessados
                          
                           do{                        
                           
                           k = k +(j-1);                   //utiliza k e z para não subistiuir i e j...
                           z = z +(i-1);                   //...com os novos valores assumidos na soma
                          
                           Soma_e = Soma_e +  A[i][k]*X[z]; //soma dos valores vao se acumular
                           
                   
                           }while ( 0<k<(N-1) || 0<z<(N-1) );      //o valor de [i] e [j] não podem exceder a matriz...
                                                                   //...como considerar o programa A[-2][-5] dentro da operacao.   
                                                                   //...se N=3, com k=2 vai sair do laço
                           do{ 
                           
                           p = p +(i+1);
                           w = w +(j+1);
                           
                           Soma_d = Soma_d +  A[i][p]*X[w];
                           
                           }while ( 0<p<(N-1) || 0<w<(N-1) );
                           
             
                           
                            X[i] = ( B[i] - Soma_e - Soma_d )/A[i][i];
                            } 
                            else
                            {}             
             cout<<"\n"<<X[0]<<"\t"<<X[1]<<"\t"<<X[2];
             }  
                           
   }
   system("pause>>null");
   return 0;
}

Abraço.

Editado Outubro 30, 2008 por kuroi
Adicionar tag CODE