Résolution Gauss Seidel, code qui foire...

[invite1301f9c9]

Bonjour, alors voilà je code actuellement Gauss Seidel, et je ne comprend pas ce qui ne va pas dans mon code. Mon vecteurX0 devrait aller en convergeant, et ce n'est pas le cas....
Je sais plus trop quoi faire.

  do
    {
        for(i = 0; i<dimension; i++)
        {
            sommeA = 0;
            sommeB = 0;
            for(j = i; j<dimension; j++)
            {
                sommeA += Matrice_A2[i][j]*VecteurX1[j];
            }
            for(j = 0; j<i; j++)
            {
                sommeB += Matrice_A2[i][j]*VecteurX2[j];
            }

            VecteurX2[i] = (Matrice_B[i]-sommeA-sommeB)/Matrice_A2[i][i];
        }

        comp_max = 0;
        for(i = 0; i<dimension; i++)
        {
             VecteurX0[i] = fabs((VecteurX2[i] - VecteurX1[i]));
             VecteurX1[i] = VecteurX2[i];
        }
        comp_max = 0;
        for(k = 0; k < dimension ; k++)
        {

            if(VecteurX0[k] > comp_max)
            {
                  comp_max = VecteurX0[i];
            }
        }
    }while(comp_max < crit_arret);

-----