import numpy as np A = np.array([ [3,-9], [2,4] ]) b = np.array([-42,2]) z = np.linalg.solve(A,b) print(z) M = np.array([ [1,-2,-1], [2,2,-1], [-1,-1,2] ]) c = np.array([6,1,1]) y = np.linalg.solve(M,c) print(y)