Solution basis of underdetermined equation set in python

Question

I have an underdetermined equation set (m equations of n variables, m smaller than n). As such, if it is solvable then the set of solutions are a linear space (if it is a homogenic set) or affine space (non-homogenic).

Is there an easy way in Python (possibly with other libraries) to obtain this space - for example, a basis of which?

Thanks.

Answer 1

Use linalg package from SciPy

Answer 2

Like the previous poster said, you'll want linalg from SciPy, but focus on the Singular Value Decomposition solution. The matrix U is the basis for the output vectors.