#include <iostream>

using namespace std;

void iter(float A[50][50], float(&x)[50][50], int ite, int n)
{
    int suma;
    for (int i = 1; i <= n; i++)
    {
        suma = 0;
        for (int j = 1; j <= n; j++)
        {
            if (i != j)
                suma = suma + ((-A[i][j]) / (A[i][i])) * (x[ite - 1][i]);

        }
        x[ite][i] = (A[i][n + 1] / A[i][i]) + suma;
        cout << "x[" << ite << "][" << i << "] = " << x[ite][i] << endl;
    }
}

int main()
{

    float A[50][50], x[50][50];
    int n, ite;

    cout << "Enter the order " << endl;
    cin >> n;

    cout << "Enter the coeeficients:" << endl << endl;
    for (int i = 1; i <= n; i++)
    {
        for (int j = 1; j <= (n + 1); j++)
        {
            cout << endl << "A[" << i << "][" << j << "] : ";
            cin >> A[i][j];
        }
    }
    system("cls");
    for (int i = 1; i <= n; i++)
    {
        x[0][i] = (A[i][n + 1]) / (A[i][i]);
    }

    for (ite = 1; ite < 10; ite++)
    {
        iter(A, x, ite, n);
    }



    system("PAUSE");
    return 0;
}