function [QC,RC]=clgs(A)

 % Classical Gram-Schmidt  (Unstable) - Alg. 7.1, page 51
 %
 %

[m,n]=size(A);
 for j=1:n
   v=A(:,j);
   for i=1:j-1
      RC(i,j)=QC(:,i)'*A(:,j);
      v=v-RC(i,j)*QC(:,i);
   end
   RC(j,j)=norm(v,2);
   QC(:,j)=v/RC(j,j);
 end