colcomp — column compression, kernel, nullspace
[W,rk]=colcomp(A [,flag] [,tol])
real or complex matrix
character string
real number
square non-singular matrix (change of basis)
integer (rank of A)
Column compression of A: Ac = A*W is
column compressed i.e
Ac=[0,Af] with Af full column rank,
rank(Af) = rank(A) = rk.
flag and tol are optional parameters: flag = 'qr'
or 'svd' (default is 'svd').
tol = tolerance parameter (of order %eps as default value).
The ma-rk first columns of W span the kernel of A
when size(A)=(na,ma)