cquiver(x, y, u, v) plots quiver for at given points (x, y) with the velocity (u, v)

quiv = cquiver(x, y, u, v) returns quiver structure with two properties,
Head and Tail. Here, Head and Tail are the patches which are used for drawing heads and tails.

quiv = cquiver(x, y, u, v, 'cdata', 'angle') set color based on angle.
quiv = cquiver(x, y, u, v, 'sampling', h) plots quiver plot
at uniform grid where grid size is h.
Default sampling method is average.
Sampling method can be changed to 'fastest' by using
quiv = cquiver(x, y, u, v, 'sampling', h, 'samplingmethod', 'fastest')
In this case, the fastest one will be chosen from the rectangle
(X(i)-h/2, X(i)+h/2) x (Y(i)-h/2, Y(i)+h/2).

If you found a bug, feel free to contact me via email in the m-file.