Pagini recente » Cod sursa (job #3206548) | Cod sursa (job #127932) | Cod sursa (job #1288802) | Cod sursa (job #1630191) | Cod sursa (job #812984)
Cod sursa(job #812984)
#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];
float 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("%f",&a[i][j]);
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];
for(k=j;k<=m+1;k++) a[i][k]/=v;
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;
}