step1;
plot(1,1,'r+','linewidth',2)
u=[1/sqrt(2);1/sqrt(2)];
quiver3(1,1,0,u(1),u(1),0,'linewidth',3,'Maxheadsize',2);
