figure(2)
u=u.*f;
v=v.*f;
quiver(x,y,u,v,10)
hold on
for i=1:rows,
   for j=1:cols,
      if (f(i,j)==0),
         plot(x(i,j),y(i,j),'o')
      end
   end
end
hold off

