orthogonal basis
Q=orth(A)
:A real or complex matrix : :Q real or complex matrix :
Q=orth(A) returns Q, an orthogonal basis for the span of A. Range( Q) = Range( A) and Q’*Q=eye.
The number of columns of Q is the rank of A as determined by the QR algorithm.
A=`rand`_(5,3)*`rand`_(3,4);
[X,dim]=`rowcomp`_(A);X=X';
`svd`_([orth(A),X(:,1:dim)])