import numpy as np import numpy.linalg as gauss A = np.array([[4, -2,1], [-2,4,-2], [1,-2,4]], dtype='float') #array:series / lists #Defining a matrix with unknown values b = np.array([[11,-16,17]]) #constant term gauss.solve(A,b.T) print(gauss.solve(A,b.T))