Yes, the very-popular NumPy package has a function to do this. Their example:
Solve the system of equations
3 * x0 + x1 = 9
andx0 + 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