#include <iostream>
#define fil 5
#define col 5
using namespace std;
int main()
{//eliminacion de gauss
double a[fil][col],b[fil],x[fil],sum,factor;
int i,j,k,n;
//lectura de a,b
//eliminacion hacia adelante
for(k-1;k<=n-1;k++)
{
for(i=k+1;i<=n;i++)
{
factor=a[k]/a[k][k];
for(j=k+1;j<=n;j++)
{
a[j]-factor*a[k][j];
}
b=b-factor*b[k];
}
}
//sustitucion hacia atras
x[n]=b[n]/a[n][n];
for(i=n-1;i>=1;i--)
{ sum=b;
for(j=i+1;j<=n;j++)
{
sum=sum-a[j]*x[j];
}
x=sum/a[j];
}
return 0;
}
#define fil 5
#define col 5
using namespace std;
int main()
{//eliminacion de gauss
double a[fil][col],b[fil],x[fil],sum,factor;
int i,j,k,n;
//lectura de a,b
//eliminacion hacia adelante
for(k-1;k<=n-1;k++)
{
for(i=k+1;i<=n;i++)
{
factor=a[k]/a[k][k];
for(j=k+1;j<=n;j++)
{
a[j]-factor*a[k][j];
}
b=b-factor*b[k];
}
}
//sustitucion hacia atras
x[n]=b[n]/a[n][n];
for(i=n-1;i>=1;i--)
{ sum=b;
for(j=i+1;j<=n;j++)
{
sum=sum-a[j]*x[j];
}
x=sum/a[j];
}
return 0;
}