==code(cpp) |
void elim(int n,int m,double s[][]) {//sistem cu n ecuatii m necunoscute
  for(int i=1,j=1,k;i<=n && j<=m;) {

    for(k=i;k<=n; ++k) if(s[k][j]!=0) break;//cautam o linie pe care sa o folosim pentru a forma zerouri pe coloana j

    for(k=i;k<=n; ++k)
      if(s[k][j]!=0) break;//cautam o linie pe care sa o folosim pentru a forma zerouri pe coloana j

    if(k>n) {//nu am gasit nicio linie pentru care s[i][j] e nenul deci trecem la urmatoarea coloana linia i nefiind finala
      ++j;
      continue;
    }
    if(k!=i)for(int l=1; l<=m+1; ++l) swap(s[i][l],s[k][l]);//interschimbam liniile pentru a avea pe linia i si coloana j un element nenul

    for(k=i+1; k<=n; ++k) for(int l=m+1; l>=j; --l)
      s[k][l]-=((s[k][j]*s[i][l])/s[i][j]);//aplicam transformarea pentru fiecare linie mai mare ca i pentru a avea 0 pe coloana j sub linia i

    for(k=i+1; k<=n; ++k)
      for(int l=m+1; l>=j; --l)
        s[k][l]-=((s[k][j]*s[i][l])/s[i][j]);//aplicam transformarea pentru fiecare linie mai mare ca i pentru a avea 0 pe coloana j sub linia i

    ++i; ++j;
  }
//aflam necunoscutele

  for(int i=n; i;--i) for(int j=1; j<=m+1; ++j) if(fabs(s[i][j])>EPS) {//pentru ca e posibil sa avem mai multe ecuatii decat necunoscute

  for(int i=n; i;--i)
    for(int j=1; j<=m+1; ++j) if(fabs(s[i][j])>EPS) {//pentru ca e posibil sa avem mai multe ecuatii decat necunoscute

//cautam pe fiecare linie primul coeficient nenul, acestia aparand de la dreapta la stanga

    if(j==m+1) {//linia nu are coeficienti nenuli deci nu avem solutie
      g<<"Imposibil";
      exit(0);

      if(j==m+1) {//linia nu are coeficienti nenuli deci nu avem solutie
        g<<"Imposibil";
        exit(0);
      }
      x[j]=s[i][m+1];
      for(int k=j+1; k<=m; ++k) x[j]-=s[i][k]*x[k];
      x[j]/=s[i][j];
      break;//trecem la linia precedenta

    }

    x[j]=s[i][m+1];
    for(int k=j+1; k<=m; ++k) x[j]-=s[i][k]*x[k];
    x[j]/=s[i][j];
    break;//trecem la linia precedenta
  }

}
==