function [Q,R]=clGS(A)
// Datoteka
// clGS.sci
// reducirana QR faktorizacija matrice A  
// klasičnim Gram-Schmidtovim postupkom

    [m,n]=size(A);            // Dimenzije matrice A
    R=zeros(n,n);             // R=0
    Q=A;
    for j=1:n
        for i=1:j-1
           R(i,j)=Q(:,i)'*Q(:,j);
           Q(:,j) = Q(:,j) - R(i,j)*Q(:,i);
        end
        R(j,j)= norm(Q(:,j));
        Q(:,j)=Q(:,j)/R(j,j);
    end
endfunction