function Q=Householder(x)
[m,n]=size(x);
tau=sign(x(1))*norm(x,2);
u=x;u(1)=u(1)+tau;
y=2/(norm(u,2)^2);
Q=eye(m,m)-y*u*u';