function [Q,R]=mGS(A)
// Datoteka:
// mGS.sci
// reducirana QR faktorizacija matrice A  
// stabiliziranim Gram-Schmidtovim postupkom

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

endfunction