9. Gauss Jordan Method

 

9. Gauss Jordan Method

by ujjal roy

#include<bits/stdc++.h>

using namespace std;

int main()

{

    double matrix[3][4];

    cout<<"Enter matrix  : ";

    for(int i=0;i<3;i++)

    {

        for(int j=0;j<4;j++)

        {

            cin>>matrix[i][j];

        }

    }

    for(int i=0;i<3;i++)

    {

        for(int j=0;j<3;j++)

        {

            if(i!=j)

            {

                double d;

            d=matrix[j][i]/matrix[i][i];

            for(int k=0;k<4;k++)

            {

                matrix[j][k]-=matrix[i][k]*d;

            }

            }

        }

    }

    cout<<"X = "<<matrix[0][3]/matrix[0][0]<<endl;

    cout<<"Y = "<<matrix[1][3]/matrix[1][1]<<endl;

    cout<<"Z = "<<matrix[2][3]/matrix[2][2]<<endl;

    return 0;

}