np.linalg.solve()

Решение системы линейных уравнений A·x = b.

Сигнатураnp.linalg.solve(a, b)

Находит x в системе A @ x = b точнее и быстрее, чем через явное обращение матрицы.

import numpy as np
A = np.array([[2.0, 1.0], [1.0, 3.0]])
b = np.array([5.0, 10.0])
print(np.linalg.solve(A, b))   # [1. 3.]