function [Q,R]=qrhouse1(A)
%Rastavlja matricu A na proizvod unitarne matrice Q i gornje-trougaone
%matrice R, primenom Hausholderove redukcije

[m,n]=size(A);


for i=1:n-1
    x=A([i:end],i);
    m1=length(x);
    sigma=norm(x); %norm(x,2) Euklidska norma
    k=-sigma*sign(x(1));
    beta=1/(sigma*(sigma+abs(x(1))));
    e=zeros(m1,1);
    e(1)=1;
    u=x-k*e;
    H=eye(m1)-beta*u*u';
    
    if (n-m1>0)
        T=[eye(n-m1) zeros(n-m1, m1); zeros(m1, n-m1) H];
        Q=Q*T;
    else
        T=H;
        Q=T;
    end
    A=T*A;
      
end
R=A;