метод Зейделя. c++

[Онечка]

В цикле do while есть ошибка, из-за неё программа не правильно считает Xold и Xnew. Помогите, пожалуйста, найти ошибку...

#pragma once

#include <iostream>//Подключение библиотеки стандартного ввода-вывода;
#include<conio.h>//Подключение библиотеки консольного ввода-вывода;
#include<iostream>//Подключение библиотеки потокового ввода-вывода;
#include<math.h>//Подключение библиотеки математических функций;

using namespace std;
void main()
{
	int i,j,m,n;
	double a[10][10],c[10][10],b[10],d[10],Xold[3], Xnew[3], eps, max;

	eps=0.00001;

	
 cout<<endl<<"Vvedite matricu A postrochno"<<endl;
  for (i=0; i<=2; i++)     
	  for (j=0; j<=2; j++)
	  {
	    cin>>a[i][j];
		  while ((i==j)&&(a[i][i]==0.0))
     	{		
	    	cout<<"Error! Enter a[i][i] again!";
		    cin>>a[i][j];
	    }
	  }

        
 
cout<<endl<<"Vasha matrica"<<endl;
   for (i=0; i<=2; i++) 
      { cout << "\n stroka" << (i+1) << ":";
       for (j=0; j<=2; j++)
       cout<<"\t"<<a[i][j];
      }

  
   
    cout<<endl<<"Vvedite matricu B postrochno"<<endl;
   for (int i=0; i<3; i++)
	   cin>>b[i];

   cout<<endl<<"Vash stolbec"<<endl;
   for (i=0; i<3; i++) 
             cout<<"\t"<<b[i]<<endl;


  
   
   for(i=0;i<=2;i++)
   {
	   for(j=0;j<=2;j++)
	   {
		   c[i][j]=-a[i][j]/a[i][i];
		   if(i==j)c[i][j]=0;
		  cout<<"c="<<c[i][j]<<endl;
		  
	   }
	      
	    d[i]=b[i]/a[i][i];
        cout<<"d="<<d[i]<<endl;
	   
   }

   Xnew[0]=1.1;   
   for(i=1;i<3;i++)Xnew[i]=0.9;
   
do 
{
  for (i=0;i<3;i++)Xold[i]=Xnew[i];
  Xnew[0]=0.0;
  for (j=0;j<3;j++)Xnew[0]+=c[0][j]*Xold[j];
  max=fabs(Xnew[0]-Xold[0]);
  for (i=1;i<3;i++)
  {
      Xnew[i]=0.0;
	  for (j=0;j<=(i-1);j++)Xnew[i]+=c[i][j]*Xnew[j];
      for (j=i;j<3;j++)Xnew[i]+=c[i][j]*Xold[j];
	  if( max<fabs(Xnew[i]-Xold[i]))
		  max=fabs(Xnew[i]-Xold[i]);
  } 
  

}
while(max>eps);

for (i=0;i<3;i++)cout<<"Xold="<<Xold[i]<<endl;
cout<<"max="<<max<<endl;

getch ();
}