    function [Q1,R]=qrh(A);
    % Beräkna QR-faktoriseringen av A med 
    % Householdertransformationer
    [m,n]=size(A); qt=eye(m); % Initiera q' till enhetsmatris
    R=A;
    for i=1:n
      u=househ(R(i:m,i));
      R(i:m,i:n)=apphouse(u,R(i:m,i:n));
      qt(i:m,:)=apphouse(u,qt(i:m,:));
    end
    R=triu(R(1:n,1:n));   % Högertriangulär och n x n 
    q=qt'; Q1=q(:,1:n);




