ЭТО САМ МЕТОД ЗДЕСЬ ТАК?

for (int j = i + 1; j<K + 1; j++)
        {
            c = -sum[j][i] / sum[i][i];
            cout<<"\n\n ! "<<c<<endl;
            for (int k = i; k<K + 1; k++)
            {
                sum[j][k] = sum[j][k] + c*sum[i][k];
            }
            b[j] = b[j] + c*b[i];
        }
    }

    for (int i = 0; i<K + 1; i++)
    {
        for (int j = 0; j<K + 1; j++)
        {
            cout << sum[i][j] << setw(15);
        }
        cout << "b" << i << "=";
        cout << b[i] << endl;
    }
    cin.get();

    a[K] = b[K] / sum[K][K];

    for (int i = K + 1 - 2; i >= 0; i--)
    {
        for (int k = i + 1; k<K + 1; k++)
        {
            b[i] = b[i] - a[k] * sum[i][k];
        }
        a[i] = b[i] / sum[i][i];
    }


    for (int i = 0; i<K + 1; i++)
        cout << "X" << i << "=" << a[i] << endl;
    cin.get();