2014-12-23 12:04:50 +0200 commented answer transform list of equations to coefficientmatrix? Thank you! this is already very helpfull! although this does not work if not all equations are explicitly in all the variables. For example [x+y = 0, x+z = 5] this code does work even with coefficients that are 0 but doesn't look pretty, can someone help me make it simpler? s = list(var('x y z')) eqs = [3*x+z == -4, 2*x-3*y+5*z==2,7*x+2*y-3*z==0] lhs = map(lambda a: a.lhs(),eqs) m = [] for leftside in lhs: for i in range(3): if leftside.coeff(s[i],1 ) == 0: m.append(0) else: m.append(leftside.coeff(s[i],1 )) M=matrix(3,m) M  2014-12-20 18:17:19 +0200 asked a question transform list of equations to coefficientmatrix? I have a list of equations in multiple variables and therefor need to know if this list has a unique solution, if I had the coefficientmatrix this would be easily done by calculating the rank. So the question is, is there a fast way to transform the list of equations to a matrix?