    function Y=apphouse(u,X);
    % Multiplicera med Householdermatris
    % Y = (I - 2 * u * u') * X
    Y=X-2*u*(u'*X);

