Pagini recente » Cod sursa (job #65763) | Cod sursa (job #1112600) | Cod sursa (job #2341459) | Cod sursa (job #2638650) | Cod sursa (job #813023)
Cod sursa(job #813023)
#include<cstdio>
#define eps 0.00001
#define zero(x) (x>=-eps && x<=eps)
#define difzero(x) (x<-eps || x>eps)
int n,m,i,j,k,p,L[301],C[301],x;
double a[301][301],v;
int main()
{
freopen("gauss.in","r",stdin);
freopen("gauss.out","w",stdout);
scanf("%d%d",&n,&m);
for(i=1;i<=n;i++)
for(j=1;j<=m+1;j++) {scanf("%d",&x);a[i][j]=(double)x;}
for(j=1;j<=m;j++)
{
for(i=1;i<=n;i++)
if(!L[i] && difzero(a[i][j])) break;
L[i]=j;C[j]=i;
v=a[i][j];
if(v)for(k=j;k<=m+1;k++) a[i][k]/=v;
else continue;
for(k=1;k<=n;k++)
{
if(k==i || zero(a[k][j])) continue;
v=a[k][j];
for(p=j;p<=m+1;p++) a[k][p]-=a[i][p]*v;
}
}
for(i=1;i<=n;i++)
if(!L[i] && difzero(a[i][m+1]))
{
printf("Imposibil");
return 0;
}
for(j=1;j<=m;j++)
{
if(C[j]) printf("%.10f ",a[C[j]][m+1]);
else printf("0 ");
}
return 0;
}