function [v] = house(x)
%Solving Householder matrix H
[nrows, ncols]=size(x);
e1=[1;zeros(nrows-1, ncols)];
v=sign(x(1))*norm(x)*e1 + x;
v=v/norm(v);
end