Cod sursa(job #1602641)

Utilizator gorni97aaa aaa gorni97 Data 16 februarie 2016 20:57:18
Problema Algoritmul lui Gauss Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 0.85 kb
#include <iostream>
#include <fstream>
#include <iomanip>
using namespace std;
int main()

{float a[301][301],v[301];
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<=n;i++)

    g<<fixed<<v[i]<<setprecision(10)<<" ";

}