plot([u(1,:);v(1,:)], ...
     [u(2,:);v(2,:)], 'g-');
text(((u(1,:)+v(1,:))/2)', ((u(2,:)+v(2,:))/2)', num2str([1:size(u,2)]'));