fork download
copy 
  1. /*
  2. Name: Anuska Alam
  3. Id: 230241021
  4. Report: Gauss Jordan Elimination
  5. */
  6.  
  7. #include<bits/stdc++.h>
  8. using namespace std;
  9. int main (){
  10.     int n; cin >> n;
  11.     double a[101][101+1], x[100], ratio;
  12.     for (int i=1; i<=n; i++){
  13.         for (int j=1; j<=n+1; j++){
  14.             cin >> a[i][j]; 
  15.         }
  16.     }
  17.     // Convert to uper triangular form using forward elimination
  18.     for (int i=1; i<n; i++){
  19.         if (a[i][i]==0){
  20.             cout << "Error!" << endl;
  21.             break;
  22.         }
  23.         for (int j=i+1; j<=n; j++){
  24.             ratio = a[j][i]/a[i][i];
  25.             for (int k=1; k<=n+1; k++){
  26.                 a[j][k] = a[j][k] - a[i][k]*ratio;
  27.             }
  28.         }
  29.     }
  30.  
  31.     for (int i=1; i<=n; i++){
  32.         for (int j=1; j<=n+1; j++){
  33.             cout << a[i][j] << " ";
  34.         }
  35.         cout << endl;
  36.     }
  37.     // Convert to lower triangular form using backward elimination
  38.     for (int i = n; i >= 2; i--) {
  39.         for (int j = i - 1; j >= 1; j--) {
  40.             double ratio = a[j][i] / a[i][i];
  41.             for (int k = 1; k <= n + 1; k++) {
  42.                 a[j][k] = a[j][k] - ratio * a[i][k];
  43.             }
  44.         }
  45.     }
  46.  
  47.     for (int i=1; i<=n; i++){
  48.         for (int j=1; j<=n+1; j++){
  49.             cout << a[i][j] << " ";
  50.         }
  51.         cout << endl;
  52.     }
  53.  
  54.     for (int i=1; i<=n; i++){
  55.         x[i]= a[i][n+1]/a[i][i];
  56.     }
  57.     for (int i=1; i<=n; i++){
  58.         cout << x[i] << endl;
  59.     }
  60.  
  61.  
  62. }
Success #stdin #stdout 0.01s 5276KB
comments (?)
stdin
copy 
3
1 1 1 5 
2 3 5 8 
4 0 5 2
stdout
copy 
1 1 1 5 
0 1 3 -2 
0 0 13 -26 
1 0 0 3 
0 1 0 4 
0 0 13 -26 
3
4
-2
https://ideone.com/KCj9q1
language:
C++ (gcc 8.3)
created:
12 days ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!