% points
x = rand(50,1) ;
y = rand(50,1) ;
z = rand(50,1) ;
% vectors
u = rand(50,1)/10 ;
v = rand(50,1)/10 ;
w = rand(50,1)/10 ;
%
O = zeros(size(x)) ;
% plot
plot3(x,y,z,'.-r')
hold on
quiver3(x,y,z,u,O,O,'r')
quiver3(x,y,z,O,v,O,'b')
quiver3(x,y,z,O,O,w,'g')