Метод простых итераций (для трехмерной матрицы) - C++

// строка 32
    for (i = 0; i < n; i++)
        for (j = 0; j < n; j++)
        {
            if (i = j)
                al[i][j] = 0;
            else
                al[i][j] = (-1) * (a[i][j] / a[i][i]);
        }
 
        for (i = 0; i < n; i++)
        {
            bt[i] = b[i] / a[i][i];
            x0[i] = bt[i];
        }
 
        for (i = 0; i < n; i++)
            for (j = 0; j < n; j++)
            {
                x1[i] = bt[i] + (al[i][j] * x0[i]);
                x0[i] = x1[i];
            }
            
        for (i = 0; i < n; i++)
            k = k + sqrt(x1[i] - x0[i]);