orth - calcul d'une base orthogonale
renvoie Q, une base orthogonale de l'image de A. Im(Q) = Im(A) et Q'*Q = I.
Le nombre de colonnes de Q est égal au rang de A, comme déterminé par l'algorithme QR.
A=rand(5,3)*rand(3,4); [X,dim]=rowcomp(A);X=X'; svd([orth(A),X(:,1:dim)])