1import numpy as np23# матрица коэффициентов4A = np.array([5 [4, -1, 1],6 [1, 5, 1],7 [9, 32, 84]8], dtype=float)910# вектор правых частей11b = np.array([12 -21,13 -12,14 -185915], dtype=float)1617# точное решение18exact = np.linalg.solve(A, b)1920print("Точное решение:")21print(exact)222324# ---------------------------25# Метод простой итерации26# ---------------------------2728x = np.zeros(3)2930print("\nМетод простой итерации")3132for i in range(3):33 x_new = np.zeros(3)3435 x_new[0] = (-21 + x[1] - x[2]) / 436 x_new[1] = (-12 - x[0] - x[2]) / 537 x_new[2] = (-1859 - 9*x[0] - 32*x[1]) / 843839 x = x_new40 print(f"Итерация {i+1}: {x}")414243# ---------------------------44# Метод Зейделя45# ---------------------------4647x = np.zeros(3)4849print("\nМетод Зейделя")5051for i in range(3):5253 x[0] = (-21 + x[1] - x[2]) / 454 x[1] = (-12 - x[0] - x[2]) / 555 x[2] = (-1859 - 9*x[0] - 32*x[1]) / 845657 print(f"Итерация {i+1}: {x}")