figure
Um=(U1+U2)/2;
Vm=(V1+V2)/2;
figure(16)
quiver(X1,Y1,Um,Vm,'-r');