Pagini recente » Cod sursa (job #1303348) | Cod sursa (job #2073613) | Cod sursa (job #1748744) | Cod sursa (job #163315) | Cod sursa (job #1602646)
#include <iostream>
#include <fstream>
#include <iomanip>
using namespace std;
int main()
{float a[305][305],v[305];
float x;
int i,j,m,n,k,y,aux,x1;
fstream f("gauss.in",ios::in);
fstream g("gauss.out",ios::out);
f>>n>>m;
for(i=1;i<=n;i++)
for(j=1;j<=m+1;j++)
f>>a[i][j];
for(i=1;i<=n;i++)
{
if(a[i][i]==0)
{k=1;
while(a[k][i]==0)
k++;
for(j=i+1;j<=m+1;j++)
{aux=a[i][j];
a[i][j]=a[k][j];
a[k][j]=aux;}
}
for(x1=i+1;x1<=n;x1++)
for(y=i+1;y<=m+1;y++)
a[x1][y]=a[x1][y]-(a[x1][i]*a[i][y])/a[i][i];
for(j=i+1;j<=n;j++)
a[j][i]=0;
}
x=0;
for(i=n;i>=1;i--)
{
for(j=m;j>i;j--)
x=x+v[j]*a[i][j];
v[i]=(a[i][m+1]-x)/a[i][i];
x=0;
}
for(i=1;i<=m;i++)
g<<fixed<<v[i]<<setprecision(10)<<" ";
}