Vettore di Householder

v=house(x) - Calcola il vettore di Householder associato al vettore x scalato secondo la prima componente. Ovviamente se $\underline{x}=\underline{0}$ restituisce $\underline{0}$

%HOUSE
%v=house(x)
%Pre: x vettore non nullo
% Questa funzione restituisce in v il vettore che genera la
% matrice di householder
% P=I -2*(vv')/(v'v) tale che P*x=a*e1 dove a è una costante 
% reale ed e1 è il % primo vettore della base canonica.
%
% See also QRHOUSE
function v=house(x)
n=length(x);
nx=norm(x);
if nx==0
   v=zeros(n,1);
elseif x(1)>0
   temp=x(1)+nx;
   v(1)=-nx;
   v(2:n)=x(2:n)/temp;
else
   temp=x(1)-nx;
   v(1)=nx;
   v(2:n)=x(2:n)/temp;
end

return

   
   

/tmp/cn/ehouse.txt