>>> A = np.array([[2, 1], [-1, 4]])
>>> b= np.array([15, -7])
>>> np.linalg.solve(A, b)
array([7.44444444, 0.11111111])