function u=myhouse(x)

u=x;
u(1)=sign(u(1))*norm(u)+u(1);
u=u/norm(u);