Pagini recente » Cod sursa (job #693169) | Cod sursa (job #677014) | Cod sursa (job #477099) | Cod sursa (job #819026) | Cod sursa (job #1147697)
#include <fstream>
#include <iomanip>
#define maxn 301
#define eps 1e-10
using namespace std;
ifstream fin("gauss.in");
ofstream fout("gauss.out");
int n,m,fix[maxn];
double a[maxn][maxn],b[maxn],sol[maxn];
bool null (double a)
{
if (a < 0)
a = -a;
return a < eps;
}
int main()
{
fin>>n>>m;
for (int i=1; i<=n; ++i)
{
for (int j=1; j<=m; ++j)
{
fin>>a[i][j];
}
fin>>b[i];
}
int i = 1, j = 1;
while (i <= n && j <= m)
{
int ok=0;
for (int ii=i; ii<=n; ++ii)
if ( !null (a[ii][j]) )
{
swap (a[i],a[ii]);
swap (b[i],b[ii]);
ok=1;
break;
}
if (!ok)
{
++j;
continue;
}
fix[i] = j;
for (int jj=j+1; jj<=m; ++jj)
{
a[i][jj] /= a[i][j];
}
b[i] /= a[i][j];
a[i][j] = 1;
for (int ii = i+1; ii <= n; ++ii)
{
for (int jj = j+1; jj <= m; ++jj)
{
a[ii][jj] -= (a[ii][j]*a[i][jj]);
}
b[ii] -= a[ii][j]*b[i];
a[ii][j] = 0;
}
++i; ++j;
}
for (i = n; i >= 1; --i)
{
if (!fix[i])
{
if ( !null(b[i]) )
{
fout<<"Imposibil";
return 0;
}
continue;
}
int j = fix[i];
sol[j] = b[i];
for (int ii=i-1; ii>=1; --ii)
{
b[ii] -= (a[ii][j]*sol[j]);
}
}
for (int i=1; i<=m; ++i)
{
fout<<fixed<<setprecision(10)<<sol[i]<<" ";
}
}