Advertisement
Guest User

Untitled

a guest
Oct 21st, 2019
99
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.85 KB | None | 0 0
  1. #include "stdafx.h"
  2. #include "iostream"
  3. #include "math.h"
  4. #include "stdlib.h"
  5. #include "locale.h"
  6. #include "conio.h"
  7. using namespace std;
  8.  
  9. int main()
  10. {
  11. int i, j, n, m;
  12. setlocale(LC_ALL, "rus");
  13. //создаем массив
  14. cout << "введи число уравнений: ";
  15. cin >> n;
  16. cout << "введи число неизвестных: ";
  17. cin >> m;
  18. m+=1;
  19. float **matrix = new float *[n];
  20. for (i=0; i<n; i++)
  21. matrix[i] = new float [m];
  22.  
  23. //инициализируем
  24.  
  25. for (i = 0; i<n; i++)
  26.  
  27. for (j = 0; j<m; j++)
  28. {
  29. cout << "Элемент " << "[" << i+1 << " , " << j+1 << "]: " ;
  30.  
  31. cin >> matrix[i][j];
  32. }
  33.  
  34. //выводим массив
  35. cout << "matrix: " << endl;
  36. for (i=0; i<n; i++)
  37. {
  38. for (j=0; j<m; j++)
  39. cout << matrix[i][j] << " ";
  40. cout << endl;
  41. }
  42. cout << endl;
  43.  
  44. //Метод Гаусса
  45. //Прямой ход, приведение к верхнетреугольному виду
  46. float tmp, xx[m];
  47. int k;
  48.  
  49. for (i=0; i<n; i++)
  50. {
  51. tmp=matrix[i][i];
  52. for (j=n;j>=i;j--)
  53. matrix[i][j]/=tmp;
  54. for (j=i+1;j<n;j++)
  55. {
  56. tmp=matrix[j][i];
  57. for (k=n;k>=i;k--)
  58. matrix[j][k]-=tmp*matrix[i][k];
  59. }
  60. }
  61. /*обратный ход*/
  62. xx[n-1] = matrix[n-1][n];
  63. for (i=n-2; i>=0; i--)
  64. {
  65. xx[i] = matrix[i][n];
  66. for (j=i+1;j<n;j++) xx[i]-=matrix[i][j]*xx[j];
  67. }
  68.  
  69. //Выводим решения
  70. for (i=0; i<n; i++)
  71. cout << xx[i] << " ";
  72. cout << endl;
  73.  
  74. delete[] matrix;
  75. getch();
  76. return 0;
  77. }
  78.  
  79. вот например, 46 строка ругается
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement