Not sure about efficiency, but perhaps simpler:
P = [p p p]
V = [u v w]
C = repeat([:red,:green,:blue], inner=4)
quiver(P[1,:], P[2,:], quiver=(V[1,:],V[2,:]), c=C)
NB: cannot get rid of warning message related to NaNs inserted by quiver()
Not sure about efficiency, but perhaps simpler:
P = [p p p]
V = [u v w]
C = repeat([:red,:green,:blue], inner=4)
quiver(P[1,:], P[2,:], quiver=(V[1,:],V[2,:]), c=C)
NB: cannot get rid of warning message related to NaNs inserted by quiver()