Pagini recente » Cod sursa (job #2273347) | Cod sursa (job #2493225) | Cod sursa (job #773227) | Cod sursa (job #1064630) | Cod sursa (job #2743125)
#include <bits/stdc++.h>
using namespace std;
ifstream fin("gauss.in");
ofstream fout("gauss.out");
const double EPS=1e-10;
const int NMAX=305;
int N,M;
double A[NMAX][NMAX],X[NMAX];
void gauss()
{
int i = 1, j = 1, k, p;
while(i <= N && j <= M)
{
bool ok = 0;
for(k = i; k <= N; ++k)
if(abs(A[k][j]) > EPS)
{
ok = 1;
for(p = j; p <= M + 1; p++)
swap(A[i][p], A[k][p]);
break;
}
if(!ok)
{
++j;
continue;
}
for(p = j + 1; p <= M + 1; p++)
A[i][p] = A[i][p] / A[i][j];
A[i][j] = 1;
for(int u = i + 1; u <= N; u++)
{
for(p = j + 1; p <= M + 1; p++)
A[u][p] -= A[u][j] * A[i][p];
A[u][j] = 0;
}
++i;
++j;
}
}
int solutie()
{
for(int i = N; i >= 1; --i)
for(int j = i; j <= M + 1; ++j)
if(abs(A[i][j]) > EPS)
{
if(j == M + 1)
return 0;
X[j] = A[i][M + 1];
for(int k = j + 1; k <= M; ++k)
X[j] -= X[k] * A[i][k];
break;
}
return 1;
}
void afisare()
{
fout.setf(ios::fixed,ios::floatfield);
fout.precision(10);
for(int i=1;i<=M;i++)
{
fout<<X[i]<<' ';
}
}
int main()
{
fin>>N>>M;
for(int i=1;i<=N;i++)
{
for(int j=1;j<=M;j++)
{
fin>>A[i][j];
}
}
gauss();
if(solutie())afisare();
else fout<<"Imposibil";
return 0;
}