Solución:
Sí, el paquete NumPy muy popular tiene una función para hacer esto. Su ejemplo:
Resuelve el sistema de ecuaciones
3 * x0 + x1 = 9
yx0 + 2 * x1 = 8
:>>> import numpy as np >>> a = np.array([[3,1], [1,2]]) >>> b = np.array([9,8]) >>> x = np.linalg.solve(a, b) >>> x array([ 2., 3.])
https://docs.scipy.org/doc/numpy-1.13.0/reference/generated/numpy.linalg.solve.html
Consulte http://sympy.org/ y http://numpy.scipy.org/.
Específicamente, http://docs.scipy.org/doc/numpy/reference/routines.linalg.html
Y http://docs.sympy.org/0.7.0/tutorial.html#algebra, http://docs.sympy.org/dev/modules/solvers/solvers.html
Editar: enlace de solucionadores agregado desde el comentario.
Puede usar el método de mínimos cuadrados en Python para resolver un sistema de ecuaciones, por ejemplo, para resolver ecuaciones 3x + 4y = 7 y 5x + 6y = 8
>>> import numpy
>>> a=[[3,4],[5,6]]
>>> b=[7,8]
>>> numpy.linalg.lstsq(a,b)
(array([-5. , 5.5]), array([], dtype=float64), 2, array([ 9.27110906, 0.21572392]))