Cod sursa(job #1777216)

Utilizator andrew_assassin789Andrei Manea andrew_assassin789 Data 12 octombrie 2016 09:51:53
Problema Algoritmul lui Gauss Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 0.96 kb
#include <fstream>
#define nmax 305
using namespace std;
double a[nmax][nmax];
double opus(double x)
{
    return -x;
}
int main()
{
    ifstream f("gauss.in");
    ofstream g("gauss.out");
    int n,m,i,j,k;
    double r,p=1;
    f>>n>>m;
    for (i=1;i<=n;i++)
    {
        for (j=1;j<=m+1;j++)
            f>>a[i][j];
    }
    for (j=1;(j<=m)&&p;j++)
    {
        r=a[j][j];
        p=p*r;
        for (k=j;(k<=m+1)&&r;k++)
        {
            a[j][k]/=r;
        }
         for (i=j+1;(i<=n)&&p;i++)
        {
            if (a[i][j]!=0)
            {
                r=-a[i][j];
                for (k=j;k<=m+1;k++)
                {
                    if (a[i][k]!=0)
                        a[i][k]+=a[j][k]*r;
                }
            }
        }
    }
    for (i=1;i<=n;i++)
    {
        for (j=1;j<=m+1;j++)
            g<<a[i][j]<<' ';
        g<<'\n';
    }
    f.close();
    g.close();
    return 0;
}