function z=HouseholderQTy(A,y);
% HOUSEHOLDERQTY applies Householder reflections transposed
%   z=HouseholderQTy(A,y); computes z=Q'y using the Householder
%   reflections Q stored as vectors in the matrix A by
%   A=HousholderQR(A)

[m,n]=size(A); z=y;
for j=1:n,
  z(j:m)=z(j:m)-A(j:m,j)*(A(j:m,j)'*z(j:m));
end;
