   function [y] = hoApp(A,bet,x) 
%% function [y] = hoApp(A,bet,x) 
%% applies householder transform to 
%% a vector x
[m,n] = size(A);
y = x(:);
v = zeros(m,1);
for k=1:n
   v(k) = 1;
   v(k+1:m) = A(k+1:m,k) ;
   t  = bet(k) * v(k:m)'*y(k:m);
   y(k:m) = y(k:m) - t*v(k:m); 
end
