Cod sursa(job #1142214)

Utilizator Vladinho97Iordan Vlad Vladinho97 Data 13 martie 2014 16:46:17
Problema Algoritmul lui Gauss Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 0.74 kb
#include <fstream>
#include <cmath>
#define eps 1e-10
using namespace std;
double r[309],a[309][309],x[309];
int n,m;
void transf(int k,int l,double coef)
{
    int i,nr_0=0;
    double ult;
    for(i=1;i<=m;++i)
    {
        a[k][i]=a[k][i]-a[l][i]*coef;
    }
    r[k]=r[k]-r[l]*coef;
}
int main()
{
    int i,j;
    double coef;
    ifstream f("gauss.in");
    ofstream g("gauss.out");
    f>>n>>m;
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=m;++j)
            f>>a[i][j];
        f>>r[i];
    }
    for(i=n;i>0;--i)
    {
        for(j=n;j>0;--j)
        if(j!=i)
        {
            coef=a[j][i]/a[i][i];
            transf(j,i,coef);
        }
    }
    for(i=1;i<=n;i++)
        g<<r[i]/a[i][i]<<" ";

}