import numpy as np A = np.array([[1.,3.],[4.,2.]]) B = np.array([1.,1.]) X = np.linalg.solve(A, B) print(X)