function H =householder(x,y)
%Solving Householder matrix H


rho=-sign(x(1))*norm(x)/norm(y);
y=rho*y;
v=(x-y)/norm(x-y);
%I=eye(length(x));
%H=I-2*v*v';
H=v
end