11.Gauss Elimination Method

11.Gauss Elimination Method

ujjal roy 

#include<bits/stdc++.h>

using namespace std;

int main()

{

   cout<<"Enter element : ";

   int i,j,n;

   n=3;

   double arr[3][4];

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

   {

       for(j=0;j<4;j++)cin>>arr[i][j];

   }

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

   {

       if(arr[i][i]==0.0)

       {

           cout<<"ERROR!!\n";

           exit(0);

       }

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

       {

           double val;

           val=arr[j][i]/arr[i][i];

           for(int k=0;k<n+1;k++)

           {

               arr[j][k]-=(arr[i][k]*val);

           }

       }

   }

double x,y,z;

z=arr[2][3]/arr[2][2];

y=(arr[1][3]-z*arr[1][2])/arr[1][1];

x=(arr[0][3]-z*arr[0][2]-y*arr[0][1])/arr[0][0];

   cout<<"X : "<<x<<endl;

   cout<<"Y : "<<y<<endl;

   cout<<"Z : "<<z<<endl;

    return 0;

}