Advertisement
Guest User

Untitled

a guest
Apr 24th, 2017
62
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 1.23 KB | None | 0 0
  1. #include<iostream>
  2. #include<fstream>
  3. #include<cstdlib>
  4.  
  5. using namespace std;
  6.  
  7. void main(){
  8.     ifstream fin("equations.txt");
  9.     ofstream fout("equations1.txt");
  10.     int n = 0;
  11.     fin >> n;
  12.     int m = n+1;
  13.     double sol[1000] = {0};
  14.  
  15.     double **matrix = (double**) malloc (m * sizeof(int*));
  16.     for(int i=0; i<m; i++){
  17.         matrix[i] = (double*)malloc(m * sizeof(int));
  18.     }
  19.  
  20.     for(int i=0; i<n; i++){
  21.         for(int j=0; j<=n; j++){
  22.             fin >> matrix[i][j];
  23.         }
  24.     }
  25.  
  26.     for (int k = 0; k < n; k++) {
  27.         for (int j = k + 1; j < n; j++) {
  28.             double d = matrix[j][k] / matrix[k][k];
  29.             for (int i = k; i < n; i++) {
  30.                 matrix[j][i] = matrix[j][i] - d * matrix[k][i];
  31.             }
  32.             matrix[j][n] -= d * matrix[k][n];
  33.         }
  34.     }
  35.  
  36.     for (int k = n - 1; k >= 0; k--) {
  37.         double d = 0.0;
  38.         for (int j = k + 1; j < n; j++) {
  39.             double s = matrix[k][j] * sol[j];
  40.             d = d + s;
  41.         }
  42.         sol[k] = (matrix[k][n] - d) / matrix[k][k];
  43.     }
  44.  
  45.     for( int i=0; i<n; i++){
  46.         fout << fixed;
  47.         fout.precision(8);
  48.         fout << sol[i];
  49.         fout << " ";
  50.     }
  51.     fin.close();
  52.     fout.close();
  53. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement