Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #include "stdafx.h"
- #include "iostream"
- #include "math.h"
- #include "stdlib.h"
- #include "locale.h"
- #include "conio.h"
- using namespace std;
- int main()
- {
- int i, j, n, m;
- setlocale(LC_ALL, "rus");
- //создаем массив
- cout << "введи число уравнений: ";
- cin >> n;
- cout << "введи число неизвестных: ";
- cin >> m;
- m+=1;
- float **matrix = new float *[n];
- for (i=0; i<n; i++)
- matrix[i] = new float [m];
- //инициализируем
- for (i = 0; i<n; i++)
- for (j = 0; j<m; j++)
- {
- cout << "Элемент " << "[" << i+1 << " , " << j+1 << "]: " ;
- cin >> matrix[i][j];
- }
- //выводим массив
- cout << "matrix: " << endl;
- for (i=0; i<n; i++)
- {
- for (j=0; j<m; j++)
- cout << matrix[i][j] << " ";
- cout << endl;
- }
- cout << endl;
- //Метод Гаусса
- //Прямой ход, приведение к верхнетреугольному виду
- float tmp, xx[m];
- int k;
- for (i=0; i<n; i++)
- {
- tmp=matrix[i][i];
- for (j=n;j>=i;j--)
- matrix[i][j]/=tmp;
- for (j=i+1;j<n;j++)
- {
- tmp=matrix[j][i];
- for (k=n;k>=i;k--)
- matrix[j][k]-=tmp*matrix[i][k];
- }
- }
- /*обратный ход*/
- xx[n-1] = matrix[n-1][n];
- for (i=n-2; i>=0; i--)
- {
- xx[i] = matrix[i][n];
- for (j=i+1;j<n;j++) xx[i]-=matrix[i][j]*xx[j];
- }
- //Выводим решения
- for (i=0; i<n; i++)
- cout << xx[i] << " ";
- cout << endl;
- delete[] matrix;
- getch();
- return 0;
- }
- вот например, 46 строка ругается
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement