function b=QMULT(A,b,beta)
% b=QMULT(A,b)
% b=QMULT(A,b,beta)
%
% Multipliziert b mit Produkt von Householdermatrizen, gespeichert im unteren
% Dreieck von A


[m,n]=size(A);
