Метод зейделя

[Praud]

Ребята, помогите найти программу решения систем алгебраических уравнений методом Зейделя. Код нужен на С++ или на С#. Переделывать с паскаля, с джавы или еще с чего-то нет не времени не сил. Нашел десяток разных реализаций и чуть ли не каждая выдавала свой уникальный ответ на разные входные данные. При этом в интернет онлайн-калькуляторе выдавался еще один уникальный ответ. Я иссяк уже

[Bugrimov]

Так вы не ищите готовую реализацию, а найдите алгоритм и сами его реализуйте... Вот тогда будет самый правильный ответ....
Шикарная ссылка: http://shpora.gixx.ru/shpory/%D0%9C%...0%BB%D1%8F.htm

[Praud]

const int n = 3;

using namespace std;

int main()
{
   double A[n][n]; // входная матрца
   double b[n]; // свободные члены
   double x[n]; // решение
   double eps; // требуемая точность вычисления
   double c = 0;
   double p = 0;
   
   cout<<"eps: ";
   cin>> eps;

   cout<<"Matrix:"<<endl;

   for(int i=0;i<n;i++)
	   for(int j=0;j<n;j++)
		   cin>>A[i][j];
   
   cout<<"Koef:"<<endl;
   for(int i=0;i<n;i++)
	   cin>>b[i];

   for(int i = 0; i < n; i++)
   {
		p = 0;
		for(int j = 0; j < n; j++)
		{
			if(i!=j)
				p = p+fabs(A[i][j]/A[i][i]);
		}
		if(c < p)
			c = p;
   }

   for(int i = 0; i < n; i++)
	   x[i] = b[i] / A[i][i];
	
   do{
		p = 0;
		for(int j = 0; j < n; j++)
		{
			c = x[j];
			x[j] = b[j];
			for(int i = 0; i < n; i++)
			{
				if(i!=j)
					x[j]=x[j]-A[j][i]*x[i];
			}
			x[j] = x[j] / A[j][j];
			if(abs(c - x[j]) > (1 - c) / c * eps)
				p+=1;
		}
   }while(p > 1);
	
   cout<<endl<<"Results:"<<endl;
   for(int i = 0; i < n; i++)
	   cout<< x[i]<< endl;

   _getch();
    return 0;
}

Вот по предоставленному выше методу решения попытался написать код и я вишу в бесконечном цикле.
Входные данные брал из того же источника.

[rrrFer]

Praud
На каких данных ты "висишь в бесконечном цикле"?
МБ ты данные будешь с файла вводить - там тестировать проще (не надо каждый раз матрицы руками вводить) и можно файл с матрицей нам скинуть... (мы тоже попробуем).

[Praud]

Интернета не было, не мог ответить.

// входная матрца
	double A[n][n] = {
						{4.2 , 2.3 , 1.3}, 
						{-1.3 , 7.2, 5.4},
						{2.3 , 5.1, 11.2}
					};
   double b[n] = {6.1 , 8.6 , 8.2 };  // свободные члены

[Praud]

Рабочий код нашелся) Помогите написать достаточное условие сходимости метода или может у кого-то имеется. Заранее спасибо.