#include <iostream>
// 012
#include <iomanip>
// 013
using namespace std;
// 014
int n=10;
// 015

// 016
int main()
// 017
{
// 018
  float tab[3][4]={{0,0,0},
                   {0,0,0},
                   {0,0,0}};

// 019
  float wynik[n];
// 020
  float wspolczynnik,pomoc;
// 021
  int i,j,k,n=3;
// 022
  float suma;
// 023

// 024
 /*Wprowadzenie do tablicy wspolczynnikow*/
// 025
  cout << "\n\tPodaj liczbe rownan w ukladzie:  n = ";
// 026
  cin>>n;
// 027
  cout<<"\n\n";
// 028
  for (i=0 ;i<n ;i++)
// 029
    for (j=0 ;j<n+1 ;j++)
// 030
    {
// 031
      cout << "Podaj tab[" << i+1 << "]" << "[" << j+1 << "] = ";
// 032
      cin >> tab[i][j] ;
// 033
    }
// 034

// 035
cout<<"\n\n";
// 036
system("pause");
// 037
system("cls");
// 038

// 039
 /* Wyswietlanie zawartosci tablicy współczynników*/
// 040

// 041
  cout<<"\n\n\tWyswietlenie zawartosci tablicy wspolczynnikow\n\n";
// 042
  for (i=0; i<n; i++)
// 043
  {
// 044
    cout << endl;
// 045
    for (j=0; j<n+1; j++)
// 046
      cout << setw(5) << tab[i][j];
// 047
  }
// 048

// 049
cout<<"\n\n";
// 050
system("pause");
// 051
system("cls");
// 052

// 053
/* Sprawdzenie czy uklad jest ukladem Cramera*/
// 054
  for (k=0; k<n; k++)
// 055
  {
// 056
    if (tab[k][k]==0)
// 057
    {
// 058
      int z = k+1;
// 059
      while (tab[z][k]==0 && z<n)
// 060
        z++;
// 061
      if (z==n)
// 062
      {
// 063
        cout << "\n\n\tTo nie jest uklad Cramera\n\n";
// 064
        system("pause");
// 065
        return 1;
// 066
      }
// 067

// 068
/* zamiana wierszy miejscami*/
// 069
      else
// 070
        for (j=k; j<n+1; j++)
// 071
        {
// 072
          pomoc = tab[z][j];
// 073
          tab[z][j] = tab[k][j];
// 074
          tab[k][j] = pomoc;
// 075
        }
// 076
    }
// 077
    for (i=k+1; i<n; i++)
// 078
      for (j=n; j>=0; j--)
// 079
        tab[i][j] = tab[i][j]-tab[i][k]*tab[k][j]/tab[k][k];
// 080
  }
// 081
  cout <<"\n\n";
// 082

// 083
  /*sprowadzenie układu do postaci trójkątnej*/
// 084
  cout<<"\n\n\tSprowadzenie ukladu do postaci trojkatnej\n\n";
// 085
  for (i=0; i<n; i++)
// 086
  {
// 087
    cout << endl;
// 088
    for (j=0; j<n+1; j++)
// 089
      cout << setw(5) << tab[i][j];
// 090
  }
// 091
 cout<<"\n\n";
// 092
 system("pause");
// 093
 system("cls");
// 094

// 095
 /*Obliczenie niewidomych ukladu i zapisanie wartosci do tablicy wynikow*/
// 096
  for (i=n-1; i>=0; i--)
// 097
  {
    suma = 0;
    for (j=i+1; j<n; j++){
      suma = suma+tab[i][j]*wynik[j];
      wynik[i] = (tab[i][n]-suma)/tab[i][i];}
  }
  cout << endl;
  cout<<"\n\n\t";
  for(i=0; i<n; i++)
    cout <<"X" << i+1 << " = " << wynik[i]<<"\t";
  cin.ignore();
  getchar();
  return 0;
}