метод релаксации

[xripunov]

Здравствуйте, форумчане!

У меня проблема возникла, требуется реализовать метод релаксации СЛАУ для конкретной системы уравнений. Так вот, прочитал тонны литературы по теме, посмотрел что люди в интернетах пишут, написал прогу. Вроде бы все как в алгоритме, но что то не работает.
Насколько я понял, в процессе работы программы ответы (корни) должны постепенно сходиться к правильным, но этого почему то не происходит.
UPD. Корни изменяются с итерациями, стремятся к значениям, но они не являются корнями => ошибка скорее всего в алгоритме...

	while (convergence==true) // пока сходится
	{for (i=0;i<n;i++)
		{	sigma = 0;
		
		for(j=0;j<i-1;j++)
			{
				sigma+=A[i][j]*xk[j];
			}
			for(j=i;j<n;j++)
			{
				sigma+=A[i][j]*xkp[j];
			}

			xk[i]=xkp[i]+w*(((b[i]-sigma)/A[i][i])-xkp[i]); //конечное приближение (корни) строится на основе начального
		}
			result(A,xk,n);
			system("pause");

                        for(int l=0;l<n;l++)
				xkp[l]=xk[l];
			
			if(converge(b,A,xkp,n,eps)==false)break;
			
	}

Может быть, кто то уже сталкивался с этим методом, буду рад за помощь или подсказку, куда копать.
Спасибо.
Весь код в приложении.