I have a set of equations with variables denoted with lowercase variables and constants with uppercase variables as such
A = a + b B = c + d
To build on ewcz's answer, both the nullspace and pseudo-inverse can be calculated using numpy.linalg.svd. See the links below:
numpy.linalg.svd
pseudo-inverse
nullspace