function out=proj(v)
out=(1/dot(v,v))*v*v';