Tristan Desrues

Perfect !!
let you plot large dataset of coordinate, colored with velocity, in 1 line:
color_line(x,y,sqrt(vx.^2+vy.^2),'.')

instead of
for i=1:numel(x)
    plot(x(i),y(i),'.','color',myColorMap(sqrt(vx(i)^2+vy(i)^2)))
end

And it's a lot faster!!
thx!!