Do not use the Maple command LinearAlgebra:-ColumnSpace for this. it gives the
columns in the RREF. The correct way is to obtain the corresponding columns of
the pivot columns in the original matrix \(A\). Hence use the command Basis like
this
A:=Matrix([[1,0,0],[1,1,1]]); LinearAlgebra:-Basis([seq(A[..,i],i=1..LinearAlgebra:-ColumnDimension(A) )]);
Which gives
If you use ColumnSpace command you’ll get this
These are different. Basis is the correct command to use, which matches the standard
definition in textbooks.