Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #include<iostream>
- #include<fstream>
- #include<cstdlib>
- using namespace std;
- void main(){
- ifstream fin("equations.txt");
- ofstream fout("equations1.txt");
- int n = 0;
- fin >> n;
- int m = n+1;
- double sol[1000] = {0};
- double **matrix = (double**) malloc (m * sizeof(int*));
- for(int i=0; i<m; i++){
- matrix[i] = (double*)malloc(m * sizeof(int));
- }
- for(int i=0; i<n; i++){
- for(int j=0; j<=n; j++){
- fin >> matrix[i][j];
- }
- }
- for (int k = 0; k < n; k++) {
- for (int j = k + 1; j < n; j++) {
- double d = matrix[j][k] / matrix[k][k];
- for (int i = k; i < n; i++) {
- matrix[j][i] = matrix[j][i] - d * matrix[k][i];
- }
- matrix[j][n] -= d * matrix[k][n];
- }
- }
- for (int k = n - 1; k >= 0; k--) {
- double d = 0.0;
- for (int j = k + 1; j < n; j++) {
- double s = matrix[k][j] * sol[j];
- d = d + s;
- }
- sol[k] = (matrix[k][n] - d) / matrix[k][k];
- }
- for( int i=0; i<n; i++){
- fout << fixed;
- fout.precision(8);
- fout << sol[i];
- fout << " ";
- }
- fin.close();
- fout.close();
- }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement