In Maple there is a command "GenerateMatrix" that will take a collection of linear equations and generate a symbolic matrix A and symbolic vector b in "Ax=b" where "x" is a vector of variables specified by the command. The help page can be found here if this isn't clear:
http://www.maplesoft.com/support/help/maple/view.aspx?path=LinearAlgebra%2FGenerateMatrix
In Maxima a similar operation can be accomplished using the "augcoefmatrix" command. Is there a similar command in SymPy? How could this command be implemented if the command does not exist?