Pagini recente » Cod sursa (job #926760) | Cod sursa (job #255927) | Cod sursa (job #3344395) | Cod sursa (job #2735626) | Cod sursa (job #3319316)
#include <fstream>
#include <iomanip>
#include <cmath>
#define INF 0.0000001
using namespace std;
ifstream fin("gauss.in");
ofstream fout("gauss.out");
int n, m, i, j;
double a[305][305], x[305];
void swap_matrice(int x, int y)
{
for(int i=1; i<=m+1; i++)
swap(a[x][i], a[y][i]);
}
int main()
{
fin>>n>>m;
for(i=1; i<=n; i++)
for(j=1; j<=m+1; j++)
fin>>a[i][j];
for(i=1; i<=n; i++)
{
for(int j=i; j<=n; j++)
if(fabs(a[j][i])>INF)
{
swap_matrice(i, j);
break;
}
if(fabs(a[i][i])>INF)
{
for(j=i+1; j<=n; j++)
{
double coef=a[j][i]/a[i][i];
for(int k=i; k<=m+1; k++)
a[j][k]-=a[i][k]*coef;
}
}
}
//daca coef mei sunt 0, dar a[..][m+1]!=0
for(int i=n; i>=1; i--)
{
bool ok=false;
for(j=1; j<=m; j++)
if(fabs(a[i][j])>INF)
{
ok=true;
break;
}
if(ok==false && fabs(a[i][m+1])>INF)
{
fout<<"Imposibil";
return 0;
}
}
for(i=n; i>=1; i--)
{
for(j=1; j<=m; j++)
{
if(fabs(a[i][j])>INF)
{
x[j]=a[i][m+1];
for(int k=j+1; k<=m; k++)
x[j]-=x[k]*a[i][k];
x[j]/=a[i][j];
break;
}
}
}
for(i=1; i<=m; i++)
fout<<fixed<<setprecision(8)<<x[i]<<" ";
return 0;
}