import numpy as np a = np.array([[6,2,-5], [3,3,-2], [7,5,-3]]) b = np.array([13,13,26]) x = np.linalg.solve(a, b) print(x)