#include <iostream>
 #include "fstream"
 using namespace std;
 void gauss (double** mas,double* korni, int n,int m)
 {
 for (int i=0;i<m-1;i++)
 {
 double a [n];
 int cj = 0;
 for (int f=0 ; f<n;f++)
 {
 a[f] = mas [n-m][f];
 }
 if (mas [i+n-m][n-m]>a[n-m])
 {
 for (int f=0 ; f<n;f++)
 {
 a[f] = mas [i+n-m][f];
 }
 cj = i;
 }
 for (int f=0 ; f<n;f++)
 {
 mas [cj][f] = mas [0][f];
 mas [0][f]= a[f];
 }
 }
 for (int i=0;i<m-1;i++)
 {
  double k = - (mas [i+n-m+1][n-m]/mas[n-m][n-m]);
 for (int j =0;j<n+1;j++)
 {
 mas [i+n-m+1][j]+=mas [i+n-m][j]*k;
 }
 }
 if (m>1)
 {
 gauss(mas,korni,n,(m-1));
 }
 for (int i=0;i < n;i++)
 {
  for (int j;j<(n);j++)
  {

   cout <<mas [i][j]<<"  ";
  }
  }
  }

  int main()
  {
  int n;
  ifstream matrix("LS1.DAT");
   matrix>>n;
   double** fad = new double* [n-1];
   for (int i=0;i < n;i++)
   {
    fad [i]= new double;
   for (int j;j<(n);j++)
   {
   matrix>> fad [i][j] ;
   cout<<fad [i][j];
   }
     }
   double* korniur = new double[n];

   gauss (fad,korniur,n,n);
   delete [] fad;
    delete [] korniur;
   return 0;
   }