Advertisement
themlgyo

Untitled

Oct 31st, 2017
89
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.25 KB | None | 0 0
  1. void koef(double **A){
  2. double P1, Q1, W1, P2, Q2, W2,P3, Q3, W3, P4, Q4, W4;
  3. P1 = A[1][1] + (A[1][0] * (-A[0][1] / A[0][0]));
  4. Q1 = -A[1][2] / P1;
  5. W1 = (A[1][5] - A[1][0] * (A[0][5] / A[0][0])) / P1;
  6. cout << "Гамма(2): " << P1 << " Альфа(2): " << Q1 << " Бетта(2): " << W1 << endl;
  7. P2 = A[2][2] + (A[2][1] * Q1);
  8. Q2 = -A[2][3] / P2;
  9. W2 = (A[2][5] - A[2][1] * (W1) )/ P2;
  10. cout << " " << endl;
  11. cout << "Гамма(3): " << P2 << " Альфа(3): " << Q2 << " Бетта(3): " << W2 << endl;
  12. P3 = A[3][3] + (A[3][2] * Q2);
  13. Q3 = -A[3][4] / P3;
  14. W3 = (A[3][5] - A[3][2] * (W2)) / P3;
  15. cout << " " << endl;
  16. cout << "Гамма(4): " << P3 << " Альфа(4): " << Q3 << " Бетта(4): " << W3 << endl;
  17. P4 = A[4][4] + (A[4][3] * Q3);
  18. W4 = (A[4][5] - A[4][3] * (W3)) / P4;
  19. cout << " " << endl;
  20. cout << "Гамма(5): " << P4 << " Бетта(5): " << W4 << endl;
  21. double x5 = ;
  22. double x4 = Q4*x5 + W4;
  23. double x3 = Q2*x4 + W2;
  24. double x2 = Q1*x3 + W1;
  25. double x1 = (-A[1][0] / A[0][0])*x2 + A[0][5] / A[0][0];
  26. cout << "Обратный ход прогонки. Находим значения неизвестных: ";
  27. cout << "x1= " << x1 << " x2= " << x2 << " x3= " << x3 << " x4= " << x4 << " x5= " << x5 << endl;
  28. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement