У меня есть несколько целочисленных матриц среднего размера (несколько сотен строк). Мне нужно решить уравнения вида Ax = b, где b - стандартный базисный вектор, а A - одна из моих матриц. Я использовал numpy.linalg.lstsq для этой цели, но ошибки округления оказываются слишком значительными.
Как мне выполнить точное символьное вычисление?
(PS Мне не нужен код, чтобы быть эффективным; меня больше беспокоит легкость кодирования.)