function u=householder(a)
% gibt den Normaleneinheitsvektor u 
% der Spiegelungsebene zureck,
% mit der man a auf ein Vielfaches 
% des ersten Einheitsvektors spiegelt
norm_a=norm(a,2);
if norm_a<eps
    error('Vektor a zu klein')
end
sg=sign(a(1,1));
v=a;
v(1,1)=v(1,1)+sg*norm_a;
norm_v=norm(v,2);
u=v/norm_v;