Regarding the base field, this is something that needs to be fixed in
the SymPy matrices in general. For starters, I would just assume that
the span is complex. Once we have strong domain support in the
matrices, we can expand this.
Regarding how to implement it, I think it would be best to
I want to add a few things to sympy's matrices class such as image
and kernel of a matrix. Both of these are represented as the span of a set
of vectors.
Spans will normally take the form of span(v1,v2,v3...,vx) where all vs are
vectors of equal length.
I'm not exactly sure how to represent