    function u=househ(x)
    % Beräkna Householdervektorn u så att
    % (I - 2 u * u')x = k*e_1, där 
    % k är lika med euklidiska normen av x
    % och e_1 är första enhetsvektorn
    n=length(x);      % Antal komponenter i x
    kap=norm(x); v=zeros(n,1);
    v(1)=x(1)+sign(x(1))*kap;
    v(2:n)=x(2:n);
    u=(1/norm(v))*v;



