Pagini recente » Cod sursa (job #2790095) | Cod sursa (job #272260) | Cod sursa (job #3201711) | Cod sursa (job #1105502) | Cod sursa (job #1724532)
#include <fstream>
#include <iomanip>
using namespace std;
ifstream fin("gauss.in");
ofstream fout("gauss.out");
int n,m,i,j,k;
double a[305][305];
double sol[305];
int main()
{
fin>>n>>m;
for(i=1;i<=n;++i)
for(j=1;j<=m+1;++j)
fin>>a[i][j];
int ok=1;
for(i=1;i<=n-1;++i)
{
j=i;
while(j<=n && a[j][i]==0) ++j;
if(j==n+1) {ok=0; break;}
else if(j!=i)
{
for(k=1;k<=m+1;++k)
swap(a[i][k],a[j][k]);
}
for(j=i+1;j<=n;++j)
{
double x=-a[j][i]/a[i][i];
a[j][i]=0;
for(k=i+1;k<=m+1;++k)
a[j][k]+=a[i][k]*x;
}
}
if(!ok) {fout<<"Imposibil"; return 0;}
sol[m]=a[n][m+1]/a[n][m];
for(i=m-1;i>=1;--i)
{
double suma=0;
for(j=i+1;j<=m;++j)
suma+=sol[j]*a[i][j];
sol[i]=(a[i][m+1]-suma)/a[i][i];
}
for(i=1;i<=m;++i)
fout<<setprecision(8)<<fixed<<sol[i]<<" ";
return 0;
}