kernel, nullspace
W=kernel(A [,tol,[,flag])
:A full real or complex matrix or real sparse matrix : :flag character string ‘svd’ (default) or ‘qr’ : :tol real number : :W full column rank matrix :
W=kernel(A) returns the kernel (nullspace) of A. If A has full column rank then an empty matrix [] is returned.
flag and tol are optional parameters: flag = ‘qr’ or ‘svd’ (default is ‘svd’).
tol = tolerance parameter (of order %eps as default value).
A=`rand`_(3,1)*`rand`_(1,3);
A*kernel(A)
A=`sparse`_(A);
`clean`_(A*kernel(A))