function a=uvec(v)
t=dot(v,v);
a=v./sqrt(t);