orthogonal basis
Q=orth(A)
real or complex matrix
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.