Метод релаксаций

void SUR(double **a,int n,double w,double *f,double *x_app/*,double eps*/)
{
/*  for(int i = 0;i < n;i++)
    {
        x_app[i] = 0;
    }*/
    //-------------------------------------------
    double * fut = new double[n];
    //-------------------------------------------
    for(int i = 0;i < n;i++)
    {
        fut[i] = x_app[i];
    }
    //-------------------------------------------
 
    bool stop = true;
    for(int k = 0;stop;k++)
    {
        for(int i = 0;i < n;i++)
        {
            double s = 0;
            for(int j = 0;j < i;j++)
            {
                s += a[i][j]*fut[j];
            }
 
            for(int j = i + 1;j < n;j++)
            {
                s += a[i][j]*x_app[j];
            }
 
            double t = 1./a[i][i];
            fut[i] = (1. - w)*x_app[i] + w*t*(f[i] - s);
        }
 
        double temp = fabs(fut[0] - x_app[0]);
 
        for(int i = 1;i < n;i++)
        {
             if(fabs(fut[i] - x_app[i]) > temp) temp = fabs(fut[i] - x_app[i]); 
        }
 
        if(temp < eps)
        {
            stop = false;
        }
 
        for(int i = 0;i <n;i++)
        {
            x_app[i] = fut[i];
        }
        //Сюда следует добавить условие останова
 
        std::cout << std::endl;
    }   
}
// x_app - первое приближение
// f - правая часть
// w - параметр
// a - Основная матрица СЛАУ