Gauss Elimination

 #include <iostream>

#include <iomanip>

using namespace std;

int main()

{

    int i,j,k,n;

    float a[10][10],x[10];

    cout<<"Enter order of equation: ";

    cin>>n;

    cout<<"Enter coefficient of equation row-wise\n";

    for(i = 0;i<n;i++)    //i = 1;i<=n;i++

    {

        for(j=0;j<=n;j++) //j = 1

        {

            cin>>a[i][j];

        }

    }

    //miain eqn

    for(i=0;i<n;i++) //i=1;i<=n

    {

        for(j=i+1;j<=n;j++) //j=1

        {

            if(j>i)

            {

                float s = a[j][i]/a[i][i];

                for(k=i;k<=n;k++) //k=1;k<=n+1

                {

                    a[j][k] = a[j][k] - s*a[i][k];

                }

            }

        }

    }

    setprecision(10);

    x[n-1] = a[n-1][n]/a[n-1][ n-1]; //x[n] = a[n][n+1]/a[n][n];

    for(i=n-2;i>=0;i--) //i=n-1;i>=1

    {

        float sum = 0.0;

        for(j=i+1;j<n;j++)

        {

            sum = sum + a[i][j] * x[j];

        }

        x[i] = (a[i][n]-sum)/a[i][i]; //x[i] = (a[i][n+1]-sum)/a[i][i];

    }

    cout<<"\nThe solutions are:\n\n";

    for(i =0;i<n;i++) //i=1;i<=n

    {

        cout<<"X"<<i<<" = "<<x[i]<<"\n";

    }

    return 0;

    

}